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INTRODUCTION 

Much  research  has  examined  the  potential  for  axon  regeneration  across  an  injury  site  to  re¬ 
establish  functional  connectivity  between  higher  motor  centers  and  spinal  networks12. 
Functional  recovery  may  remain  elusive,  however,  should  the  networks  that  carry  out  movement 
be  so  altered  after  injury  that  they  no  longer  respond  to  a  newly  reconnected  source  of  input. 
The  overall  goal  of  this  project  was  to  address  an  understudied  aspect  of  regeneration  -  how 
does  removal  of  descending  neural  input  alter  motor  network  activity  below  the  site  of  injury, 
and  what  are  its  consequences  should  we  become  able  to  reconnect  higher  motor  centers  to 
motor  networks?  The  study  objective  was  to  begin  to  understand  the  extent  and  mechanisms  by 
which  neural  networks  change  after  removing  descending  input,  and  to  characterize  activity  that 
results  from  re-establishing  connectivity.  Our  over-arching  hypothesis  was  that  neural  networks 
are  fundamentally  altered  when  descending  inputs  are  removed,  therefore  making  the  effects  of 
reconnecting  these  networks  unpredictable.  Ultimately  the  goal  is  to  understand  the  full  suite  of 
medical  intervention  that  may  be  necessary  to  allow  for  spinal  cord  injury  recovery,  including 
preventing  or  rectifying  changes  in  neural  networks  below  the  injury  site. 

BODY 

From  the  point  of  view  of  the  nervous  system,  the  cure  for  spinal  cord  injury  must  progress  on  at 
least  two  complementary  fronts:  re-establishing  connectivity  between  the  brain  and  spinal 
networks  (regeneration)2,  and  understanding  how  to  safeguard  or  recover  functional  viability  of 
isolated  spinal  networks  in  order  to  receive  input  from  regenerated  tracts3"5.  Our  work  is  aimed 
at  elucidating  the  biological  basis  by  which  isolated  spinal  networks  may  remain  functionally 
viable,  which  is  the  foundation  for  developments  that  aim  to  increase  function  and  recovery 
among  those  who  have  suffered  a  spinal  cord  injury.  This  basic  research  approach  to  uncover 
mechanism  of  how  networks  are  altered  by  denervation  (disrupting  nerve  input)  is  a  critical  first 
step  towards  this  goal,  as  simply  describing  these  phenomena  will  fail  to  add  insight  into 
possible  future  therapeutic  approaches.  Only  by  understanding  what  is  changing,  as  opposed 
to  simply  understanding  that  networks  change,  will  we  progress  towards  viable  treatment 
options. 

The  CDMRP  was  will 

neural  networks  in  an  unlikely  system  for  spinal  cord  injury:  namely  crustacean  motor  networks 
such  as  the  stomatogastric  ganglion  (STG)6.  While  the  STG  is  an  invertebrate  preparation,  it 
provides  many  parallels  to  the  organization  of  motor  networks  in  mammalian  systems,  yet  offers 
a  substantial  level  of  manipulation  that  is  currently  unachievable  in  mammalian  systems.  Our 
overall  model  system  for  this  proposal  was  the  transient  removal  and  restoration  of  descending 

subsequent  function  with  restored  inputs.  In  addition,  we  focus  on  the  underlying  cellular 
mechanisms  involved  in  these  kinds  of  changes, 

disruption  in  that  activity.  Our  hypothesis  is  that  the  onset  of  these  homeostatic  mechanisms  is 
actually  detrimental  to  the  restoration  of  activity  through  restored  connectivity  of  descending 
inputs. 
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We  made  substantial  progress  with  this  proposal.  Here  we  summarize  our  progress  in 
understanding  the  functional  impacts  of  removal  of  descending  inputs  to  motor  networks.  For 
clarity  and  succinctness,  I  will  report  these  results  in  a  bulleted  list  format.  All  of  the  following 
descriptions  of  results  from  this  proposal  can  be  supported  by  data  in  the  attached  publications 
found  in  the  Appendices. 

THE  FOLLOWING  SUMMARY  IS  PUBLISHED  AND  FOUND  IN  THE  FOLLOWING  ARTICLE 
ATTACHED  IN  THE  APPENDICES: 

Nahar  J,  Lett  KM,  Schulz  DJ  (2012)  Restoration  of  descending  inputs  fails  to  rescue  activity  following 
deafferentation  of  a  motor  network.  Journal  of  Neurophysiology  108:  871-881. 

•  Removal  of  modulatory,  descending  projections  (stomatogastric  nerve  STN)  to  the 
motor  network  (the  STG)  causes  an  immediate  and  complete  loss  of  the  motor  output. 

blocks  all  descending  neural  activity  without  injuring  the  axons  of  the  blockade.  The 
block  is  completely  reversible  and  under  control  of  the  experimenter  at  all  times,  so  we 
can  decentralize  for  any  time  period  we  like  and  then  reconnect  the  descending  control 
centers  to  the  network  at  will. 

•  Over  the  first  ~1 2  hours  of  blockade,  subsequent  restoration  of  descending  inputs  is  able 
to  re-establish  full  functional  motor  output.  Subsequent  periods  of  48h,  72h,  and  96h  of 
blockade  lead  to  an  inability  to  restore  full  functionality  in  the  motor  network,  with  the 
effect  becoming  more  pronounced  over  longer  blockade  schemes.  We  performed 
multiple  experiments  to  ensure  the  health  of  the  blocked  axons  in  these  experiments, 
and  determined  that  sucrose  blockade  is  not  damaging  these  fibers.  This  effect  seems 
to  be  some  physiological  change  in  the  system  that  alters  the  ability  of  the  descending 
projections  and  motor  network  to  communicate  effectively  after  extended  blockade. 

•  The  effects  above  are  seemingly  not  due  to  a  loss  of  the  descending  projections  to 
release  neuromodulator.  After  3  days  of  blockade,  we  can  electrically  stimulate  the 
descending  nerve  both  above  (with  the  block  removed)  and  below  the  block  (with  the 
block  intact)  and  see  an  increase  in  motor  activity  of  equal  magnitude.  But  neither  of 
these  stimulation  regimes  is  able  to  restore  full  functionality  to  the  motor  network. 

•  There  is  a  significant  decrease  in  the  activity  of  the  descending  projections  over  the 
course  of  the  blockade  that  is  highly  correlated  to  the  decrease  in  the  ability  of  restored 
projections  to  drive  appropriate  motor  output.  This  was  perhaps  our  least  expected 
result.  We  had  initiated  these  experiments  to  look  at  changes  in  the  motor  network  itself 
as  a  result  of  it  losing  inputs,  but  apparently  altering  feedback  interactions  between 
higher  centers  and  the  motor  network  alter  activity  of  the  descending  projections 
themselves. 

•  There  is  a  significant  change  in  the  sensitivity  of  the  STG  to  exogenously  applied 
neuromodulation  following  3  days  of  blockade.  However,  this  is  not  across  all 
modulators.  A  neuropeptide  (proctolin)  causes  enhanced  responses  in  blocked 
preparations,  but  responses  to  a  biogenic  amine  analog  of  noradrenaline  (octopamine) 
show  no  differences  between  blocked  and  unblocked  preparations. 
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THE  FOLLOWING  SUMMARY  OF  RESEARCH  IS  PUBLISHED  AND  FOUND  IN  THE 
FOLLOWING  ARTICLE  ATTACHED  IN  THE  APPENDICES: 

Ransdell  JL,  Nair  SS,  Schulz  DJ  (2012)  Rapid  homeostatic  plasticity  of  intrinsic  excitability  in  a  central 
pattern  generator  network  stabilizes  functional  neural  network  output.  Journal  of  Neuroscience  32:  9649- 
9658. 


•  We  determined  that  mRNA  levels  in  single  identified  motor  neurons  for  the  channels 
which  encode  U  and  kca  are  positively  correlated,  yet  the  ionic  currents  themselves  are 
negatively  correlated,  across  a  population  of  motor  neurons.  We  then  determined  that 
these  currents  are  functionally  coupled;  decreasing  levels  of  either  current  within  a 
neuron  causes  a  rapid  increase  in  the  other.  These  results  suggest  that  compensatory 
interactions  may  exist  that  are  utilized  in  times  of  perturbation  to  cell  output  that  could 
ultimately  be  initiated  by,  and  influence  recovery  from,  injury. 

•  This  functional  interdependence  results  in  homeostatic  stabilization  of  both  the  individual 
neuronal  and  the  network  output.  Perturbations  that  alter  ionic  currents  and  disrupt 

create  behaviorally  appropriate  output.  However,  after  a  relatively  short  time  (within  one 
hour)  there  is  stabilizing  compensation  that  occurs  to  help  restore  normal  network 
activity.  This  restoration  causes  a  significant  change  in  the  underlying  physiology  of  the 
neuron,  resulting  in  a  compensatory  increase  in  lA  when  lKca  is  blocked,  and  vice  versa. 

•  These  compensatory  increases  are  mechanistically  independent  (proceeding  by  distinct 
intracellular  signaling  mechanisms),  suggesting  robustness  in  the  maintenance  of  neural 
network  output  that  is  critical  for  survival.  Taken  together,  we  generate  a  complete 
model  for  homeostatic  plasticity  from  mRNA  to  network  output  where  rapid  post- 
translational  compensatory  mechanisms  acting  on  a  reservoir  of  channels  proteins 
regulated  at  the  level  of  gene  expression  provide  homeostatic  stabilization  of  both 
cellular  and  network  activity.  We  feel  these  compensatory  mechanisms,  while  initially 
beneficial,  may  ultimately  impact  long-term  recovery  from  perturbation  or  cause  tradeoffs 
in  the  n 

THE  FOLLOWING  SUMMARY  OF  RESEARCH  IS  PUBLISHED  AND  FOUND  IN  THE 
FOLLOWING  ARTICLE  ATTACHED  IN  THE  APPENDICES: 

Temporal,  S.,  et  al.,  (2012)  Neuromodulation  independently  determines  correlated  channel  expression  and 
conductance  levels  in  motor  neurons  of  the  stomatogastric  ganglion.  J  Neurophysiology  107:  718-727. 

•  Neurons  that  undergo  loss  of  modulatory  inputs  do  not  show  gross  overall  changes  in 
expression  level  for  ion  channels.  However,  this  does  not  indicate  that  gene  expression 
for  these  critical  components  is  not  changing  (see  next  point). 

•  Neurons  that  undergo  a  loss  of  connectivity  show  a  change  in  the  relationship  of  their 
expressed  ionic  membrane  currents.  Loss  of  neuromodulation  results  in  changes  in  the 
ratios  of  Ga,  Ghtk  (high  threshold  K+),  and  Gh,  which  leads  to  changes  in  output  (or  lack 
of  appropriate  output).  In  part  these  relationships  remain  intact  with  neuromodulation, 
despite  activity-dependent  manipulations,  suggesting  in  part  that  neuromodulation 
directly  influences  ionic  current  level  relationships. 
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•  Neurons  that  undergo  a  loss  of  connectivity  show  changes  in  their  underlying  ion 
channel  gene  expression  patterns.  For  example,  under  control  conditions,  levels  of 
mRNA  levels  for  ion  channels  show  correlated  patterns  of  expression.  In  other  words, 
there  are  correlations  in  mRNA  copy  number  between  channel  genes  for  voltage-gated 
channels  such  as  shal,  BKKCa,  and  IH.  Following  decentralization,  these  correlated 
patters  of  conductance  and  channel  expression  are  substantially  disrupted. 

THE  FOLLOWING  SUMMARY  OF  RESEARCH  IS  PUBLISHED  AND  FOUND  IN  THE 
FOLLOWING  ARTICLE  ATTACHED  IN  THE  APPENDICES: 

Ransdell  JL,  Temporal  S,  West  NL,  Leyrer  ML,  Schulz  DJ  (2013)  Characterization  of  inward  currents 
and  channels  underlying  burst  activity  in  motor  neurons  of  the  crab  cardiac  ganglion.  Journal  of 
Neurophysiology  110:42-54. 

•  As  a  result  of  the  above  work  on  ionic  currents,  and  as  a  companion  to  some  of  the 
forthcoming  work  that  has  not  yet  been  published  (see  below),  we  were  able  to 
determine  that  voltage-gated  calcium  currents  in  particular  were  likely  to  be  significantly 
involved  in  the  loss  of  function  following  loss  of  neuromodulation.  Yet  not  thorough 
characterization  of  these  calcium  currents,  or  the  channels  that  carry  them,  had  ever 
been  carried  out  in  our  model  system.  This  became  an  impediment  to  our  further 
understanding  of  the  impacts  of  injury  in  this  system. 

•  We  were  able  to  carry  out  a  full  and  thorough  characterization  of  these  channels  and 
currents  in  our  system  to  move  this  research  forward.  We  determined  that  in  particular 
two  voltage-gated  calcium  currents  are  responsible  for  the  excitability  of  the  cells  in  our 
model  system.  Furthermore,  we  identified  a  persistent  sodium  current  and  a  calcium- 
activated  inward  current  that  were  previously  uncharacterized. 

•  We  performed  a  molecular  characterization  of  the  channels  that  also  carry  these 
currents  so  that  we  could  employ  them  in  our  molecular  analyses  as  described 
throughout  the  proposal  and  this  report  (see  Figs  2  and  3  below).  Because  these 
currents  are  the  major  driving  force  for  activity  in  these  cells,  this  became  a  crucial 
addition  to  our  understanding  of  the  system. 

THE  FOLLOWING  SUMMARY  OF  RESEARCH  IS  PUBLISHED  AND  FOUND  IN  THE 
FOLLOWING  ARTICLE  ATTACHED  IN  THE  APPENDICES: 

Ransdell  JL,  Nair  SS,  Schulz  DJ  (2013)  Neurons  within  the  same  network  independently  achieve 
conserved  output  by  differentially  balancing  variable  conductance  magnitudes.  Journal  of  Neuroscience 
33:  9950-9956. 

•  In  the  process  of  this  work,  we  also  identified  a  critical  and  novel  phenomenon  regarding 
the  variability  of  ionic  currents  and  channel  expression  that  underlie  activity  in  these 
motor  networks.  Specifically,  motor  neurons  with  similar  outputs  have  significantly 
different  underlying  physiology  and  molecular  expression  profiles.  The  implications  of 
these  findings  are  that  injury  may  differentially  affect  motor  neurons  in  networks,  and 
therefore  the  potential  application  of  therapeutics  may  be  unpredicatable. 
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•  For  the  first  time,  we  definitively  showed  that  compensation  (see  Ransdell  et  al.  2012) 
and  variability  result  in  fundamentally  different  mechanisms  for  generating  convergent 
output.  When  combined  with  our  overall  results  that  injury  and  loss  of  input  results  in 
changing  of  ionic  currents,  and  channel  expression,  this  suggests  that  an  injury 
response  will  cause  disparate  effects  across  individuals.  This  has  substantial 
implications  for  therapy  and  our  understanding  of  injury  effects  going  forward. 


THE  FOLLOWING  SUMMARY  OF  RESEARCH  IS  NOT  YET  PUBLISHED  BUT  WILL  BE 
PUBLISHED  AS  COMPLETE  DATA  ANALYSIS  AND  WRITING  MOVES  FORWARD. 
FIGURES  PROVIDED  AS  SUPPORTING  DATA  AND  INCLUDED  IN  THE  BODY  OF  THIS 
REPORT  AS  CITED: 

•  Neurons  that  undergo  loss  of  modulatory  inputs  do  not  show  gross  overall  changes  in 
expression  level  for  ion  channels  (Fig  1).  However,  this  does  not  indicate  that  gene 
expression  for  these  critical  components  is  not  changing  (see  next  point). 

•  Neurons  that  undergo  a  loss  of  connectivity  show  changes  in  their  underlying  ion 
channel  gene  expression  patterns.  For  example,  under  control  conditions,  levels  of 
mRNA  levels  for  ion  channels  show  correlated  patterns  of  expression  (Fig  2).  In  other 
words,  there  are  correlations  in  mRNA  copy  number  between  channel  genes  for  voltage¬ 
gated  channels  such  as  shal  (Kv4.1),  shab  (Kv2.1),  Ca(v)1  (Cav1.1),  Ca(v)2( Cav1.2),  and 
Ca(v)3  (Cav3.1).  Following  decentralization,  these  correlated  patters  of  conductance 
and  channel  expression  are  completely  disrupted  (Fig  2). 

•  This  loss  of  correlated  channel  expression  is  in  part  due  to  the  loss  of  neurochemical 
modulation  directly,  and  not  just  the  activity  that  elicited  by  them.  If  we  disrupt  activity, 

correlated  patterns  of  channel  expression  (Fig  3). 


Taken  together,  our  work  has  set  the  stage  for  an  improved  understanding  and  investigation 
of  the  underlying  mechanisms  of  the  changes  in  motor  networks  that  are  deprived  of  their 
descending  inputs.  Overall,  the  changes  are  widespread,  seemingly  affecting  cells  both  above 
and  below  the  injury/block  site.  Also  quite  striking  is  the  fact  that,  over  the  first  set  of  time 
scales  we  examined,  that  the  effect  was  not  reversible  simply  by  restoring  the  connections. 


KEY  RESEARCH  ACCOMPLISHMENTS 


We  have  established  that  motor  networks  do  not  respond  appropriately  to  reconnected 
descending  inputs  within  24  hours  of  the  initial  injury  or  blockade 
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•  Descending  projections  are  still  capable  of  releasing  neurotransmitters/neuromodulators, 
suggesting  that  transient  block  does  not  cause  a  loss  of  ability  to  release  chemical 
transmitters. 

•  Targets  (motor  neurons)  of  descending  projections  have  altered  sensitivity  to 
neuromodulatory  compounds. 

•  Targets  (motor  neurons)  undergo  changes  to  their  intrinsic  properties  and  gene 
expression,  including  loss  of  coordinated  expression  of  ionic  conductances  and  mRNA 
levels  for  the  channels  underlying  those  conductances. 

•  The  activity  of  descending  inputs  is  also  seemingly  altered  by  this  blockade,  presumably 
as  a  result  of  loss  of  feedback  interactions  with  the  motor  network  itself. 

•  Restoring  inputs  fails  to  restore  both  activity  and  the  underlying  physiology  over  a  time 
course  of  multiple  days.  That  is,  the  effects  listed  above  are  not  ameliorated  by 
reconnecting  the  inputs,  at  least  over  the  time  scale  we  have  investigated. 

•  Functional  interdependence  of  ionic  currents  results  in  homeostatic  stabilization  of  both 
the  individual  neuronal  and  the  network  output.  Perturbations  that  alter  ionic  currents 

ability  to  create  behaviorally  appropriate  output. 

•  After  a  relatively  short  time  (within  one  hour)  there  is  stabilizing  compensation  that 
occurs  to  help  restore  normal  network  activity.  This  restoration  causes  a  significant 
change  in  the  underlying  physiology  of  the  neuron,  resulting  in  a  compensatory  increase 
in  lA  when  lKca  is  blocked,  and  vice  versa. 

•  Neurons  that  lose  descending  modulatory  input  show  altered  patterns  of  gene 
expression  that  correlate  with  disrupted  activity.  However,  these  patterns  are  only  seen 
by  a  more  thorough  analysis  of  relationships  among  gene  products,  as  mean  amounts 
are  not  seen  to  significantly  change  across  individuals. 

•  These  changes  in  gene  expression  patterns  can  in  part  be  prevented  by  exogenous 
application  of  neuromodulators  that  mimic  that  actions  of  endogenous  neuromodulatory 
compounds. 

•  Further  analysis  shows  that  changes  in  gene  expression  are  likely  a  result  of  loss  of 
activity  of  these  neurons  through  an  activity  dependent  feedback  mechanism. 

•  Changes  in  these  physiological  parameters  as  a  result  of  injury  are  likely  to  be  variable 
and  unpredictable  across  individuals  as  a  result  of  natural  inherent  variability  in  the 
properties  of  their  motor  neurons. 
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REPORTABLE  OUTCOMES 

The  work  funded  by  CDMRP  has  resulted  in  5  publications  thus  far: 

1.  Temporal,  S.,  et  al.,  (2012)  Neuromodulation  independently  determines  correlated 

channel  expression  and  conductance  levels  in  motor  neurons  of  the  stomatogastric 
ganglion.  J  Neurophysiology  -727.  See  Appendix  1 . 

2.  Ransdell  JL,  Nair  SS,  Schulz  DJ  (2012)  Rapid  homeostatic  plasticity  of  intrinsic 
excitability  in  a  central  pattern  generator  network  stabilizes  functional  neural  network 
output.  Journal  of  Neuroscience  32:  9649-9658.  See  Appendix  2. 

3.  Nahar  J,  Lett  KM,  Schulz  DJ  (2012)  Restoration  of  descending  inputs  fails  to  rescue 
activity  following  deafferentation  of  a  motor  network.  Journal  of  Neurophysiology  1 08: 
871-881.  See  Appendix  3. 

4.  Ransdell  JL,  Temporal  S,  West  NL,  Leyrer  ML,  Schulz  DJ  (2013)  Characterization  of 
inward  currents  and  channels  underlying  burst  activity  in  motor  neurons  of  the  crab 
cardiac  ganglion.  Journal  of  Neurophysiology  1 1 0:42-54.  See  Appendix  4. 

5.  Ransdell  JL,  Nair  SS,  Schulz  DJ  (2013)  Neurons  within  the  same  network  independently 
achieve  conserved  output  by  differentially  balancing  variable  conductance  magnitudes. 
Journal  of  Neuroscience  33:  9950-9956.  See  Appendix  5. 


In  addition,  we  have  presented  this  work  as  a  total  of  9  contributed  or  invited  presentations: 

1.  Temporal  S,  Lett  KM,  Schulz  DJ  (2011)  Recovery  of  rhythmic  activity  in  decentralized 
cells  of  the  stomatogastric  ganglion  is  independent  of  correlated  levels  of  channel 
expression.  Society  for  Neuroscience  Annual  Meeting  201 1 ,  Washington  DC. 

2.  Lett  KM,  Nahar  J,  Lane  B,  Schulz  DJ  (2011)  Decentralization  alters  neuromodulator 
sensitivity  in  the  pyloric  network  of  the  stomatogastric  ganglion.  Society  for 
Neuroscience  Annual  Meeting  201 1 ,  Washington  DC. 

3.  Ransdell  J,  Schulz  DJ.  Generation  and  Maintenance  of  Motor  Neuron  Bursts  and 

Network  Behavior  by  an  IA  -  IHTK  Compensatory  Relationship.  Society  for 

Neuroscience,  Washington  DC,  November  201 1 . 

4.  S  Individual 

and  Species  Variation  in  Neuroendocrine  Mechanisms  th  annual  Society  for 

Behavioral  Endocrinology  meeting,  Queretaro,  Mexico. 

5.  Lett  KM,  Nahar  J,  Schulz  DJ  (2012)  Isolated  pattern  generating  network  fails  to  respond 
to  restored  modulatory  input.  Society  for  Neuroscience  Annual  Meeting  2012,  New 
Orleans,  LA. 

6.  Temporal  S,  Schulz  DJ  (2012)  Neuromodulation  maintains  channel  mRNA  correlations 
in  neurons  of  the  stomatogastric  ganglion.  Society  for  Neuroscience  Annual  Meeting 
2012,  New  Orleans,  LA. 
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7.  Ransdell  J,  Nair  SS,  Schulz  DJ  (2012)  Rapid  homeostatic  plasticity  of  intrinsic  excitability 
and  electrical  synapse  strength  in  a  central  pattern  generator  network  stabilizes  neuron 
and  network  output.  Society  for  Neuroscience  Annual  Meeting  2012,  New  Orleans,  LA. 

8.  Lane  B,  Lett  K,  Schulz  DJ  (2012)  The  Effect  of  Neuromodulation  on  the  Synchronization 
of  Motor  Neurons  in  the  Cardiac  Ganglion  of  C.  borealis.  Satellite  Meeting  -  Society  for 
Neuroscience  Annual  Meeting  2012,  New  Orleans,  LA. 

9.  Ransdell  J,  Nair  SS,  Schulz  DJ.  Rapid  homeostatic  plasticity  of  intrinsic  excitability  and 
electrical  synapse  strength  in  a  central  pattern  generator  network  stabilizes  neuron  and 
network  output.  Society  for  Neuroscience,  New  Orleans,  LA,  October  2012. 

The  data  generated  by  this  research  has  allowed  me  to  apply  for  other  sources  of  funding  to 
continue  this  work,  as  well  as  transition  the  hypotheses  being  developed  into  higher 
mammalian  systems.  While  a  subsequent  DoD  CDMRP  proposal  based  on  these  findings 
was  declined  for  further  support,  the  Missouri  Spinal  Cord  Injuries  Research  Program  has 
funded  new  mammalian  work  in  my  laboratory  that  builds  off  of  some  of  the  hypotheses 
generated  by  the  research  in  this  report. 
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CONCLUSION 

Taken  together,  our  work  has  set  the  stage  for  an  improved  understanding  and  investigation 
of  the  underlying  mechanisms  of  the  changes  in  motor  networks  that  are  deprived  of  their 
descending  inputs.  Overall,  the  changes  are  widespread,  seemingly  affecting  cells  both  above 
and  below  the  injury/block  site.  Also  quite  striking  is  the  fact  that,  over  the  first  set  of  time 
scales  we  examined,  that  the  effect  was  not  reversible  simply  by  restoring  the  connections.  Our 
work  will  continue  to  characterize  the  nature  of  changes  both  above  and  below  block/injury 
sites,  as  well  as  determine  the  time  course  over  which  these  changes  occur.  We  will  attempt  to 
determine  whether  the  effect  is  reversible  with  restored  connectivity  over  time,  and  also 
determine  whether  intervention  can  occur  during  the  blockade  time  to  prevent  the  changes  that 
are  causing  loss  of  motor  function.  Ultimately  the  goal  is  to  determine  whether  these  results  can 
be  applied  to  spinal  cord  injury  to  begin  to  inform  the  full  suite  of  therapeutic  intervention 
necessary  to  carry  out  regeneration  of  injured  fibers  (not  in  the  scope  of  this  proposal),  but  also 
to  ensure  that  reconnected  networks  appropriately  communicate. 
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Figure  1.  Comparison  of  mean  mRNA  (+/-  SD)  expression  levels  for  different  channels  of  the  STG- 
pacemaker  cell,  PD  under  various  network  conditions  that  had  either  neuromodulation,  activity, 
neither  or  both  inputs  removed. 
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Figure  2.  A  Correlations  among  ion  channel  expression  level  from  intact  control  and  decentralized 
conditions  for  PD,  the  STG  pacemaker  cell.  Open  circles  represent  the  mRNA  levels  of  an  individual  cell. 
Regression  lines  are  shown  for  significant  correlations  (Pearson’s  test  P<0.05).  B  Representative  Ivn 
extracellular  recording  to  show  typical  pyloric  rhythm  prior  to  pulling  individual  cell. 
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Figure  3.  Correlations  among  ion  channel  expression  levels  of  the  STG  pacemaker  cell,  PD  under  various 
conditions.  Open  circles  represent  the  mRNA  levels  of  an  individual  cell.  Regression  lines  are  shown  for 
significant  correlations  (Pearson’s  test  P<0.05).  Under  each  group  of  channel  correlations  is  an 
extracellular  recording  of  the  Ivn  to  show  the  typical  pyloric  activity  under  each  condition.  To  determine 
factors  that  PD  channel  correlations  may  be  dependent  on,  neuromodulation  and  activity  were 
decoupled.  Time-control  are  intact  networks  incubated  for  3  days  to  take  into  account  biological  changes 
due  to  long  incubation  periods.  Recovered  are  decentralized  networks  that  regained  the  pyloric  network 
function  following  several  days  of  incubation.  Desynapsed  have  neuromodulatory  input  but  activity  was 
disrupted  by  blocking  glutamatergic  synapses  using  picrotoxin  (PTX).  Pilocarpine  artificially  maintains 
activity  through  one  modulator,  the  muscarinic  agonist  pilocarpine,  not  the  normal  descending  inputs. 
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correlated  channel  expression  and  conductance  levels  in  motor  neurons  of  the  stomatogastric 
ganglion.  J  Neurophysiology  1 07:  71 8-727 
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ronal  identity  depends  on  the  regulated  expression  of  numerous 
molecular  components,  especially  ionic  channels,  which  determine 
the  electrical  signature  of  a  neuron.  Such  regulation  depends  on  at 
least  two  key  factors,  activity  itself  and  neuromodulatory  input. 
Neuronal  electrical  activity  can  modify  the  expression  of  ionic  cur¬ 
rents  in  homeostatic  or  nonhomeostatic  fashion.  Neuromodulators 
typically  modify  activity  by  regulating  the  properties  or  expression 
levels  of  subsets  of  ionic  channels.  In  the  stomatogastric  system  of 
crustaceans,  both  types  of  regulation  have  been  demonstrated.  Fur¬ 
thermore,  the  regulation  of  the  coordinated  expression  of  ionic  cur¬ 
rents  and  the  channels  that  carry  these  currents  has  been  recently 
reported  in  diverse  neuronal  systems,  with  neuromodulators  not  only 
controlling  the  absolute  levels  of  ionic  current  expression  but  also, 
over  long  periods  of  time,  appearing  to  modify  their  correlated 
expression.  We  hypothesize  that  neuromodulators  may  regulate  the 
correlated  expression  of  ion  channels  at  multiple  levels  and  in  a 
cell-type-dependent  fashion.  We  report  that  in  two  identified  neuronal 
types,  three  ionic  currents  are  linearly  correlated  in  a  pairwise  manner, 
suggesting  their  coexpression  or  direct  interactions,  under  normal 
neuromodulatory  conditions.  In  each  cell,  some  currents  remain  cor¬ 
related  after  neuromodulatory  input  is  removed,  whereas  the  correla¬ 
tions  between  the  other  pairs  are  either  lost  or  altered.  Interestingly,  in 
each  cell,  a  different  suite  of  currents  change  their  correlation.  At  the 
transcript  level  we  observe  distinct  alterations  in  correlations  between 
channel  mRNA  amounts,  including  one  of  the  cell  types  lacking  a 
correlation  under  normal  neuromodulatory  conditions  and  then  gain¬ 
ing  the  correlation  when  neuromodulators  are  removed.  Synaptic 
activity  does  not  appear  to  contribute,  with  one  possible  exception,  to 
the  correlated  expression  of  either  ionic  currents  or  of  the  transcripts 
that  code  for  the  respective  channels.  We  conclude  that  neuromodu¬ 
lators  regulate  the  correlated  expression  of  ion  channels  at  both  the 
transcript  and  the  protein  levels. 

voltage  clamp;  ion  channel;  mRNA;  correlation;  crab 


neuronal  identity  depends  on  the  molecular  composition  that 
determines  morphology  and  function,  and  it  is  regulated  at 
multiple  levels  (Hobert  et  al.  2010).  The  distinct  activity 
pattern  of  individual  neurons  is  due  to  the  regulation  of 
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neuronal  excitability  as  determined  principally  by  the  activity 
of  voltage-gated  ionic  currents.  The  generation  of  these  cur¬ 
rents  is  the  result  of  the  pattern  of  gene  expression  of  voltage¬ 
gated  channels  as  well  as  the  posttranscriptional  and  posttrans- 
lational  events  that  influence  the  biophysical  properties  of  the 
channels  inserted  in  the  membrane  (Hille  2001).  Furthermore, 
extensive  plasticity  in  neuronal  excitability  results  from  both 
activity-dependent  (Cudmore  and  Turrigiano  2004;  Desai  et  al. 
1999;  Li  et  al.  2004;  Loebrich  and  Nedivi  2009;  Zhang  and 
Linden  2003)  and  activity-independent  mechanisms  (Khorkova  and 
Golowasch  2007;  MacLean  et  al.  2003,  2005).  These  mecha¬ 
nisms  include  forms  of  plasticity  of  intrinsic  excitability  (in¬ 
trinsic  plasticity)  that  either  alter  or  homeostatic  ally  maintain 
the  output  of  these  cells  depending  on  conditions  (Zhang  and 
Linden  2003).  Recently,  the  coordinated  expression  of  conduc¬ 
tances  (or  conductance  properties)  and  the  contribution  of  such 
coordination  to  neuronal  functions  have  been  reported  in  an 
increasing  number  of  systems  (Bergquist  et  al.  2010;  Coggan  et 
al.  2010;  Del  Negro  et  al.  2002;  Khorkova  and  Golowasch 
2007;  MacLean  et  al.  2003;  McAnelly  and  Zakon  2000;  Schulz 
et  al.  2006,  2007).  At  this  point,  however,  neither  the  role  of 
this  coregulation  of  ionic  currents  nor  the  level  at  which 
correlations  are  controlled  is  well  known.  It  also  is  not  known 
if,  within  a  functional  network,  ionic  currents  are  globally 
regulated  among  all  functionally  related  neurons  or  if  they  are 
regulated  at  the  individual  neuron  level.  In  this  study,  we 
examined  at  what  level  these  current  correlations  are  con¬ 
trolled.  We  did  this  by  looking  at  the  correlated  distributions  of 
ionic  currents  and  of  the  mRNA  that  code  for  the  channels  that 
carry  these  currents  in  two  identified  neurons.  We  used  the 
well-defined  pyloric  neuronal  network  of  the  crab  Cancer 
borealis,  in  which  individual  neurons  can  unambiguously  be 
identified. 

The  stomatogastric  ganglion  (STG)  of  C.  borealis  contains 
25-26  neurons  (Kilman  and  Marder  1996)  that  are  identifiable 
on  the  basis  of  their  anatomic  projections  (Maynard  and  Dando 
1974).  These  neurons  belong  to  two  distinct  neuronal  networks 
that  form  central  pattern  generators  (CPGs)  producing  different 
rhythmic  motor  patterns.  The  faster  pyloric  rhythm  (cycle 
frequency  of  0. 5-2.0  Hz)  is  composed  of  six  neuronal  subtypes 
that  produce  rhythmic  activity  driven  by  the  pacemaker  activity 
of  one  neuron  (the  anterior  burster,  AB  neuron),  which  is 
strongly  electrically  coupled  to  two  pyloric  dilator  (PD)  neu¬ 
rons.  Thus  this  group  of  three  neurons  acts  as  a  single  pace¬ 
maker  unit.  This  rhythmic  activity  pattern  is  conditionally 
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dependent  on  the  release  of  neuromodulators  from  neurons 
located  in  adjacent  ganglia  that  project  into  the  STG  (Luther  et 
al.  2003;  Thoby-Brisson  and  Simmers  2002).  These  substances 
have  paracrine  actions  on  pyloric  neurons  and  are  thought  to 
regulate  activity  via  short-  and  long-term  actions  (Khorkova 
and  Golowasch  2007;  Swensen  and  Marder  2001). 

Individual  neurons  of  the  STG  CPGs  have  distinct  biophys¬ 
ical  properties  as  well  as  gene  expression  profiles.  For  exam¬ 
ple,  individually  identified  neurons  of  the  pyloric  network  have 
highly  variable  kinetics  and  conductance  density  levels  of  the 
hyperpolarization-activated  current,  /H  (Peck  et  al.  2006),  and 
the  transient  K+  current,  IA  (Goldman  et  al.  2001;  Golowasch 
et  al.  1999;  Khorkova  and  Golowasch  2007;  Peck  et  al.  2001). 
Particular  pyloric  neurons  also  differ  in  the  quantitative  rela¬ 
tionships  of  expression  of  voltage-gated  channel  genes  (Baro  et 
al.  1997;  Schulz  et  al.  2006,  2007),  with  different  cell  types 
expressing  different  relative  amounts  of  voltage-gated  channel 
mRNA  and  distinct  correlations  of  channel  mRNA  levels 
(Schulz  et  al.  2007).  These  results  suggest  a  differential  regu¬ 
lation  (and  coregulation)  of  ion  channel  gene  transcription  in 
distinct  STG  neuron  types.  Together,  the  overall  regulation  of 
ion  channel  gene  expression  and  posttranscriptional  processing 
must  play  a  substantial  role  in  generating  the  unique  output  of 
each  cell  type. 

In  this  study,  we  characterized  the  correlations  of  ionic 
currents  expressed  by  two  key  neuronal  subtypes  of  the  pyloric 
network,  the  PD  and  lateral  pyloric  (LP)  neurons,  under  dif¬ 
ferent  neuromodulatory  conditions  using  a  combination  of 
electrophysiological  and  molecular  analyses.  We  also  exam¬ 
ined  the  role  of  synaptic  interactions  and  altered  activity  in  the 
modification  of  these  correlations. 

METHODS 

Adult  male  Jonah  crabs  (C.  borealis )  were  purchased  from  local 
fishermen  (Newark,  NJ)  and  kept  in  cooled  seawater  aquaria  at 
~12°C.  Animals  were  anesthetized  by  cooling  for  15-30  min  on  ice. 
The  foregut  was  removed,  and  the  STG,  with  a  portion  of  the  nerves 
attached,  were  isolated  as  previously  described  (Harris- Warrick  1992; 
Selverston  et  al.  1976). 

The  removal  of  neuromodulatory  input  (i.e.,  decentralization)  was 
accomplished  by  either  transection  of  the  stomatogastric  nerve  (stn)  or 


by  blocking  action  potential  transmission  along  the  stn,  which  carries 
all  known  neuromodulatory  inputs  from  adjacent  ganglia  (esophageal 
and  commissural  ganglia)  to  the  STG  (Fig.  IA).  To  block  transmis¬ 
sion,  a  Vaseline  well  was  built  around  the  desheathed  nerve,  and 
either  the  well  was  filled  with  isotonic  (750  mM)  sucrose  plus  10~6  M 
tetrodotoxin  (Biotium),  or  the  stn  was  cut  as  close  to  the  STG  as 
possible.  Preparations  in  which  experiments  lasted  more  than  12  h 
(such  as  decentralization  experiments)  were  maintained  in  organo¬ 
typic  culture  at  6°C  between  recording  sessions  during  which  the 
preparation  was  bathed  in  physiological  saline. 

Electrophysiology.  Neurons  were  identified  and  intracellularly  re¬ 
corded  from  exactly  as  described  previously  (Zhao  et  al.  2010). 
Extracellular  recordings  were  performed  with  differential  recordings 
using  an  AC  amplifier  (model  1700;  A-M  Systems)  and  two  tungsten 
wires,  one  placed  inside  a  Vaseline  well  built  around  a  motor  nerve 
(either  the  lateral  ventricular  nerve,  Ivn,  or  the  pyloric  dilator  nerve, 
pdn,  or  both)  and  the  other  placed  in  the  bath  outside  the  well.  Ionic 
currents  were  measured  at  11-13°C  in  two-electrode  voltage  clamp, 
with  current  injection  electrode  resistances  of  18-22  Mfl  and  voltage 
recording  electrode  resistances  of  20-28  MO.  K+  currents  were 
measured  in  standard  C.  borealis  physiological  saline.  The  high- 
threshold,  voltage-gated  outward  potassium  currents,  /HTK,  were  ac¬ 
tivated  from  a  holding  voltage  of  —40  mV.  The  voltage-gated  tran¬ 
sient  current,  IA.  was  activated  with  depolarizing  steps  from  a  holding 
voltage  of  —80  mV.  The  high-threshold  currents  activated  during  the 
/A  activation  protocol  were  removed  by  subtracting  the  currents 
previously  measured  from  a  holding  potential  of  —40  mV.  The 
hyperpolarization-activated  current,  was  measured  using  8-s-long 
hyperpolarizing  pulses  from  a  holding  potential  of  —40  mV. 

Ionic  conductance  of  /HTK  and  IA  was  determined  by  dividing  the 
current  levels  measured  at  +20  mV  by  the  driving  force,  assuming  a 
K+  reversal  potential  (EK)  of  —80  mV,  and  ionic  conductance  of  In 
was  similarly  determined  using  a  current  measured  at  —110  mV  and 
a  reversal  potential  (EH)  of  — 10  mV  (Haedo  and  Golowasch  2006; 
Khorkova  and  Golowasch  2007). 

After  ionic  currents  were  measured  under  control  conditions,  the 
preparations  were  either  treated  with  10~5  M  picrotoxin  or  decentral¬ 
ized  and  placed  at  6°C  for  24  h.  The  currents  were  measured  again 
using  the  same  protocols,  and  the  neurons  were  subsequently  pro¬ 
cessed  for  RNA  extraction  (see  below). 

mRNA  quantification.  Quantitative  PCR  was  performed  as  de¬ 
scribed  by  Schulz  et  al.  (2006).  Briefly,  total  RNA  was  isolated  using 
the  RNeasy  microcolumn-based  RNA  extraction  kit  (Qiagen,  Valen¬ 
cia,  CA),  reverse  transcribed  using  Superscript  III  reverse  transcrip¬ 
tase  (Invitrogen),  and  used  as  a  template  in  real-time  RT-PCR  with 


2s 


Pdn 


Fig.  1.  Stomatogastric  nervous  system  and  ac¬ 
tivity  changes  under  different  modulatory  con¬ 
ditions.  A\  the  stomatogastric  nervous  system. 
OG,  esophageal  ganglion;  CoG,  commissural 
ganglion;  STG,  stomatogastric  ganglion;  stn, 
stomatogastric  nerve;  Ivn,  lateral  ventricular 
nerve;  pdn,  pyloric  dilator  nerve;  mvn,  medial 
ventricular  nerve.  The  circle  around  stn  repre¬ 
sents  a  Vaseline  well  used  to  block  action  po¬ 
tential  transmission  (methods).  B:  effect  of  the 
removal  of  neuromodulatory  input  to  the  STG. 
Pyloric  activity  is  temporarily  lost  after  decen¬ 
tralization  (4  h).  Activity  is  slower  but  can  be 
identified  as  pyloric  24  h  after  decentralization. 
C:  effect  of  the  removal  of  glutamatergic  inhib¬ 
itory  synapses  with  picrotoxin  (PTX).  The  lat¬ 
eral  pyloric  (LP)  neuron  becomes,  and  remains, 
tonic  over  24  h. 
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SYBR  green  (SABiosciences,  Frederick,  MA).  Primers  specific  for 
real-time  PCR  detection  of  Shal,  BK-KCa,  H ,  and  18S  rRNA  with 
SYBR  green  were  developed  and  designed  using  Primer3  software 
and  are  the  same  as  previously  reported  (Schulz  et  al.  2006,  2007). 

Statistics.  Comparisons  of  overall  levels  of  mRNA  and  conduc¬ 
tance  density  were  made  using  a  one-way  ANOVA  followed  by  post 
hoc  Tukey’s  f-test  pairwise  comparisons  (Systat;  Aspire  Software 
International,  Ashburn,  VA).  Correlation  (parametric  Pearson’s  prod¬ 
uct  moment  and  nonparametric  Spearman’s  rank  correlation  test)  and 
linear  regression  analyses  were  carried  out  using  Sigmastat  (Systat). 
We  compared  correlation  coefficients  between  experimental  groups 
for  a  given  pairwise  correlation  by  performing  a  Fisher  transformation 
of  p  using 


P' 


(0.5)loge 


1  +  P 

1  -  p 


and  then  calculating  a  Fisher’s  z  statistic  as  follows: 


z  = 


Pi  Pi 


1 


n2  —  3 


Finally,  slopes  of  relationships  between  conductances  and  between 
mRNA  levels  were  compared  between  experimental  groups  with  a 
modified  ANCOVA  (GraphPad  Prism  version  5;  GraphPad  Software, 
San  Diego,  CA). 


RESULTS 

As  has  been  shown  previously  (Luther  et  al.  2003;  Thoby- 
Brisson  and  Simmers  2002),  decentralization  of  the  STG  leads 
to  the  rapid  but  reversible  loss  of  rhythmic  activity  of  the 
pyloric  network  (Fig.  IB).  The  rhythmic  pyloric  activity  re¬ 
covers  to  levels  and  characteristics  similar  to  those  observed 
under  control  conditions  with  a  delay  of  hours  to  days  (Luther 
et  al.  2003).  The  removal  of  most  chemical  synapses  in  the 
ganglion,  which  are  primarily  inhibitory  glutamatergic,  was 
accomplished  using  10~5  M  picrotoxin,  which  disrupts  the 
triphasic  pyloric  rhythm  (Bidaut  1980;  Marder  and  Eisen 
1984),  making  the  follower  neurons  (LP  and  other  neurons)  fire 
tonically,  but  leaves  the  pacemaker’s  activity  (monitored  by 
the  PD  neurons’  activity  on  both  the  Ivn  and  pdn )  intact  over 
the  course  of  24  h  (Fig.  1C).  Picrotoxin  blocks  inhibitory 
chemical  synapses,  which  produces  drastic  changes  in  activity 
of  LP  but  not  PD  neurons.  This  allows  us  to  decouple  activity- 
dependent  mechanisms  from  neuromodulator-dependent  corre¬ 
lations  of  mRNA  and  conductance  levels. 

Ion  channel  conductances  are  correlated  across  channel 
types.  We  hypothesized  that,  as  shown  previously  for  PD 
neurons  (Khorkova  and  Golowasch  2007),  ionic  conductances 
are  expressed  in  other  pyloric  neurons  in  a  correlated  fashion. 
To  test  this,  we  measured  three  ionic  currents,  /HTK,  IA,  and  Ih 
in  PD  and  LP  neurons  and  estimated  the  ionic  conductances 
(<?htk’  8 A’  and  8h)-  Control  PD  neurons  expressed  pairwise 
correlations  (Fig.  2,  top  row)  very  similar  to  those  described 
previously  (Khorkova  and  Golowasch  2007).  Table  1  lists  all 
coefficients  of  determination  (R2)  and  Pearson  product  moment 
coefficients  (p)  for  each  pairwise  conductance  comparison.  We 
performed  the  same  measurements  in  another  identified  neu¬ 
ron,  the  LP  neuron,  and  found  that  these  three  conductances 
were  all  similarly  correlated  in  a  pairwise  manner  (Fig.  3,  top 
row,  and  Table  1).  The  nonparametric  Spearman’s  rank  corre¬ 
lation  test  was  also  used  and  gave  nearly  identical  results  to  the 


Pearson  product  moment  correlation  test;  thus  we  report  only 
the  latter  tests  in  this  report. 

To  determine  what  factors  may  control  these  correlations, 
we  treated  these  preparations  in  two  different  ways  designed  to 
disrupt  the  major  inputs  these  cells  receive  within  the  STG:  in 
one  we  used  picrotoxin  to  block  the  synaptic  input  these  two 
cells  receive,  and  in  the  other  we  blocked  all  neuromodulatory 
input  to  the  cells  (decentralization).  For  the  picrotoxin  group, 
we  treated  the  entire  stomatogastric  nervous  system  (Fig.  1A) 
with  10~5  M  picrotoxin. 

In  this  process  of  removing  synaptic  inhibition  to  LP,  we 
greatly  alter  its  firing  pattern  (Fig.  1C).  Despite  these  manip¬ 
ulations,  we  observed  that  all  pairwise  correlations  of  the  three 
ionic  currents  were  maintained  with  similar  correlation  coeffi¬ 
cients  in  both  PD  (Fig.  2)  and  LP  neurons  (Fig.  3,  Table  1). 
Even  though  the  gHTK't?H  Pair  f°r  the  LP  neurons  shows  a 
marginally  higher  P  value  (0.056)  than  the  cutoff  of  0.05,  we 
accepted  it  as  indicative  of  a  significant  correlation. 

For  the  “decentralized”  group,  significant  differences  be¬ 
tween  the  two  cell  types  were  observed  after  24  h:  in  PD 
neurons,  only  the  pairwise  gA  vs.  gH  correlation  was  main¬ 
tained,  whereas  the  gHTK  vs.  gA  and  the  gHTK  vs.  gH  correla¬ 
tions  were  lost  (Fig.  2,  Table  1);  in  contrast,  in  LP  neurons,  all 
of  the  correlations  were  maintained,  but  the  correlation  coef¬ 
ficient  for  gA  vs.  g 1 1  was  significantly  lower  in  decentralized 
relative  to  control  groups  (Fig.  3,  Table  1). 

In  addition  to  the  presence  or  absence  of  a  significant 
correlation  between  conductances,  we  determined  whether  cor¬ 
relations  that  were  maintained  between  conditions  (e.g.,  the 
same  correlation  present  in  control  and  picrotoxin  groups) 
differed  in  the  slope  of  the  relationship.  Despite  the  fact  that 
correlations  were  maintained  among  conductance  pairs  across 
experimental  groups,  there  were  significant  changes  in  a  subset 
of  these  relationships  in  both  PD  and  LP  cells  (Table  2).  There 
was  a  significant  decrease  in  the  slope  of  gHTK  vs.  gA  between 
the  control  and  picrotoxin-treated  groups,  as  well  as  a  signif¬ 
icant  increase  in  the  gA  vs.  gH  slope  of  the  decentralized 
relative  to  the  control  group  in  PD  cells.  However,  there  was  no 
change  in  the  slope  of  the  relationship  between  gHTK  and  gH  in 
control  and  picrotoxin-treated  PD  cells.  For  LP  cells,  the  slope 
of  the  gHTK  vs.  gA  relationship  was  lower  in  picrotoxin-treated 
LP  cells,  whereas  the  slope  of  all  three  relationships  was  lower 
in  decentralized  LP  cells  relative  to  controls  (Table  2). 

Ion  channel  correlations  at  the  mRNA  level:  distinct  from 
conductance  correlations.  Pairwise  comparison  of  channel 
mRNA  levels  under  control  conditions  revealed  a  three-way 
correlation  between  BK-KCa,  H,  and  Shal  in  PD  neurons  (Fig.  4, 
Table  3),  consistent  with  the  correlations  observed  at  the  ionic 
conductance  level  described  above.  Surprisingly,  LP  neurons 
showed  no  correlation  of  the  BK-KCa-Shal  pair,  whereas  the 
other  two  channel  mRNA  pairs  showed  strong  and  statistically 
significant  correlation  levels  (Fig.  5,  Table  3).  Again,  as  before, 
to  determine  what  factors  may  control  these  correlations,  we 
treated  one  group  of  preparations  with  1 0  5  M  picrotoxin  for 
24  h,  decentralized  a  second  group  for  24  h,  and  measured 
channel  mRNA  levels  in  both.  We  observed  that  PD  neurons 
treated  with  picrotoxin  showed  similarly  strong  and  statisti¬ 
cally  significant  correlations  in  all  three  pairs  as  in  the 
control  preparations  (Fig.  4,  Table  3).  Consistent  with  the 
correlations  observed  under  control  conditions,  in  LP  neu¬ 
rons  treated  with  picrotoxin,  the  BK-KCa-Shal  pair  showed 
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Fig.  2.  Ionic  conductance  correlations  in  pyloric  dilator  (PD)  neurons.  Conductances  of  the  high-threshold  K+  current  (yHTK).  transient  K+  current  (gA),  and 
hyperpolarization-activated  inward  current  (gH)  are  graphed  in  all  pairwise  combinations.  Each  point  represents  a  different  cell.  Regression  lines  are  shown  only 
for  significant  correlations  (Pearson  correlation  analysis,  P  <  0.05).  Measurements  were  obtained  under  the  3  conditions  indicated:  control  ( top  row),  after  24 
h  in  HU5  M  PTX  (middle  row),  and  24  h  after  decentralization  ( bottom  row).  Ionic  conductances  are  reported  in  p.S. 


no  correlation,  whereas  the  other  two  channel  mRNA  pairs 
did  (Fig.  5,  Table  3). 

Finally,  we  evaluated  the  pairwise  channel  mRNA  correla¬ 
tion  levels  in  the  decentralized  preparation  group  and  found 
another  surprise:  in  the  LP  neurons,  the  BK-KCa-Shal  pair  that 


had  shown  no  statistically  significant  correlations  in  either  the 
control  or  picrotoxin  group  now  showed  a  moderate  but  sta¬ 
tistically  significant  correlation.  The  other  two  channel  mRNA 
pairs  ( BK-KCa-H  and  Shal-H)  on  the  other  hand,  were  mixed. 
There  was  no  significant  change  in  BK-KCa-H  correlation,  but 


Table  1.  Ion  channel  conductance-level  correlations 


Shtk  vs.  gA 

gHTK  vs.  gH 

Sa  vs.  g„ 

R2 

p  {P  value) 

n 

R2 

p  ( P  value) 

n 

R2 

p  (P  value) 

n 

Controls 

0.481 

0.693"  (<  0.001) 

95 

0.376 

PD  neuron  conductances 

0.613a(  <0.001) 

94 

0.409 

0.640“  (<  0.001) 

89 

Picrotoxin 

0.397 

0.630"  (0.028) 

12 

0.517 

0.719"  (0.019) 

10 

0.422 

0.722"  (0.050) 

9 

Decentralized 

0.055 

0.234  (0.102) 

50 

0.009 

0.094  (0.562) 

40 

0.517 

0.719"  K0.001) 

41 

Controls 

0.631 

0.794"  (<  0.001) 

76 

0.405 

LP  neuron  conductances 

0.637a(<  0.001) 

69 

0.514 

0.717"b  (<0.001) 

70 

Picrotoxin 

0.396 

0.630"  (0.038) 

11 

0.383 

0.619“  (0.056) 

10 

0.795 

0.892“  (<  0.001) 

10 

Decentralized 

0.443 

0.665"  (<  0.001) 

31 

0.213 

0.462“  (0.022) 

24 

0.276 

0.526b  (0.006) 

26 

Coefficients  of  determination  ( R 2)  and  Pearson  product  moment  correlation  coefficients  (p)  of  pyloric  dilator  (PD)  and  lateral  pyloric  (LP)  neurons  (n  =  no. 
of  neurons)  were  determined  under  control,  picrotoxin,  and  decentralized  conditions  for  pairs  of  ion  channel  conductance  levels:  gHTK,  high-threshold  K+  current 
conductance;  gA,  transient  K+  current  conductance;  gH,  hyperpolarization-activated  inward  current  conductance.  Boldface  indicates  a  significant  (P  <  0.05) 
correlation  in  the  Pearson’s  correlation  analysis.  a,bDifferent  superscript  letters  between  2  groups  indicate  significant  differences  among  correlation  coefficients 
of  these  groups.  Correlation  coefficient  analysis  was  not  performed  with  groups  in  which  a  significant  correlation  was  not  found. 
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Fig.  3.  Ionic  conductance  correlations  in  LP  neurons.  Conductances  gHTK,  8a>  and  £h  are  graphed  in  all  pairwise  combinations.  Each  point  represents  a  different 
cell.  Regression  lines  are  shown  only  for  significant  correlations  (Pearson  correlation  analysis,  P  <  0.05).  Measurements  were  obtained  under  the  3  conditions 
indicated:  control  ( top  row),  after  24  h  in  10-5  M  PTX  {middle  row),  and  24  h  after  decentralization  {bottom  row).  Ionic  conductances  are  reported  in  fx S. 


there  was  a  significant  decrease  in  the  correlation  coefficient  of 
the  Shal-H  relationship  (Fig.  5,  Table  3). 

We  also  determined  whether  the  slopes  of  the  relationships 
changed  for  mRNA  correlations  across  treatment  groups 
(Table  4).  There  was  a  significant  decrease  in  both  the  BK-KCa 
vs.  H  and  Shal  vs.  H  relationships  from  control  to  the  picrotoxin 

Table  2.  Slopes  for  pairwise  conductance-level  relationships 


group,  whereas  the  slope  of  the  relationship  between  BK-KCa  and 
Shal  was  maintained  in  PD  cells.  In  decentralized  PD  cells,  there 
was  a  significant  decrease  in  the  slope  of  the  BK-KCa  vs.  Shal  and 
BK-KCa  vs.  H  relationships  relative  to  control,  and  the  Shal  vs.  H 
relationship  was  lost.  For  LP  cells,  there  was  a  significant  increase 
in  the  Shal  vs.  H  relationship  in  picrotoxin,  whereas  the  BK-KCa 


Shtk  vs.  gA 

gHTK  VS.  gH 

g A  VS.  gH 

Slope 

SD 

n 

Slope 

SD 

n 

Slope 

SD 

n 

PD  neuron 

conductances 

Controls 

0.8875“ 

0.10 

95 

0.087“ 

0.011 

94 

0.013“ 

0.012 

91 

Picrotoxin 

0.2765b 

0.11 

12 

0.048“ 

0.017 

10 

0.064“'b 

0.028 

9 

Decentralized 

NS 

NS 

50 

NS 

NS 

40 

0.088b 

0.014 

41 

LP  neuron 

conductances 

Controls 

0.68" 

0.06 

76 

0.076“ 

0.011 

69 

0.104“ 

0.012 

70 

Picrotoxin 

0.23b 

0.09 

11 

0.042",b 

0.019 

10 

0.16“ 

0.029 

10 

Decentralized 

0.35b 

0.07 

31 

0.027b 

0.011 

24 

0.055b 

0.018 

26 

Conductances  were  measured  under  control,  picrotoxin,  and  decentralized  conditions,  and  slopes  for  significantly  correlated  conductances  were  analyzed  across 
groups  via  modified  analysis  of  covariance  (ANCOVA).  Boldface  indicates  comparisons  in  which  the  conditions  result  in  significantly  different  slopes  from  one 
another  {P  <  0.05).  ^Different  superscript  letters  for  slope  values  indicate  pairs  for  which  the  difference  is  significant  (thus  pairs  with  the  same  superscript  are 
not  different).  NS,  no  significant  pairwise  correlation  was  found  between  the  2  conductance  levels  in  a  given  treatment  group. 
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Fig.  4.  Channel  mRNA  correlations  in  PD  neurons.  mRNA  levels  of  BK-KCa,  Shell,  and  II  are  graphed  in  all  pairwise  combinations.  Each  point  represents  a 
different  cell.  Regression  lines  are  shown  only  for  significant  correlations  (Pearson  correlation  analysis,  P  <  0.05).  Measurements  were  obtained  under  the  3 
conditions  indicated:  control  (top  row),  after  24  h  in  1 0  5  M  PTX  ( middle  row),  and  24  h  after  decentralization  ( bottom  row).  Numbers  indicate  mRNA  copy 
numbers. 


vs.  H  relationship  was  unchanged.  No  significant  changes  were 
seen  in  the  slopes  of  mRNA  relationships  relative  to  control  in 
decentralized  LP  cells,  although  a  novel  BK-KCa  vs.  Shal  rela¬ 
tionship  was  detected. 

DISCUSSION 

Neurons  have  the  ability  to  generate  a  functional  output  that 
is  stable  over  long  time  scales.  To  ensure  this,  neurons  express 


plasticity  in  their  intrinsic  properties  and  excitability  (Daoudal 
and  Debanne  2003;  Davis  2006;  Frick  and  Johnston  2005; 
Haedo  and  Golowasch  2006;  Thoby-Brisson  and  Simmers 
2002;  Turrigiano  and  Nelson  2004;  Zhang  and  Linden  2003), 
which  manifests  itself  via  changes  in  ionic  channel  properties 
(Zhang  and  Linden  2003),  but  the  levels  at  which  this  is 
regulated,  and  the  molecular  mechanisms  involved,  remain  to 
be  understood.  It  has  been  shown  previously  that  ionic  chan- 


Table  3.  Ion  channel  mRNA-level  correlations 


BK-KCa  vs.  Shal 

BK-KCa  vs.  H 

Shal  vs.  H 

R2 

p  ( P  value) 

n 

R2 

p  (P  value) 

n 

R2 

p  ( P  value) 

n 

Controls 

0.848 

0.921°  (0.001) 

9 

0.755 

PD  neuron  mRNA 

0.869a  (0.005) 

8 

0.847 

0.920°  (<  0.001) 

9 

Picrotoxin 

0.886 

0.941°  (C0.001) 

12 

0.791 

0.889°  (<  0.001) 

11 

0.728 

0.853°  (<  0.001) 

12 

Decentralized 

0.238 

0.4871’  (0.021) 

22 

0.303 

0.550b  (0.006) 

23 

0.014 

0.118  (0.600) 

22 

Controls 

0.089 

0.298  (0.216) 

19 

0.602 

LP  neuron  mRNA 

0.7051  (0.001) 

16 

0.782 

0.884°  (<  0.001) 

18 

Picrotoxin 

0.034 

0.186  (0.585) 

11 

0.649 

0.806°  (0.005) 

10 

0.932 

0.965"  (<  0.001) 

9 

Decentralized 

0.357 

0.598  (0.015) 

16 

0.328 

0.573“  (0.020) 

16 

0.328 

0.573b  (0.020) 

16 

R 2  and  p  of  PD  and  LP  neurons  were  determined  under  control,  picrotoxin,  and  decentralized  conditions  for  pairs  of  ion  channel  mRNA  levels.  Boldface 
indicates  a  significant  (P  <  0.05)  correlation  in  the  Pearson’s  correlation  analysis.  abDifferent  superscript  letters  between  2  groups  indicate  significant 
differences  among  correlation  coefficients  of  these  groups.  Correlation  coefficient  analysis  was  not  performed  with  groups  in  which  a  significant  correlation  was 
not  found. 
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Fig.  5.  Channel  mRNA  correlations  in  LP  neurons.  mRNA  levels  of  BK-KCa ,  Shal,  and  H  are  graphed  in  all  pairwise  combinations.  Each  point  represents  a 
different  cell.  Regression  lines  are  shown  only  for  significant  correlations  (Pearson  correlation  analysis,  P  <  0.05).  Measurements  were  obtained  under  the  three 
conditions  indicated:  control  ( top  row),  after  24  h  in  1 0  5  M  PTX  ( middle  row),  and  24  h  after  decentralization  ( bottom  row).  Numbers  indicate  mRNA  copy 
numbers. 


nels  may  be  codependent  at  both  the  mRNA  (Schulz  et  al. 
2007)  and  the  conductance  levels  (Khorkova  and  Golowasch 
2007;  MacLean  et  al.  2003)  and  that  such  codependence  may 
be  related  to  the  homeostatic  regulation  of  function  in  the 
pyloric  network  (Ball  et  al.  2010;  Burdakov  2005;  Franklin  et 
al.  2010;  Khorkova  and  Golowasch  2007;  MacLean  et  al. 
2005).  The  purpose  of  this  study  was  twofold:  first,  to  deter¬ 


mine  whether  the  apparent  coregulation  of  ionic  channel  ex¬ 
pression  is  cell-type  specific,  and  second,  to  determine  how 
neuromodulatory  input  and  synaptic  activity  affects  the  corre¬ 
lated  expression  of  ionic  channels  at  both  the  ionic  conduc¬ 
tance  and  mRNA  levels.  We  have  shown  that  1 )  neurons 
comodify  ion  channel  mRNA  ( BK-KCa ,  Slial,  and  If)  levels  in 
a  cell-specific  manner:  under  control  conditions  (i.e.,  neuro- 


Table  4.  Slopes  for  pairwise  ion  channel  mRNA-level  relationships 


BK-KCa  vs.  Shal 

BK-KCa  vs.  H 

Shal  vs.  H 

Slope  SD 

n 

Slope 

SD 

n 

Slope 

SD 

n 

Controls 

0.39“ 

0.07 

Picrotoxin 

0.47“ 

0.05 

Decentralized 

0.14b 

0.06 

Controls 

NS 

NS 

Picrotoxin 

NS 

NS 

Decentralized 

0.43 

0.05 

PD  neuron  mRNA 


8 

0.43" 

0.10 

12 

0.15b 

0.03 

22 

0.14b 

0.05 

LP  neuron  mRNA 

19 

0.35" 

0.88 

11 

0.42" 

0.11 

16 

0.21" 

0.08 

8 

1.10" 

0.18 

9 

11 

0.28b 

0.05 

12 

23 

NS 

NS 

22 

18 

0.57" 

0.08 

18 

10 

1.74b 

0.18 

9 

16 

1.1 4"b 

0.44 

16 

mRNA  levels  were  measured  under  control,  picrotoxin,  and  decentralized  conditions,  and  slopes  for  significantly  correlated  conductances  were  analyzed  across 
groups  via  modified  ANCOVA.  Boldface  indicates  comparisons  in  which  the  conditions  result  in  significantly  different  slopes  from  one  another  ( P  <  0.05). 
a'bDifferent  superscript  letters  for  slope  values  indicate  pairs  for  which  the  difference  is  significant  (thus  pairs  with  the  same  superscript  are  not  different).  NS, 
no  significant  pairwise  correlation  was  found  between  the  2  mRNA  levels  in  a  given  treatment  group. 
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modulatory  and  synaptic  input  present),  PD  neurons  express 
these  three  mRNA  species  in  a  three-way  correlated  manner, 
whereas  LP  neurons  appear  to  do  so  only  for  the  pairs  Shal  vs. 
H  and  BK-KCa  vs.  H  and  not  for  the  BK-KCa  vs.  Shal  pair.  2) 
Removal  of  neuromodulatory  input  can  destroy,  induce,  or 
alter  the  strength  or  slope  of  the  correlated  expression  of  ion 
channel  mRNA  in  a  cell-specific  manner:  in  PD  neurons,  all 
correlated  pairs  either  lose  or  have  significantly  weakened 
correlations,  whereas  LP  neurons  show  a  gain  in  the  correlation 
of  BK-KCa  vs.  Shal  but  weakening  of  the  correlations  of  the 
other  two  pairs  examined.  3)  At  the  conductance  level,  three- 
way  pairwise  correlations  appear  to  be  common  to  both  cell 
types  examined,  but  the  response  to  the  loss  of  neuromodula¬ 
tory  input  is  cell-type  specific  (i.e.,  different  between  PD  and 
LP  neurons)  and  is  also  distinct  from  that  observed  at  the 
mRNA  level.  4)  Finally,  chemical  synaptic  activity  does  not 
seem  to  play  a  role,  at  least  for  the  two  cell  types  examined,  in 
controlling  the  correlated  expression  of  the  ionic  currents  we 
measured  at  either  the  mRNA  or  the  conductance  levels.  Of 
these  results,  only  the  three-way  pairwise  correlations  and  the 
absence  of  activity-dependent  regulation  of  ionic  conductances 
in  PD  neurons  have  been  shown  previously  (Khorkova  and 
Golowasch  2007).  Yet,  it  was  previously  unclear  whether  these 
effects  acted  solely  at  the  level  of  membrane  conductances  or 
whether  they  extended  to  multiple  levels  of  cellular  organiza¬ 
tion  such  as  mRNA  for  ion  channels.  Furthermore,  it  was  not 
known  whether  these  relationships  would  exist  in  multiple 
neurons  in  the  network,  including  a  “follower”  neuron  such  as 
LP,  or  be  limited  to  critical  pacemaker  cells  such  as  PD. 

There  is  a  substantial  overlap  in  correlations  among  conduc¬ 
tances  and  those  among  mRNA  levels;  in  PD  cells,  correlations 
were  consistent  in  six  of  nine  relationships  between  conduc¬ 
tances  and  mRNA,  whereas  in  LP  cells,  seven  of  nine  were 
consistent.  These  results  therefore  implicate  a  functional  link 
between  regulation  at  the  mRNA  level  and  membrane  conduc¬ 
tance,  which  has  previously  been  seen  for  LP  but  not  PD 
neurons  (Schulz  et  al.  2006).  Yet,  one  of  the  more  striking 
outcomes  of  our  study  is  that  the  correlations  seen  among 
channel  mRNA  types  and  those  among  ionic  conductances, 
although  both  seemingly  influenced  strongly  by  neuromodula¬ 
tion,  can  themselves  be  distinct.  These  results  indicate  that 
neuromodulation  plays  a  role,  at  least  in  part,  in  coordinating 
regulation  across  multiple  levels  of  cellular  organization.  Nev¬ 
ertheless,  our  work  suggests  that  modulators  can  have  distinct 
and  independent  effects  on,  and  perhaps  independent  pathways 
for,  regulation  of  transcription  and  posttranslational  mecha¬ 
nisms  involved  in  generating  membrane  currents.  This  is  con¬ 
sistent  with  recent  studies  demonstrating  variation  in  similar 
kinds  of  relationships  between  mRNA  and  protein,  presumably 
as  a  result  of  transcriptional  and  nontranscription-level  varia¬ 
tion  (Foss  et  al.  2011).  Although  it  is  an  obvious  notion  that 
mRNA  abundance  does  not  simply  equate  to  mature  protein  or 
channel  activity,  this  is,  to  our  knowledge,  the  first  demonstra¬ 
tion  of  separate  control  by  neuromodulators  of  membrane  ion 
channel  conductances  and  their  mRNA.  Given  that  our  manip¬ 
ulations  remove  the  input  to  the  STG  of  a  large  number  of 
neuromodulators  [at  last  count,  at  least  27  (Marder  and  Bucher 
2007)],  it  is  conceivable  that  different  modulators  have  distinct 
actions  at  different  levels.  However,  it  is  interesting  to  note  that 
a  single  peptidergic  modulator  (proctolin)  has  been  shown  to 
revert  fully,  in  an  activity-independent  manner,  all  the  effects 


of  decentralization  on  conductance  correlations  in  PD  neurons 
(Khorkova  and  Golowasch  2007).  Our  results  are  significant  in 
that  they  reveal  important  differences  in  how  the  correlated 
expression  of  ionic  currents  ultimately  comes  about,  suggest¬ 
ing  multiple  levels  of  potential  regulation:  at  the  very  least,  at 
the  levels  of  transcription,  and  also  at  some  level  between  the 
beginning  of  translation  and  the  full  integration  of  the  channel 
proteins  into  the  plasma  membrane. 

What  are  the  possible  mechanisms  underlying  this  phenom¬ 
enon?  Cotranscription  is  a  process  by  which  the  close  proxim¬ 
ity  of  molecular  elements  and  genes  involved  in  transcription 
increases  the  likelihood  of  them  being  simultaneously  ex¬ 
pressed  (Sutherland  and  Bickmore  2009).  This  would  result  in 
ion  channels  appearing  to  be  “tracking”  each  other’s  expres¬ 
sion  levels.  If  transcription  is  the  point  at  which  this  process  of 
coexpression  occurs,  the  many  modifications  and  variations 
that  can  occur  along  the  way  as  the  ion  channel  transcripts  are 
processed,  transported,  and  translated,  and  the  expressed  chan¬ 
nel  proteins  are  then  modified  and  exported  into  the  plasma 
membrane,  could  explain  the  lack  of  perfect  inter-ion  channel 
correlation  across  preparations  (Tables  1  and  3).  Furthermore, 
this  could  provide  an  explanation  for  the  lower  level  of  corre¬ 
lation  observed  in  most  pairs  at  the  conductance  rather  than  at 
the  mRNA  level  (when  significant  correlations  exist).  Never¬ 
theless,  it  appears  that  correlations  in  expression  may  also  be 
implemented  downstream  of  transcription.  For  instance,  in  LP 
neurons  (control  conditions),  the  BK-KCa  vs.  Shal  pair,  and  in 
PD  neurons  (decentralized  conditions),  the  Shal  vs.  H  pair, 
show  no  apparent  correlation  at  the  transcription  level  but  a 
high  level  of  correlation  at  the  conductance  level.  Mechanisms 
to  explain  this  could  be  the  cotranslational  interaction  of  ion 
channels  similar  to  what  occurs  during  ion  channel  subunit 
assembly  (Shi  et  al.  1996),  the  cotrafficking  into  the  plasma 
membrane  (Vanoye  et  al.  2010),  the  coassembly  of  ion  chan¬ 
nels  into  macromolecular  complexes  (e.g.,  Arcangeli  2011; 
Frank  2011;  Zhang  et  al.  2011),  or  the  direct  interactions 
between  ion  channels  via  nonconducting  properties  (Kacz- 
marek  2006).  The  results  observed  by  MacLean  et  al.  (2003)  in 
which  the  transfection  of  lobster  PD  neurons  with  Shal  mRNA 
resulted  in  the  activity-independent  enhanced  coexpression  of 
the  H  current,  suggest  that  the  control  of  coexpression  does 
occur  at  a  translational  or  posttranslation  level.  We  did  not 
measure  transcription  per  se,  but  only  steady-state  mRNA 
levels,  and  thus  the  correlations  of  mRNA  types  could  be  the 
result  of  posttranscriptional  regulation  of  mRNA  stability. 
Therefore,  our  results  indicate  that  neuromodulators  may  reg¬ 
ulate  either  transcription  directly  or  the  stability  of  mRNA  in  a 
concerted  manner. 

The  loss  of  cotranscription  observed  in  PD  neurons  after 
decentralization  could  be  attributed  to  neuromodulators  acti¬ 
vating  a  transcription  regulation  pathway  that  affects  the  co¬ 
expression  of  ion  channels  at  any  of  the  steps  mentioned  above. 
Interestingly,  our  data  suggest  that  in  different  neuronal  cell 
types,  this  process  is  different  (i.e.,  inverse,  such  as  in  LP 
neurons),  suggesting  that  neuromodulators  may  also  act  to 
repress  cotranscription.  Although  activity  is  known  to  influ¬ 
ence  the  expression  of  ion  channels  in  a  wide  variety  of 
systems,  including  the  pyloric  network  (Cudmore  and  Turri- 
giano  2004;  Desai  et  al.  1999;  Golowasch  et  al.  1999;  Haedo 
and  Golowasch  2006;  Li  et  al.  2004;  Loebrich  and  Nedivi 
2009;  Turrigiano  et  al.  1994;  Zhang  and  Linden  2003),  the 
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coexpression  of  the  ion  channels  examined  in  this  work  ap¬ 
pears  not  to  be  influenced  by  synaptic  activity.  Our  results  do 
not  rigorously  rule  out  activity  as  a  factor,  but  STG  neurons  are 
known  to  be  sensitive  to  patterned  activity  (Golowasch  et  al. 
1999;  Haedo  and  Golowasch  2006;  Turrigiano  et  al.  1994),  and 
synaptic  activity  is  known  to  dramatically  influence  these 
patterns  of  activity  in  certain  neurons  (LP  neurons,  for  exam¬ 
ple,  see  Fig.  1).  Thus  our  results  are  consistent  with  a  lack  of 
effect  of  activity  on  the  correlated  expression  of  ion  channels, 
which  was  also  demonstrated  previously  for  PD  neurons 
(Khorkova  and  Golowasch  2007). 

What  functional  role  does  coexpression  of  ion  channels 
serve?  Given  the  high  level  of  variability  in  the  conductances 
of  distinct  ionic  currents  between  neurons  of  the  same  type  (see 
range  along  the  axes  of  any  of  our  graphs  herein),  it  appears 
that  coexpression  of  ionic  currents  is  one  way  to  reduce  the 
global  variability  of  the  system  by  linking  multiple  highly 
variable  currents  to  each  other.  It  has  been  shown  theoretically 
that  this  can  indeed  stabilize  specific  activity  features  (Ball  et 
al.  2010;  Burdakov  2005;  Franklin  et  al.  2010;  MacLean  et  al. 
2005).  The  loss  of  correlations  between  ionic  conductances 
may  be  an  evolutionary  adaptation  of  the  neurons  in  this 
system  to  permit  the  independent  modification  of  their  levels  as 
the  pyloric  network  recovers  its  lost  activity.  By  modifying  the 
process  that  keeps  these  conductances  linked  to  each  other, 
constraints  on  the  system  may  be  eased  to  make  it  better  able 
to  visit  alternative  activity  states  that  may  result  in  the  resto¬ 
ration  of  rhythmic  neuronal  activity  and  behavioral  output. 
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Neurons  and  networks  undergo  a  process  of  homeostatic  plasticity  that  stabilizes  output  by  integrating  activity  levels  with  network  and 
cellular  properties  to  counter  longer-term  perturbations.  Here  we  describe  a  rapid  compensatory  interaction  among  a  pair  of  potassium 
currents,  IA  and  IKCa,  that  stabilizes  both  intrinsic  excitability  and  network  function  in  the  cardiac  ganglion  of  the  crab,  Cancer  borealis. 
We  determined  that  mRNA  levels  in  single  identified  neurons  for  the  channels  which  encode  IA  and  IKCa  are  positively  correlated,  yet  the 
ionic  currents  themselves  are  negatively  correlated,  across  a  population  of  motor  neurons.  We  then  determined  that  these  currents  are 
functionally  coupled;  decreasing  levels  of  either  current  within  a  neuron  causes  a  rapid  increase  in  the  other.  This  functional  interdepen¬ 
dence  results  in  homeostatic  stabilization  of  both  the  individual  neuronal  and  the  network  output.  Furthermore,  these  compensatory 
increases  are  mechanistically  independent,  suggesting  robustness  in  the  maintenance  of  neural  network  output  that  is  critical  for 
survival.  Together,  we  generate  a  complete  model  for  homeostatic  plasticity  from  mRNA  to  network  output  where  rapid  post- 
translational  compensatory  mechanisms  acting  on  a  reservoir  of  channels  proteins  regulated  at  the  level  of  gene  expression  provide 
homeostatic  stabilization  of  both  cellular  and  network  activity. 


Introduction 

The  balance  of  plasticity  and  stability  in  generating  appropriate 
output  is  a  matter  of  fundamental  importance  in  the  nervous 
system  across  all  functional  levels.  These  processes  occur  even  at 
the  most  fundamental  level,  the  excitability  of  individual  neu¬ 
rons,  and  yet  little  is  known  about  mechanisms  governing  these 
processes  (Marder,  2011;  Turrigiano,  2011).  Early  work  identi¬ 
fied  such  processes  of  “homeostatic  plasticity”  of  intrinsic  excit¬ 
ability  (LeMasson  et  al.,  1993;  Turrigiano  et  al.,  1994,  1995; 
Golowasch  et  al.,  1999),  but  subsequent  focus  more  intensely 
shifted  to  determining  how  stabilization  of  synapses  is  accom¬ 
plished  through  synaptic  scaling  (Turrigiano,  2012).  Recently,  a 
resurgence  of  interest  in  plasticity  of  intrinsic  excitability  has  ac¬ 
companied  work  on  synaptic  scaling  (Debanne  and  Poo,  2010; 
Misonou,  2010;  Turrigiano,  2011). 

Fewer  studies  on  homeostatic  plasticity  have  considered  func¬ 
tional  compensation  in  the  context  of  endogenous  network  ac¬ 
tivity.  The  most  dramatic  example  may  be  complete  recovery  of 
motor  network  output  following  loss  of  central  inputs  as  a  result 
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of  changes  in  conductances  in  the  crustacean  stomatogastric  gan¬ 
glion  (STG)  (Thoby-Brisson  and  Simmers,  1998,  2002).  Knock¬ 
outs  of  K  +  channels  in  mice  have  been  shown  to  have  modest 
effects  on  phenotype  and  cellular  output  as  a  result  of  compen¬ 
sation  by  other  K+  channels  (Guo  et  al.,  2005;  Nerbonne  et  al., 
2008).  Additionally,  overexpression  of  A-type  K+  channels  in 
STG  neurons  results  in  little  change  in  neuronal  output  as  a  result 
of  compensatory  increases  in  H-current  (MacLean  et  al.,  2003, 
2005).  However,  these  examples  feature  mechanisms  that  act 
over  longer  time  scales  of  days  to  weeks.  While  initial  reports  of 
plasticity  in  intrinsic  excitability  were  found  over  shorter  time 
scales  (Desai  et  al.,  1999;  Golowasch  et  al.,  1999),  surprisingly 
little  is  known  of  the  role  these  mechanisms  may  play  in  short¬ 
term  ongoing  activity  of  biologically  intact  networks,  specifically 
where  an  expectation  for  rapid  conservation  of  output  could  be 
argued,  such  as  in  central  pattern  generators  (CPGs). 

Compensation  in  CPG  circuits  may  be  inferred  from  the  fact 
that  normal  populations  of  unmanipulated  motor  neurons  of 
two  invertebrate  CPGs,  the  cardiac  and  stomatogastric  ganglia, 
show  correlations  in  expression  levels  of  mRNAs  for  ion  channels 
(Schulz  et  al.,  2007;  Tobin  et  al.,  2009)  and  membrane  conduc¬ 
tances  (Khorkova  and  Golowasch,  2007;  Temporal  et  al.,  2012). 
One  relationship  detected  in  previous  work  (Tobin  et  al.,  2009)  is 
a  positive  correlation  between  BKKCa  and  shaker  mRNA  levels  in 
neurons  of  the  cardiac  ganglion.  These  channels  encode  calcium- 
activated  and  A-type  K+  currents,  respectively.  However,  it  is 
unclear  how  a  neuron  would  use  two  similar  hyperpolarizing 
conductances  additively  to  generate  or  maintain  its  output.  In 
this  study  we  focused  on  elucidating  the  functional  relationship 
between  these  two  conductances  in  motor  neurons  of  the  crab 
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cardiac  ganglion.  We  discovered  a  striking  discrepancy  in  the 
relationship  of  these  channels  across  functional  levels:  mRNAs 
for  these  channels  were  positively  correlated,  while  their  conduc¬ 
tances  were  negatively  correlated.  Therefore,  we  propose  a  com¬ 
prehensive  hypothesis  for  plasticity  of  excitability  from  mRNA  to 
network  output  whereby  rapid  compensation  provides  stabiliza¬ 
tion  of  cellular  and  network  activity. 

Materials  and  Methods 

Preparations.  Cancer  borealis  crabs  of  either  sex  were  purchased  and 
shipped  overnight  from  The  Fresh  Lobster  Company  (Gloucester,  MA). 
Crabs  were  kept  between  24  h  and  2  weeks  in  artificial  sea  water  at  12°C 
before  use.  Crabs  were  anesthetized  in  ice  for  15  min  before  the  dissec¬ 
tion.  The  dissection  took  place  in  chilled  physiological  saline  comprised 
of  440  mM  NaCl,  26  mM  MgCl2  13  mM  CaCl2,  1 1.2  mM  Trizma  base,  1 1 
mM  KC1,  and  5  mM  maleic  acid  (pH  =  7.4).  When  we  wanted  to  isolate 
individual  large  cells,  individual  strands  of  bulking  nylon  were  used  to 
ligate  the  nerve  on  both  sides  of  a  large  cell  soma.  To  impale  large  cells, 
each  cell  was  individually  desheathed  using  a  tungsten  needle  (Fine  Sci¬ 
ence  Tools). 

Quantitative  single-cell  RT-PCR.  Quantitative  RT-PCR  was  performed 
as  previously  described  (Schulz  et  al.,  2006a;  Tobin  et  al.,  2009).  Primers 
specific  for  real-time  PCR  detection  of  shal ,  BKKCa,  shah,  and  shaker 
using  Sybr  Green  were  developed  and  designed  using  Primer3  software 
and  are  the  same  as  previously  reported  (Schulz  et  al.,  2006a;  Tobin  et  al., 
2009).  Briefly,  total  RNA  was  isolated  using  RNeasy  micro  column-based 
RNA  extraction  kit  (Qiagen),  reverse  transcribed  using  Superscript  III 
reverse  transcriptase  (Invitrogen),  and  used  as  a  template  in  real-time 
RT -PCR  with  Sybr  Green  ( SABiosciences)  in  a  RotorGene  3000  real-time 
PCR  machine  (Corbett  Research).  Previous  studies  have  determined  that 
in  LC  motor  neurons,  correlations  can  be  equally  well  detected  among 
channel  mRNA  levels  with  and  without  normalization  of  real-time  re¬ 
sults  to  18S  rRNA  (Tobin  et  al.,  2009).  Values  reported  here  are  total  copy 
numbers  from  a  single  neuron,  and  are  not  normalized  with  respect  to 
18S  levels. 

Pharmacology.  Pharmacological  blockers  were  dissolved  in  physiolog¬ 
ical  saline  and  perfused  onto  the  cardiac  ganglion  using  a  Rabbit  peristal¬ 
tic  pump  (Rainin  Instruments)  at  a  rate  of  1.5  ml/min  or  added  to  the 
preparation  from  a  stock  solution  via  pipette.  The  following  pharmaco¬ 
logical  agents  were  used:  tetraethylammonium  dissolved  in  saline  at  25 
mM,  4-aminopyridine  dissolved  in  saline  at  1  mM,  cadmium  chloride 
dissolved  in  saline  at  250  /cm  (Acros  Organics),  tetrodotoxin  dissolved  in 
saline  at  1  /cm  (Alomone  Laboratories),  BAPTA-AM  dissolved  DMSO 
and  applied  at  30  /cm  in  saline,  ryanodine  dissolved  DMSO  and  applied  at 
100  /cm  in  saline,  staurosporine  dissolved  in  DMSO  and  applied  at  5  /cm 
in  saline,  okadaic  acid  dissolved  in  DMSO  and  applied  at  500  nM  in  saline 
(Ascent  Scientific),  cyclosporine  A  dissolved  in  DMSO  and  applied  at  2 
/cm  in  saline  (Tocris  Biosciences).  All  DMSO  applications  resulted  in  a 
final  concentration  of  DMSO  that  was  <1%  (range:  0.000025%  to  0.5%). 

Pharmacological  agents  used  to  investigate  intracellular  mechanisms 
involved  in  the  compensatory  response  (BAPTA-AM,  ryanodine,  oka¬ 
daic  acid,  cyclosporine  A,  staurosporine)  were  applied  to  the  cardiac 
ganglion  1  h  prior  (2  h  prior  with  ryanodine)  to  the  application  of  the 
blocker  which  caused  the  compensation  (TEA  or  4AP).  A  cell  or  prepa¬ 
ration  was  exposed  only  to  one  channel  blocker  type  (TEA  or  4AP)  for  a 
given  experiment. 

Current  measurements.  All  experiments  were  performed  in  physiolog¬ 
ical  saline  cooled  to  12°C.  To  measure  current  magnitudes  and  activation 
properties,  two-electrode  voltage-clamp  (TEVC)  experiments  were  per¬ 
formed  by  impaling  a  large  cell  with  two  glass  electrodes  filled  with  3  M 
KC1  (8-17  Mil  resistance)  and  an  Axoclamp  2A  amplifier  (Molecular 
Devices).  All  recordings  were  made  from  anterior  large  cell  somata;  ac¬ 
tion  potential  conductances  were  blocked  (unless  noted  otherwise)  by 
tightening  thread  ligatures  on  both  sides  of  the  large  cell  soma,  preserving 
space  clamp.  TEVC  protocols  were  created,  driven  and  recorded  with 
clampex  9.2  software  (Molecular  Devices).  Current  recordings  were  an¬ 
alyzed  with  Clampfit  9.2  software  (Molecular  Devices).  Current  and  volt¬ 
age  traces  were  sometimes  filtered  with  a  lowpass  boxcar  filter  using  7 


smoothing  points.  Most  voltage  clamps  were  modified  from  those  used 
previously  in  STG  preparations  (Golowasch  and  Marder,  1992; 
Khorkova  and  Golowasch,  2007;  Temporal  et  al.,  2012).  High  threshold 
potassium  current  (IHTK)  magnitude  was  measured  using  a  leak  sub¬ 
tracted  TEVC  protocol  with  a  holding  potential  of  —40  mV  and  16  volt¬ 
age  steps  from  —55  mV  to  +20  mV  (5  mV  intervals).  A- type  potassium 
current  (IA)  magnitude  was  measured  by  subtracting  the  IHtk  current 
traces  from  a  TEVC  protocol  that  is  identical  except  for  a  holding  poten¬ 
tial  of —80  mV.  Calcium-activated  potassium  current  (IKCa)  was  isolated 
by  subtracting  postcadmium  (250  /cm  CdCl2,  1  h)  IHTK  current  traces 
from  precadmium  IHTK  current  traces  (isolating  the  cadmium-sensitive 
outward  current).  Delayed  rectifier  potassium  current  (IKd)  was  isolated 
using  the  IHXK  TEVC  protocol  after  cadmium  exposure  (250  /jlm  CdCl2, 
1  h).  All  current  magnitude  measurements  were  taken  at  0  mV  on  an  I-V 
plot  made  from  the  current  traces. 

Large  cell  excitability  and  network  output.  Cardiac  network  output  was 
monitored  with  a  single  intracellular  recording  (using  same  equipment 
as  TEV  C  protocols)  taken  from  one  of  the  three  anterior  large  cells  and  an 
extracellular  differential  recording  made  with  a  model  1700  A-M  Systems 
AC  amplifier  and  two  stainless  steel  wires;  one  placed  inside  and  one 
outside  a  vaseline  well  located  around  the  central  nerve  of  the  CG  (see 
Fig.  1).  Pharmacological  blockers  (TEA  or  4AP)  were  perfused  on  the 
entirety  of  the  CG  with  the  exception  of  the  four  small  cell  and  two 
posterior  large  cell  somata.  These  cells  were  isolated  from  the  perfusion 
by  a  vaseline  well  placed  around  these  cells  and  the  posterior  branch 
point  containing  regular  physiological  saline  (see  Fig.  1).  Using  this  ex¬ 
perimental  configuration  we  monitored  the  effect  of  the  LC  compensa¬ 
tory  response  (isolated  to  the  best  of  our  ability  from  small  cells)  in  the 
context  of  the  functioning  network.  Network  activity  was  recorded  in 
10  min  intervals  before  and  during  blocker  perfusion.  These  record¬ 
ings  were  analyzed  using  Spike2  v6.00  software  (Cambridge  Elec¬ 
tronic  Design). 

LC  intrinsic  excitability  was  examined  under  similar  conditions  except 
10  6  M  TTX  saline  was  placed  in  the  vaseline  well  around  the  four  small 
cell  and  two  posterior  LC  somata.  This  eliminated  spontaneous  network 
activity  and  small  cell  excitatory  input  into  the  anterior  LCs.  Excitability 
in  the  anterior  LCs  was  then  monitored  using  two-electrode  current- 
clamp  (TECC)  protocols  run  before  and  every  5  min  after  TEA  or  4AP 
perfusion.  TECC  protocol  was  a  six  step  depolarizing  current  injection 
(from  1  to  6  nA)  lasting  six  seconds  per  step  and  six  seconds  between 
steps. 

Statistics.  All  statistical  tests  were  performed  with  SigmaPlot  vll.O 
(Systat,  Aspire  Software  International).  All  data  were  confirmed  to  be  of 
normal  distribution  as  required  by  statistical  analyses  used.  Relationships 
between  channel  mRNAs  and  ionic  currents  were  analyzed  using  Pear¬ 
son’s  correlation  test,  and  coefficients  of  determination  were  calculated 
from  the  resulting  correlation  coefficients.  In  the  case  of  Figure  2,  C  and 
D,  a  potential  outlier  was  identified  that  could  be  anchoring  a  false  pos¬ 
itive  for  the  Pearson’s  test  (see  arrows).  Analyses  on  these  datasets  were 
performed  both  with  and  without  the  data  point  in  question  and  both 
results  reported.  Bonferroni  corrections  were  used  for  multiple  compar¬ 
isons  in  the  correlation  analyses,  and  the  p  value  adjusted  to  0.017  for 
statistical  significance  (three  comparisons  each  for  mRNA  and  ionic  cur¬ 
rents).  Changes  in  current  magnitude  before  and  after  pharmacological 
block  were  analyzed  in  one  of  two  ways.  Raw  currents  were  analyzed 
before  and  after  blockade  via  paired  t  tests,  and  these  are  reported  in 
Figure  3.  In  Figure  5,  significant  changes  in  a  current  relative  to  baseline 
were  expressed  as  a  percent  change  from  zero,  and  analyzed  via  one- 
sample  t  test  with  the  hypothesized  population  mean  set  to  0.  Overall 
changes  in  burst  duration,  spikes  per  burst,  and  spike  frequency  within 
the  burst  reported  in  Figure  4  were  analyzed  with  repeated-measures 
ANOVA. 

Results 

The  crustacean  cardiac  ganglion  as  a  model  for  central 
pattern  generator  network  activity 

The  rhythmic  pumping  of  the  heart  in  decapod  crustaceans  such 
as  the  crab,  Cancer  borealis  (our  model  organism),  is  neurogenic 
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Figure  1.  The  crustacean  cardiac  ganglion.  When  dissected  from  the  animal,  the  cardiac  ganglion  can  be  pinned  flat  and 
simultaneous  recordings  of  network  activity  and  intracellular  recordings  can  be  obtained.  Five  large  cell  somata  (dark  ovals  marked 
with  arrowheads)  are  distributed  throughout  the  ganglion.  Individual  or  pairs  of  large  cells  can  be  isolated  for  voltage  clamp  with 
thread  ligature,  or  for  pharmacological  treatment  with  a  vaseline  well.  The  recording  shows  a  simultaneous  extracellular  and 
intracellular  recording  of  network  activity.  The  bursting  of  the  network  is  seen  at  the  site  indicated  by  the  dotted  circle  and 
corresponds  to  the  extracellular  (bottom)  trace.  Burst  output  of  a  single  LC  motor  neuron  is  seen  in  the  top  trace,  taken  from  an 
intracellular  recording  of  a  LC  soma.  This  rhythmic  motor  activity  corresponds  directly  to  the  contraction  of  the  single-chambered 
crab  heart  muscle,  showing  an  in  vitro  maintenance  of  biological  network  activity.  The  gray  oval  represents  a  vaseline  well  used  in 
pharmacological  blocker  experiments.  During  measurements  of  individual  LC  excitability  (Fig.  4  A),  this  well  was  used  to  apply  TTX 
to  the  small  cell  pacemaker  neurons  (SCs,  small  gray  circles)  to  silence  the  network  activity  but  preserve  output  capability  in 
anterior  LCs.  During  measurements  of  compensation  during  ongoing  network  activity,  this  cell  was  used  to  protect  SCs  from  blocker 
application  in  the  bath  (Fig.  46). 


in  nature,  and  under  the  control  of  a  simple  central  pattern  gen¬ 
erator  network  called  the  cardiac  ganglion  (Alexandrowicz,  1932) 
(Fig.  1).  The  ganglion  consists  of  only  nine  neurons:  four  “small 
cell”  interneurons  (SCs)  that  generate  the  pacemaker  activity  and 
five  “large  cell”  motor  neurons  (LCs)  that  innervate  the  heart 
musculature  (Hartline,  1967;  Tazaki  and  Cooke,  1983c).  The  SCs 
of  the  cardiac  ganglion  are  endogenous  oscillators,  i.e.,  they  un¬ 
dergo  spontaneous  and  rhythmic  generation  of  a  depolarizing 
wave  of  membrane  potential  that  leads  to  a  bursting  phenotype  of 
multiple  spikes  per  burst  (Tazaki  and  Cooke,  1983b;  Cooke, 
2002).  The  LCs  of  the  CG  produce  bursts  of  action  potentials  as  a 
result  of  synaptic  pacemaker  input  from  the  SCs,  ultimately  lead¬ 
ing  to  muscle  contraction  (Hartline,  1967;  Tazaki  and  Cooke, 
1979,  1983a, c;  Berlind,  1989).  The  behavioral  output  of  the  gan¬ 
glion  represents  a  direct  correlation  of  the  influence  of  LCs  on 
heart  muscle  (Sakurai  and  Wilkens,  2003;  Garcia-Crescioni  et  al., 
2010),  and  thus  a  direct  measure  of  heart  activity.  Yet  the  entire 
network  can  be  dissected  out  intact,  and  maintained  in  physio¬ 
logical  saline  for  extended  periods  of  recording  while  continuing 
to  produce  its  endogenous  rhythmic  output.  The  motor  neurons 
are  all  individually  identifiable,  and  due  to  their  distributed  na¬ 
ture  within  the  ganglion  (Fig.  1),  we  can  perform  pharmacolog¬ 
ical  manipulations  on  one  or  multiple  motor  neurons,  either  in 
isolation  or  in  the  intact,  functional  network  (Tazaki  and  Cooke, 
1983a;  Cooke,  2002).  The  underlying  burst  potentials  of  the  LCs 
represent  functional  output  at  the  motor  neuron  level,  and  our 
preliminary  modeling  studies  show  how  the  simplicity  of  this 
model  system  can  be  used  to  study  functional  implications  of  the 
relationships  between  mRNAs  and  ionic  conductances  on  cellu¬ 
lar  output  (Ball  et  al.,  2010;  Franklin  et  al.,  2010). 

Relationship  between  IA  and  IKCa  in  a  population  of  large  cell 
motor  neurons 

Potassium  currents  were  measured  using  protocols  developed  for 
the  STG  cells  of  the  same  species  (Golowasch  and  Marder,  1992). 


The  total  outward  current  of  LC  motor 
neurons  consists  primarily  of  three  K  + 
currents  (Fig.  2 B)  (Golowasch  and  Marder, 
1992):  A-type  transient  K+  current  (IA), 
calcium-activated  K+  current  (IKCa),  and 
delayed  rectifier  K+  current  (IKd).  IKCa 
and  IKd  can  be  found  in  one  combined 
current  trace  elicited  from  holding  poten¬ 
tials  at  —40  mV  or  higher,  and  is  termed 
here  as  the  high-threshold  K+  current 
(IHTk)>  while  the  A-type  current  can  be 
measured  by  subtracting  IHXK  from  the 
total  outward  current  elicited  from  a 
holding  potential  of  —80  mV. 

We  first  examined  the  relationships 
among  mRNA  levels  in  single  identified 
LC  motor  neurons  for  three  channel  genes 
that  correspond  to  these  three  K+  cur¬ 
rents  (Atkinson  et  al.,  1991;  Tsunoda  and 
Salkoff,  1995;  Kim  et  al.,  1998):  BKKCa 
(Ircs)*  shaker  (IA),  and  shab  (IKd).  We  de¬ 
tected  a  significant  correlation  between 
BKKCa  and  shaker  mRNA  levels  across  a 
population  of  20  LC  motor  neurons  (Fig. 
2 A,  left),  but  no  correlations  among  any 
other  channel  mRNAs  (Fig.  2A,  middle, 
right).  A  fourth  channel  mRNA  encoding 
an  A-type  K  +  current,  shal,  also  was  sig¬ 
nificantly  positively  correlated  to  BKKCa 
mRNA  levels  ( p  <  0.005;  R2  =  0.58),  as  well  as  to  shaker  ( p  < 
0.002;  R 2  =  0.42)  but  not  to  shab,  suggesting  an  overall  relation¬ 
ship  between  IKCa  and  IA,  but  not  with  IKd,  in  these  cells. 

We  next  examined  the  relationships  among  the  ionic  currents 
encoded  by  these  channel  genes  across  a  population  of  LC  motor 
neurons.  IA,  IKCa,  and  IKd  were  all  measured  in  each  cell  across  a 
population  of  LC  motor  neurons.  A  similar  pattern  of  correlated 
current  levels  was  seen  as  with  the  mRNA  with  one  striking  dis¬ 
tinction:  only  IA  and  IKCa  showed  a  significant  correlation  (Fig. 
2  A),  but  the  correlation  was  strongly  negative  as  opposed  to  pos¬ 
itive  as  seen  in  the  mRNA  measurements.  No  significant  correla¬ 
tions  were  found  between  IA  and  IKd,  or  IKCa  and  IKd  (Fig.  2A). 
Because  of  the  striking  difference  between  the  mRNA  correlation 
(positive),  and  the  current  correlation  (negative),  we  were  most 
interested  in  pursuing  the  relationship  between  IA  and  IKCa,  but 
had  concerns  about  the  effects  of  pharmacological  blockers  used 
to  measure  IKCa  (Fig.  3 F);  we  use  Cd2+  to  block  calcium  currents 
that  evoke  IKCa  to  measure  this  current  via  subtraction  from  IHtk> 
as  we  have  not  identified  any  blockers  specific  to  IKCa  in  our 
preparation  despite  many  attempts.  Therefore,  we  decided  to 
use  the  peak  IHtk  as  an  indicator  of  IKCa  abundance  in  these 
experiments. 

IHTk  (HighThreshold  K  +  )  is  known  from  previous  work 
(Golowasch  and  Marder,  1992;  Haedo  and  Golowasch,  2006; 
Khorkova  and  Golowasch,  2007)  to,  in  large  part,  consist  of  IKCa, 
particularly  the  transient  portion,  as  IKd  shows  no  transient  peak 
that  could  account  for  the  peak  in  IHXK  (Fig.  2  B).  We  also  deter¬ 
mined  that  measurements  of  the  peak  IHXK  were  likely  sufficient 
to  reveal  the  relationship  between  IKCa  and  IA  as  the  only  tran¬ 
sient  peak  seen  in  the  currents  that  underlie  IHXK  belongs  to  IKCa 
(Fig.  2 B,  right).  Peak  IHXK  also  shows  the  same  negative  relation¬ 
ship  with  IA  as  does  IKCa  (Fig.  2C,  left),  and  peak  IKCa  itself  very 
strongly  correlates  with  peak  IHXK  (Fig.  2C,  middle),  but  not  IKd 
(Fig.  2C,  right).  Accordingly,  for  the  remainder  of  the  study,  we 
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used  IHTK  peak  current  from  baseline  as  a 
measure  of  IKCa  abundance,  to  allow  for 
measurement  of  all  K  +  currents  under  the 
least  manipulative  pharmacological  blocker 
regime,  i.e.,  without  the  need  to  block 
voltage-gated  calcium  currents. 
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Figure  2.  Relationships  among  channel  gene  mRNA  levels  and  among  ionic  currents  in  a  normal  population  of  large  cell  motor 
neurons.  A,  Correlations  between  mRNA  levels  from  single  identified  LC  motor  neurons  for  BKKCa,  shaker,  and  shab  channel  genes 
(top  row)  and  the  corresponding  K  +  current  correlations  from  single  identified  LC  motor  neuronsfor  lKca,  lfl,  and  lKd  (bottom  row). 
Each  point  represents  values  from  a  single  neuron.  Regression  lines  shown  only  for  significant  correlations  as  revealed  by  Pearson's 
test,  for  which  thep  value  and  R2  values  are  reported  in  each  plot.  When  two  regression  lines  are  present,  the  solid  line  and  top  set 
of  statistics  refers  to  the  entire  dataset,  while  the  dotted  line  and  the  bottom  set  of  statistics  represents  the  correlation  with  the 
single  point  indicated  by  the  arrow  removed.  Despite  having  this  outlier  removed,  the  correlation  is  still  significant  and  largely 
unchanged.  fi,  The  total  K  +  current  (lK(Tot))  in  LC  motor  neurons  consists  primarily  of  three  distinct  K  +  currents:  A-type  K  +  (lA), 
and  the  high-threshold  K  +  current  (lHTK)  that  can  be  subdivided  into  calcium-activated  K  +  (lK[a)  and  delayed  rectifier  K  +  (lKd).  C, 
Correlations  between  currents  reveal  that  lHTK  peak  is  highly  representative  of  lKCa  levels,  but  not  of  lKd.  With  these  results,  we  use 
lHTK  peak  levels  as  an  indicator  of  lKCa  for  any  given  cell  for  the  remainder  of  the  study.  Statistics  as  in  A. 


IKCa  and  IA  are  rapidly  upregulated 
within  the  same  LC  motor  neuron 

Because  the  focus  of  the  study  was  on  the 
negative  relationship  between  IA  and  IKCa, 
we  decided  to  repeat  the  initial  measure¬ 
ments  in  a  new  population  of  LC  motor 
neurons  to  confirm  the  original  finding 
before  examining  the  functional  impact  of 
this  relationship.  In  a  distinct  population 
of  20  LC  motor  neurons,  we  once  again 
found  a  strong  negative  relationship  be¬ 
tween  IA  and  IKCa  (as  revealed  by  IHXK 
measurements;  Fig.  3 A),  confirming  our 
original  finding.  This  strongly  suggested  a 
causal  relationship  between  IA  and  IKCa. 

Specifically,  we  hypothesized  that  if  these 
two  currents  are  functionally  interrelated, 
then  a  decreased  level  of  IA  should  result 
in  an  increased  level  of  IKCa  and  vice  versa. 

We  tested  this  hypothesis  by  blocking  one 
current  in  this  pair  and  determining 
whether  there  was  an  effect  on  the  magni¬ 
tude  of  the  other.  We  measured  baseline 
levels  of  IA  in  a  given  cell,  then  used  tetra- 
ethylammonium  (TEA)  to  block  the 
HTK-current  for  60  min,  and  then  mea¬ 
sured  IA  again  post-TEA.  60  min  of  expo¬ 
sure  to  TEA  significantly  increased  peak 
( p  <  0.001;  paired  t  test)  A-current  (Fig. 

3  B,C).  The  converse  experiment  was  per¬ 
formed  using  4-aminopyridine  (4AP)  for 
60  min  to  block  IA.  The  peak  (p  <  0.003)  of  IHTK  was  significantly 
increased  following  60  min  of  4-AP  blockade  (Fig.  3 B,C).  Addi¬ 
tionally,  these  results  appear  to  be  due  to  an  overall  increase  in 
conductance,  since  no  change  was  seen  in  the  voltages  of  activa¬ 
tion  of  IA  and  IHXK  concomitant  with  the  changes  in  total  current 
as  a  result  of  the  block  experiments  (Fig.  3D).  No  consistent  or 
significant  effects  of  TEA  and  4-AP  on  neuronal  input  resistance 
were  observed  in  these  cells  (A  RIN  TEA:  — 1.87  ±  2.5  Mfl,  A  R1N 
4-AP:  -0.52  ±  1.7  Mil). 

Our  initial  data  strongly  implicate  IKCa  and  IA  as  the  key  pair  of 
currents  involved  in  this  response.  Since  IHtk  is  a  mixed  current 
comprised  predominantly  of  IKCa  and  IKd,  and  both  of  these  cur¬ 
rents  are  blocked  by  TEA,  we  determined  whether  the  changes 
seen  in  IHXK  as  a  result  of  4AP  blockade  were  attributable  to  just 
one  or  to  both  of  these  currents.  Most  convincingly,  while  4AP 
block  significantly  increases  IHXK  (Fig.  3C),  1  h  of  exposure  to 
4AP  had  no  significant  effect  on  IKd  magnitude  itself  (Fig.  3 E, 
left).  Unfortunately,  we  have  not  identified  a  blocker  specific  to 
IKCa  in  our  preparation  to  directly  test  the  effects  of  blocking  this 
current  on  changes  in  IA.  However,  we  can  indirectly  block  IKCa 
by  using  CdCl2  to  block  voltage-gated  Ca2+  channels  that  trigger 
IKCa.  Blocking  with  Cd2+  completely  eliminates  IKCa  in  these 
cells,  and  causes  a  significant  increase  in  IA  after  1  h  of  exposure 
that  is  very  similar  to  that  seen  with  TEA  (Fig.  3 E,  right).  To¬ 
gether,  these  data  strongly  suggest  that  the  relationship  between 


IA  and  IHXK  identified  in  the  experiments  (Fig.  3)  are  the  result  of 
a  causal  relationship  between  the  levels  of  IA  and  IKCa  in  a  given 
motor  neuron. 

IA  and  IKCa  act  in  a  compensatory  fashion  to  stabilize  cellular 
excitability  and  network  output 

The  negative  correlation  we  measured  between  IA  and  IKCa  in  the 
normal  population,  together  with  the  rapid  change  in  currents 
seen  within  a  cell  after  TEA  and  4AP  block,  led  us  to  hypothesize 
that  levels  of  IA  and  IKCa  may  be  acting  in  a  compensatory  fashion 
to  stabilize  both  cellular  and  network  output.  To  determine 
whether  such  a  functional  compensation  exists,  we  followed  over 
time  the  activity  of  both  isolated  LCs,  as  well  as  LCs  in  an  intact 
network,  with  either  4AP  or  TEA  blockade,  and  then  determined 
whether  the  output  of  the  cells  and  of  the  network  changed  over 
the  course  of  the  blockade  and  the  potential  compensation. 

As  seen  in  Figure  4,  A  and  B,  TEA  blockade  caused  a  substan¬ 
tial  increase  in  the  excitability  of  the  LC  motor  neuron  in  the 
initial  10-20  min.  However,  over  the  time  course  of  the  measured 
increase  in  IA  seen  in  previous  experiments  (Fig.  2C)  there  was  a 
compensatory  change  in  the  output  of  the  cell  and  of  the  network 
toward  the  baseline  level  of  activity.  In  isolated  LCs  (Fig.  4A), 
both  TEA  and  4AP  application  shifted  the  cells  from  a  less  excit¬ 
able  state  to  a  state  characterized  by  large  sustained  burst  poten¬ 
tials.  However,  over  the  time  course  of  compensation,  the 
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Figure  3.  Reciprocal  interaction  of  lA  and  lHTK  levels  in  LC  motor  neurons.  A,  Correlation  between  peak  lA  and  lHTK  across  a  second  population  of  LC  motor  neurons  distinct  from  those  in  Figure  2. 
Statistics  represent  results  of  a  Pearson's  test  (n  =  20).  B,  Representative  recordings  of  the  increase  in  lA  after  TEA  exposure  (left)  and  the  increase  in  lHTK  after  4AP  exposure  (right).  Black  trace  is  the 
control,  gray  trace  is  the  same  cell  after  60  min  exposure  to  blocker.  C,  Quantified  effects  of  TEA  (left)  and  4AP  (right)  on  the  peak  lA  and  lHTK  respectively  in  N  =  6  cells.  The  same  cell  before  (black) 
and  after  (white)  is  connected  with  a  solid  line.  Statistics  represent  paired  (test  on  mean  current  before  and  after  blocker.  D,  Activation  characteristics,  expressed  as  l/lmal(  for  an  N  =  5  cells,  of  lfl  (left) 
and  lHTK  (right)  before  and  after  block  with  TEA  and  4AP,  respectively.  Each  point  represents  mean  ±  SD.f,  Effects  of  4AP  on  lKd  (left;  A/  =  5)andCd2+  on  lA  (right;  N  =  6).  Insert  in  the  left  panel 
shows  representative  traces  of  lKd  before  (black)  and  after  (gray)  exposure  to  4AP.  Legends  and  statistics  as  in  C. 


excitability  of  isolated  cells  returned  to  control  levels,  losing  the 
large  burst  potentials  and  showing  firing  patterns  similar  to  those 
before  the  application  of  the  blocker  (Fig.  4A). 

In  the  intact  network,  TEA  often  caused  a  complex  multiphase 
bursting  output  in  treated  LCs  (Fig.  4B;  10  min  recording).  Given 
this  somewhat  complex  pattern  of  bursting,  we  determined  that 
measurements  of  burst  duration  across  a  population  of  these  cells 
represented  the  clearest  means  to  quantify  changes  in  output  over¬ 
time.  Over  the  course  of  60-90  min,  TEA  blockade  resulted  in  a 
significant  change  in  burst  duration  ( p  <0.01;  ANOVA;  Fig.  4C, 
left),  initially  increasing  to  nearly  twice  as  long  as  control  before 
re-establishing  a  stable  burst  duration  indistinguishable  from 
that  at  the  control  level  (Fig.  4B). 

In  the  companion  experiment  with  4AP  blockade  in  the  intact 
network,  we  see  a  similar  effect  of  compensation  as  in  the  LC 
motor  neuron  output  (seen  in  Fig.  4B,C).  4AP  block  does  not 


cause  the  same  change  in  burst  duration  as  with  TEA.  Rather,  4AP 
blockade  initially  causes  a  significant  increase  in  spike  frequency 
and  the  number  of  spikes  fired  per  burst  ( p  <  0.01;  ANOVA). 
Over  the  time  course  of  the  blockade,  this  spike  frequency  and 
number  of  spikes  per  burst  returns  to  control  levels  (Fig.  4 D). 
These  data  indicate  that  the  changes  seen  in  the  level  of  IA  and 
IKCa  as  a  result  of  TEA  and  4AP  block,  respectively,  are  a  func¬ 
tional  compensatory  response  that  stabilizes  the  output  of  the 
cells  and  of  the  network  as  a  whole.  However,  the  roles  of  IA  and 
IKCa  are  n°t  simply  functionally  redundant;  blocking  either  cur¬ 
rent  causes  an  initial  change  in  the  output  of  the  cell  that  is  not  the 
same.  These  results  implicate  distinct  roles  for  these  two  currents 
in  addition  to  their  ability  to  partially  compensate  for  one  an¬ 
other.  However,  we  do  not  rule  out  the  possibility  that  there 
could  be  other  currents  or  mechanisms  involved  in  the  compen¬ 
satory  effect. 
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Compensatory  changes  in  IA  and  IKCa 
follow  distinct  regulatory  pathways 

We  then  set  out  to  determine  the  under¬ 
lying  mechanisms  implementing  the 
compensation  in  these  currents.  We  first  in¬ 
vestigated  whether  the  change  in  IKCa  elic¬ 
ited  by  4AP  blockade  was  dependent  on 
intracellular  calcium.  To  check  this,  we  co¬ 
applied  the  calcium  chelator  BAPTA  with 
the  blocker.  When  BAPTA  was  co-applied 
with  4AP,  there  was  no  change  in  IHtk>  in 
contrast  to  the  characteristically  significant 
increase  in  IHTK  with  4AP  blockade  (Fig. 

5A).  We  also  confirmed  that  this  effect  was 
not  simply  due  to  an  effect  of  BAPTA  on  the 
calcium-dependence  of  IKCa;  BAPTA  alone 
did  not  result  in  a  significant  reduction  of 
IHtk  (data  not  shown;  see  also  Turrigiano  et 
al.,  1994).  More  specifically,  the  calcium  de¬ 
pendence  of  the  increase  in  IHTK  was  attrib¬ 
utable  to  the  release  from  intracellular 
calcium  stores  since  co-application  of  4AP 
and  the  intracellular  calcium  release  blocker 
ryanodine  also  prevented  the  compensatory 
increase  in  IHTK  normally  seen  with  4AP 
blockade  (Fig.  5A).  These  results  suggest 
that  the  compensatory  increase  in  IKCa  is  de¬ 
pendent  on  intracellular  calcium  signaling 
mechanisms  that  depend  on  the  release  of 
intracellular  calcium  stores. 

Because  relatively  rapid  changes  in 
current  magnitude  can  often  be  attributed 
to  changes  in  phosphorylation  states  of 
ion  channels  and  associated  regulatory 
proteins,  we  examined  the  effects  of  phos¬ 
phatase  and  kinase  inhibitors  on  the  ability 
of  the  blockers  to  elicit  the  compensatory 
changes  in  these  currents.  Okadaic  acid  is 
a  broad  spectrum  inhibitor  of  serine/thre¬ 
onine  protein  phosphatases  (Cohen  et  al., 

1990).  While  application  of  okadaic  acid 
alone  did  not  affect  baseline  levels  of  IHTK 
(data  not  shown),  co-application  of  oka¬ 
daic  acid  with  4AP  also  abolished  the 
compensatory  increase  in  IHTK  (Fig.  5A), 
suggesting  the  action  of  a  serine/threo¬ 
nine  protein  phosphatase  in  this  compen¬ 
satory  effect.  Because  calcineurin  is  a  well 
known  Ca2+-dependent  protein  serine/ 
threonine  phosphatase  (Klee  et  al.,  1998), 
we  used  cyclosporine  A  to  inhibit  cal¬ 
cineurin  activity  in  our  system.  Cyclosporine  alone  did  not  affect 
baseline  levels  of  IHTK  (data  not  shown),  but  combined  applica¬ 
tion  of  4AP  with  cyclosporine  eliminated  the  compensatory  in¬ 
crease  in  IHTK  (Fig.  5A),  indicating  that  the  increase  in  IHXK  is 
dependent,  at  least  in  part,  on  the  activity  of  calcineurin. 

In  complete  contrast,  the  compensatory  influence  of  TEA  on 
IA  appears  to  be  implemented  via  mechanistically  distinct  path¬ 
ways  compared  with  those  we  found  for  the  4AP  block  on  IHTK. 
The  increase  of  IA  in  response  to  TEA  blockade  is  neither 
calcium-dependent  (Fig.  5 B)  nor  dependent  on  the  activity  of 
serine/threonine  phosphatases  (Fig.  5B).  Parallel  experiments 
showed  that  BAPTA,  ryanodine,  okadaic  acid  and  cyclosporine, 
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Figure  4. 

Representative  recordings  of  the  effects  on  the  excitability  of  individual  LC  motor  neurons  isolated  from  network  activity  after 
exposure  to  TEA  (left)  and  4AP  (right).  Because  the  pacemaker  neurons  of  the  network  were  silenced  with  TTX  in  this  experiment, 
DC  current  injections  were  used  to  elicit  activity  in  LC  motor  neurons  (current  traces  shown  at  top).  8,  Representative  recordings  of 
the  effects  of  TEA  (left)  and  4AP  (right)  exposure  on  network  output  overtime.  Pacemaker  cells  were  shielded  from  exposure  to  the 
blocker  (see  Fig.  1),  while  the  3  upper  motor  neurons  were  exposed  to  blocker.  Paired  intracellular  (top  traces)  and  extracellular 
recordings  were  made  over  the  time  course  of  the  experiment.  C,  Quantified  effects  of  TEA  and  4AP  effects  on  burst  characteristics 
of  LC  motor  neurons  in  the  intact  network.  Bars  represent  mean  ±  SD  of  A/  =  5  preparations  at  10  min  intervals.  The  first  10  min 
interval  is  baseline  without  blocker,  and  subsequent  bins  represent  cumulative  time  after  blocker  exposure. 


when  co-applied  with  TEA,  failed  to  prevent  a  significant  increase 
in  IA  (Fig.  5 B).  The  only  exception  was  that  the  application  of 
ryanodine  in  conjunction  with  TEA  appeared  to  reduce  the  mag¬ 
nitude  of  increase  in  IA,  although  this  effect  was  not  statistically 
significant  with  respect  to  TEA  alone  ( p  =  0.103;  t  test),  or  TEA 
plus  BAPTA  ( p  =  0.08). 

We  also  examined  the  effects  of  inhibiting  kinase  activity  on 
the  compensatory  effects  of  4AP  and  TEA  blockade.  Staurospor- 
ine  is  a  potent,  cell-permeable  protein  kinase  C  inhibitor  which 
also  partially  inhibits  other  kinases  such  as  PKA,  PKG,  and  CaM- 
KII  (Riiegg  and  Burgess,  1989).  Co-application  of  staurosporine 
with  4AP  appears  to  cause  a  significant  decrease  in  IHTK  (Fig.  5A). 
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Figure  5.  Use  of  pharmacological  inhibitors  to  probe  mechanisms  of  changes  in  lA  and  lKCa. 
A,  Effects  of  BAPTA,  ryanodine,  okadaic  acid,  cyclosporine,  and  staurosporine  (referred  to  as  the 
treatment)  on  the  induction  of  increased  lHTK  with  4AP  exposure.  All  bars  are  mean  ±  SD, 
sample  sizes  as  shown  in  bars.  All  results  are  percentage  change  in  peak  lHTK  with  4AP  + 
treatment  relative  to  treatment  alone  (see  Materials  and  Methods).  None  of  the  treatments 
caused  a  change  in  baseline  levels  of  lHTK  with  the  exception  of  staurosporine  (dotted  box;  see 
£].***  represents  a  significant  (p  <  0.001 )  difference  from  0  via  one-sample  f  test.  8,  Effects  of 
BAPTA,  ryanodine,  okadaic  acid,  cyclosporine,  and  staurosporine  on  the  induction  of  increased 
lA  with  TEA  exposure.  All  results  are  percentage  change  in  peak  lA  with  TEA  +  treatment  relative 
to  treatment  alone  (see  Materials  and  Methods).  None  of  the  treatments  caused  a  change  in 
baseline  levels  of  lA.  ***p  <  0.001  and  *p  <  0.01  represent  a  significant  difference  from  0  via 
one-sample  f  test.  C,  Effects  of  staurosporine  alone  on  baseline  levels  of  lA  and  lHTK.  ***  repre¬ 
sents  a  significant  (p  <  0.001)  difference  from  0  via  one-sample  t  test.  D,  Overall  effects  of 
staurosporine  alone  and  then  4AP  +  staurosporine  on  lHTK.  Each  individual  cell  is  connected 
with  a  solid  line  across  the  treatments. 


dent  of  all  of  the  regulatory  pathways  implicated  in  the  IKCa  re¬ 
sponse.  These  effects  are  also  fairly  rapid,  acting  over  the  course  of 
60-90  min  to  stabilize  the  activity  of  a  critical  CPG  network 
responsible  for  cardiac  muscle  contraction  and  heart  beat  gener¬ 
ation  in  the  animal.  While  such  homeostatic  responses  and  their 
role  in  stabilizing  synaptic  function  have  been  well  studied 
(Bergquist  et  al.,  2010;  Turrigiano,  2011),  we  know  much  less 
about  the  mechanisms  underlying  homeostatic  plasticity  of  in¬ 
trinsic  excitability  and  the  role  this  form  of  plasticity  plays  in  the 
stabilization  of  neuronal  and  motor  network  output  (Turrigiano, 
2011). 

Even  neurons  and  networks  with  extremely  robust  output 
display  highly  variable  underlying  physiological  parameters  re¬ 
sponsible  for  neuronal  output,  particularly  membrane  conduc¬ 
tances  (Schulz  et  al.,  2006a,  2007;  Khorkova  and  Golowasch, 
2007;  Goaillard  et  al.,  2009;  Temporal  et  al.,  2012).  Our  work 
reveals  that  embedded  within  this  variability  are  coregulatory 
relationships  that  act  to  stabilize  the  excitability  of  the  cell,  in  part 
by  balancing  the  sum  total  of  major  transient  outward  currents, 
IA  and  IKCa.  A  similar  relationship  between  IA  and  IKCa  was  re¬ 
ported  in  the  STG.  Artificial  depolarization  of  inferior  cardiac 
neurons  results  in  increased  IA  and  decreased  IHXK,  which  is  abol¬ 
ished  by  Cd2+  blockade  of  calcium  channels  (Golowasch  et  al., 
1999).  These  experiments  did  not  reveal  a  naturally  existing  cor¬ 
relation  between  these  currents  (Golowasch  et  al.,  1999),  sug¬ 
gesting  that  excitability  in  the  STG  may  be  a  more  complex 
interaction  among  multiple  conductances,  or  perhaps  between 
conductances  and  their  constitutive  neuromodulation  (Harris- 
Warrick,  2011).  Although  it  is  not  known  what  effect  this  has  on 
network  output  in  the  STG,  our  complementary  results  suggest 
that  such  homeostatic  mechanisms  may  be  common  among  mo¬ 
tor  neurons  in  different  CPGs. 


However,  this  is  a  more  complex  result  than  that  for  the  other 
pharmacological  blockers.  Unlike  any  of  the  other  pharmacolog¬ 
ical  blockers  used  in  this  study,  staurosporine  was  the  only  one  to 
cause  a  change  in  baseline  levels  of  IHTK  when  applied  alone  (Fig. 
5C).  Therefore,  combined  application  of  4AP  and  staurosporine 
results  in  a  relative  decrease  in  IHTK  compared  to  that  with  stau¬ 
rosporine  alone,  but  the  net  effect  of  4AP  +  staurosporine  is  a 
restoration  of  IHXK  to  baseline  levels  (Fig.  5D). 

Once  again,  the  effects  of  staurosporine  on  IA  in  the  context  of 
TEA  blockade  contrast  with  those  for  4AP  and  IHTK.  There  is  no 
effect  of  staurosporine  on  baseline  levels  of  IA  (Fig.  5C),  and  while 
staurosporine  co-applied  with  TEA  results  in  an  apparent  de¬ 
crease  in  the  magnitude  of  the  effect  on  IA,  there  is  still  a  signifi¬ 
cant  increase  in  IA  with  staurosporine  present  (Fig.  5 B),  which  is 
not  statistically  significant  from  the  effect  of  TEA  alone  ( p  =  0.08; 
t  test). 

Discussion 

Rapid  compensation  between  IA  and  IKCa  preserves  both 
cellular  and  network  outputs 

We  have  identified  a  naturally  occurring  coregulatory  relation¬ 
ship  between  potassium  currents  (IA  and  IKCa)  in  an  intact  CPG 
network  that  results  in  homeostatic  compensation  of  neuronal 
excitability  as  well  as  network  function.  Furthermore,  we  have 
determined  that  these  compensatory  changes  in  K  +  current  mag¬ 
nitudes  are  independently  regulated  by  distinct  mechanisms. 
Compensatory  increases  in  IKCa  are  calcium-dependent  and  due, 
at  least  in  part,  to  the  activity  of  calcineurin-based  phosphatase 
activity.  Conversely,  compensatory  increases  in  IA  are  indepen¬ 


Distinct  intracellular  pathways  mediate  compensation 

The  mechanisms  involved  in  the  regulation  of  the  compensatory 
response  mediated  by  IKCa  after  4AP  block  in  the  CG  network  are 
consistent  with  work  supporting  changes  in  K  1  current  density 
in  a  homeostatic  fashion  (Desai  et  al.,  1999;  Schulz  et  al.,  2006b; 
Debanne  and  Poo,  2010;  Misonou,  2010).  In  particular,  the 
mechanisms  we  see  for  the  compensatory  increase  in  IKCa  are 
similar  to  those  found  in  regulation  of  excitability  in  cultured 
hippocampal  neurons  (Misonou  et  al.,  2004,  2005;  Misonou  and 
Trimmer,  2004).  Increased  excitability  in  these  cultured  neurons 
results  in  calcineurin-dependent  dephosphorylation  of  Kv2.1 
channels,  leading  to  a  functional  potentiation  of  these  channels 
via  a  shift  in  activation  voltage,  and  to  restoration  of  excitability 
(Misonou  et  al.,  2004).  While  our  data  strongly  suggest  that  IKCa 
is  primarily  responsible  for  the  observed  change  in  IHXK,  we  can¬ 
not  conclusively  rule  out  a  role  for  IKd  (similar  to  Kv2.1).  How¬ 
ever,  cyclosporine  has  been  shown  to  alter  IKCa  channels  directly 
(Hay,  1998),  and  calcineurin  is  known  to  have  widespread  effects 
on  neuronal  plasticity  and  excitability  beyond  its  effects  on  Kv2. 1 
activity  (Groth  et  al.,  2003).  The  effects  seen  in  our  experiments 
likely  can  be  attributed  to  an  overall  increase  in  the  maximal 
conductance  of  IKCa,  as  we  see  no  changes  in  either  activation 
curves  or  any  consistent  or  significant  changes  in  input  resis¬ 
tances  of  these  cells  (O’Leary  et  al.,  2010). 

The  mechanisms  underlying  compensatory  increases  in  IA  are 
less  clear,  but  distinct  from  the  pathways  involved  in  the  upregu- 
lation  of  IKCa.  Unlike  for  IKCa,  the  increase  in  IA  is  neither 
calcium-dependent,  nor  influenced  by  the  activity  of  calcineurin 
or  any  other  phosphatases  that  would  be  affected  by  treatment 
with  okadaic  acid  (Cohen  et  al.,  1990).  Indeed,  blocking  neither 
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Figure  6.  A  working  model  for  the  induction  of  homeostatic  plasticity  of  intrinsic  excitability  via  rapid  compensation  between  lA  and  lKCa.  This  model  summarizes  all  of  thedata  in  the  current  study 
and  provides  a  framework  for  interpretation  and  future  experimentation.  The  two  halves  of  the  cell  represent  the  distinct  mechanisms  induced  by  a  block  of  lA  with  4AP  (left,  blue  pathway)  and  a 
block  of  lKCa  with  TEA  (right,  green  pathway).  Blue  Pathway  (4AP  blockade):  The  baseline  activity  (1 )  is  disrupted  by  the  block  of  A-type  channels  with  4AP  (2),  causing  an  alteration  in  the  excitability 
of  the  cell  and  its  output  in  the  network,  in  this  case  rapid  initial  firing  frequency  (3).  This  change  presumably  results  in  the  release  of  calcium  from  intracellular  stores  (4)  that  increases  activity  of  the 
calcium-dependent  phosphatase  calcineurin  (5).  Calcineurin  activity  inducesan  increase  in  lKCa  (6),  restoring  the  excitability  of  the  cell  and  re-establishing  motor  neuron  activity  in  the  intact  network 
(7).  Green  Pathway  (TEA  blockade):  Baseline  network  activity  (1)  is  disrupted  via  TEA  blockade  (2),  causing  substantial  changes  in  the  output  properties  of  the  motor  neuron  and  the  network  (3). 
While  the  intracellular  mechanisms  remain  unclear  atthis  point,  a  subsequent  increase  in  lA  that  is  independent  of  the  pathway  invoked  by  4AP  blockade  compensatesfor  loss  of  lKCa,  restoring  output 
of  the  neurons  and  the  network  (5).  This  compensation  is  rapid  (on  the  order  of  1  h),  suggesting  thatthere  isan  already  existing  "reservoir"  or  channel  protein  available  that  is  maintained  via  positive 
coregulation  at  the  transcription/mRNA  level  of  the  channel  genes  responsible  for  lKCa  and  lA  (Gray  Oval). 


phosphatase  nor  kinase  activity  by  two  broad  spectrum  blockers 
was  able  to  prevent  an  increase  in  IA  as  a  result  of  TEA  blockade, 
possibly  precluding  the  contribution  of  phosphorylation  state  to 
this  half  of  the  compensation  story.  However,  we  were  able  to 
dampen  the  effect  by  blocking  release  of  intracellular  calcium 
with  ryanodine,  as  well  as  by  inhibiting  kinase  activity,  suggesting 
a  more  complex  path  of  regulation  in  this  response.  Given  the 
potential  nonspecific  effects  of  pharmacological  treatment,  we 
cannot  rule  out  at  least  some  role  for  a  kinase  in  this  pathway.  For 
example,  Kvl.2  potassium  channels  are  known  to  be  affected  by 
cAMP/protein  kinase  A  pathways  which  enhance  their  conduc¬ 
tance,  in  part  by  altering  trafficking  of  these  channels  (Connors  et 
al.,  2008).  Our  staurosporine  treatment  may  result  in  only  partial 


inhibition  of  PKA  (Riiegg  and  Burgess,  1989),  which  could  ex¬ 
plain  the  intermediate  effects  seen  on  the  A-type  current  com¬ 
pensation.  Regardless,  the  fact  that  IA  and  IKCa  are  regulated  by 
distinct  pathways  demonstrates  an  inherent  robustness  to  this 
homeostatic  response,  which  may  be  characteristic  of  networks 
responsible  for  critical  functions,  such  as  CPGs. 

Reservoirs  of  channel  proteins  might  implement 
rapid  compensation 

Our  data  also  begin  to  shed  light  on  a  striking  disparity  between 
levels  of  regulation  in  this  system.  Namely,  we  see  an  entirely 
different  relationship  at  the  level  of  mRNA  for  the  channel  genes 
BKKCA  and  shaker  (positive  correlation),  than  we  do  for  the 
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currents  IKCa  and  IA  that  these  channels  most  directly  encode 
(negative  correlation).  Indeed,  to  date  every  correlation  we  have 
found  between  channel  mRNAs  at  the  single  cell  level  has  been  a 
positive  correlation  (Schulz  et  al.,  2007;  Tobin  et  al.,  2009;  Tem¬ 
poral  et  al.,  2012).  These  data  include  measurements  of  at  least 
nine  different  channel  genes  in  seven  different  motor  neuron 
types  in  the  crustacean  stomatogastric  and  cardiac  ganglia  (S. 
Temporal  and  D.  J.  Schulz,  unpublished  data),  for  a  total  of  no 
fewer  than  60  correlations  in  channel  expression  detected  among 
channel  transcript  levels,  and  every  one  has  thus  far  been  positive. 
In  addition,  our  data  impact  the  present,  although  admittedly 
limited,  understanding  of  the  relationship  between  channel 
mRNA  and  ionic  current  in  a  given  neuron.  Namely,  that  mRNA 
for  a  given  K  +  channel  correlates  positively  with  its  correspond¬ 
ing  ionic  conductance  (Schulz  et  al.,  2006).  Although  both  of 
these  measurements  were  not  made  in  the  same  cells  in  this  study, 
logic  dictates  that  there  cannot  simply  be  a  positive  relationship 
between  mRNA  and  conductance  for  both  K  +  currents  in  this 
study.  Clearly,  and  not  particularly  surprisingly,  our  data  dem¬ 
onstrate  that  more  complexity  lies  within  the  relationship  of 
mRNA  and  ionic  conductance  than  a  one-for-one  tracking  be¬ 
tween  these  disparate  levels  of  function  in  all  cell  types.  However, 
our  data  continue  to  support  the  concept  that  correlation  among 
mRNA  levels  (Schulz  et  al.,  2007;  Tobin  et  al.,  2009)  underlies  a 
functional  relationship  at  the  conductance  level. 

So  then  why  do  we  see  only  positively  correlated  mRNA  levels 
for  channel  genes  in  these  motor  networks?  The  linkage  to  cellu¬ 
lar  and  network  output  in  our  study  enables  us  to  propose  a  new 
hypothesis  in  this  regard.  We  hypothesize  that  rapid  compensa¬ 
tion  of  intrinsic  excitability,  such  as  that  reported  in  this  study,  is 
implemented  via  an  already  existing  pool  of  protein.  That  is, 
protein  pools  for  the  two  opposing  channels  must  already  be 
available  in  relatively  equivalent  numbers,  even  if  they  are  not 
functionally  equally  represented.  This  requires  that  mechanisms 
must  occur  at  the  level  of  gene  regulation  to  ensure  a  sufficient 
“reservoir”  of  protein  capable  of  compensating  for  the  loss  of  a 
given  channel.  Thus,  regardless  of  whether  channels  act  in  con¬ 
cert  to  set  neuronal  output  (Ball  et  al.,  2010;  Franklin  et  al.,  2010), 
or  act  in  a  compensatory  fashion,  the  overall  gene  regulation 
between  the  two  must  be  relatively  balanced.  Such  a  reservoir 
perspective  is  reminiscent  of  mechanisms  seen  for  synaptic  vesi¬ 
cle  protein  dynamics  (Fernandez-Alfonso  et  al.,  2006),  or  extra- 
synaptic  AMPA  receptors  that  serve  as  a  reservoir  during  synaptic 
plasticity  (Hayashi  et  al.,  2000;  Zhu  et  al.,  2000).  Interestingly, 
insertion  of  AMPA  receptor  during  synaptic  plasticity  is  also  un¬ 
der  the  control  of  phosphorylation-dependent  processes  (Lin  et 
al.,  2009).  Therefore,  this  hypothesis  provides  a  compelling 
framework  for  future  investigation  of  the  relationship  between 
mRNA  levels  and  ionic  conductance,  as  well  as  the  interplay  be¬ 
tween  gene  expression  and  post-transcriptional  mechanisms  in¬ 
volved  in  homeostatic  plasticity  of  intrinsic  excitability  over 
different  time  scales. 

Putting  it  all  together 

Our  findings  represent  one  of  the  first  comprehensive  demon¬ 
strations  of  rapid  homeostatic  plasticity  of  intrinsic  excitability 
that  results  in  a  stabilization  of  output  in  a  mature,  intact  network 
of  a  CPG.  These  rapid  compensatory  increases  are  mechanisti¬ 
cally  independent,  suggesting  robustness  in  the  maintenance  of 
neural  network  output  that  is  critical  for  survival.  Furthermore, 
this  study  reveals  a  distinct  mechanism  for  compensation  that 
leads  us  to  at  least  one  working  model  of  homeostatic  plasticity  in 
this  system  (Fig.  6).  We  hypothesize  that  one  pathway  to  func¬ 


tional  compensation  in  this  system  relies  on  intracellular  calcium 
concentration  as  a  measure  of  cell  excitability  (Kennedy,  1989; 
Ross,  1989;  LeMasson  etal.,  1993).  A  block  of  A-type  K+  chan¬ 
nels  leads  to  an  increase  in  the  excitability  of  the  cell,  causing 
release  of  intracellular  calcium  stores.  This  calcium  influx  alters 
calcineurin  activity,  resulting  in  dephosphorylation  of  targets, 
perhaps  the  KCa-channels  themselves,  that  increases  IKCa.  This 
ultimately  restores  the  outward  current  balance  that  was  lost  dur¬ 
ing  decreased  IA  and  thus  restores  excitability.  The  corresponding 
mechanism  responsible  for  compensatory  increases  in  IA  (Fig.  6) 
is  presently  unknown,  but  is  likely  to  involve  distinct  mechanisms 
for  monitoring  excitability  (Dirnagl  et  al.,  2003),  as  well  as  other 
candidate  mechanisms  that  alter  A-type  conductances  (Connors 
et  al.,  2008).  Finally,  we  hypothesize  that  positive  coregulation  of 
mRNA  numbers  for  channel  genes  may  ultimately  provide  a  res¬ 
ervoir  of  channel  protein  (Fig.  6)  for  implementing  compensa¬ 
tion  over  rapid  time  scales  during  which  de  novo  synthesis  maybe 
insufficient  to  implement  full  compensation. 
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Nahar  J,  Lett  KM,  Schulz  DJ.  Restoration  of  descending  inputs 
fails  to  rescue  activity  following  deafferentation  of  a  motor  network. 
j  Neurophysiol  108:  871-881,  2012.  First  published  May  2,  2012; 
doi  :10.1152/j n.00183.2012. —  M  otor  networks  such  as  the  pyloric  net¬ 
work  of  the  stomatogastric  ganglion  often  require  descending  neuro- 
modulatory  inputs  to  initiate,  regulate,  and  modulate  their  activity  and 
their  synaptic  connectivity  to  manifest  physiologically  appropriate 
output.  Prolonged  removal  of  these  descending  inputs  often  results  in 
a  compensatory  response  that  alters  the  inputs  themselves,  their 
targets,  or  both.  Using  the  pyloric  network  of  the  crab,  Cancer 
borealis,  we  investigated  whether  isolation  of  motor  networks  would 
result  in  alterations  that  change  the  responses  of  these  networks  to 
restored  modulatory  input.  We  used  a  reversible  block  with  isotonic 
sucrose  to  transiently  alter  descending  inputs  into  the  pyloric  network 
of  the  crab  stomatogastric  ganglion.  Using  this  method,  we  found  that 
blocking  neuromodulatory  inputs  caused  a  reduced  ability  for  subse¬ 
quently  restored  modulatory  projections  to  appropriately  generate 
network  output.  Our  results  suggest  that  this  could  be  due  to  changes 
in  activity  of  descending  projection  neurons  as  well  as  changes  in 
sensitivity  to  neuromodulators  of  the  target  neurons  that  develop  over 
the  time  course  of  the  blockade.  These  findings  suggest  that  although 
homeostatic  plasticity  may  play  a  critical  role  in  recovery  of  func¬ 
tional  output  in  a  deafferented  motor  network,  the  results  of  these 
compensatory  changes  may  alter  the  network  such  that  restored  inputs 
no  longer  function  appropriately. 

neuromodulation;  homeostatic  plasticity;  stomatogastric  ganglion;  de¬ 
centralization 


motor  networks  often  consist  of  central  pattern  generator 
circuits  (CPGs)  that  endogenously  encode  network  output 
features  as  a  result  of  the  intrinsic  properties  of  the  constituent 
neurons  as  well  as  their  pattern  of  synaptic  connectivity  (Grill- 
ner  2006;  Selverston  2010).  M  oreover,  the  ability  to  generate 
this  endogenous  activity  often  requires  both  excitatory  drive  as 
well  as  neuromodulatory  inputs  that  tune  synaptic  and  intrinsic 
properties  of  the  constituent  neurons  to  enable  pattern  gener¬ 
ation  (Harris-Warrick  2011),  such  that  loss  of  these  inputs 
often  results  in  a  loss  of  function,  even  though  the  motor 
networks  themselves  remain  intact. 

More  recently,  a  substantial  amount  of  plasticity  in  this 
dependence  on  descending  drive  and  motor  output  has  been 
observed.  For  example,  in  the  crustacean  stomatogastric  gan¬ 
glion  (STG),  the  pyloric  motor  network  is  a  CPG  that  is 
dependent  on  descending  neuromodulatory  inputs  to  produce  a 
rhythmic  motor  output  (Moulins  and  Cournil  1982).  Yet  al¬ 
though  the  pyloric  motor  output  ceases  soon  after  the  modu- 
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latory  inputs  are  removed,  over  a  period  of  3-7  days  (depend¬ 
ing  on  species  studied)  the  pyloric  rhythm  can  recover  its 
functional  output  in  the  absence  of  all  descending  neuromodu¬ 
lation  (Golowasch  et  al.  1999;  Thoby-Brisson  and  Simmers 
1998,  2002).  These  results  suggest  that  in  the  absence  of  the 
neuromodulatory  inputs  that  produce  ongoing  activity,  the 
intrinsic  motor  circuits  "retune"  themselves  so  that  they  are 
able  to  produce  rhythmic  activity  by  input-independent  mech¬ 
anisms.  These  mechanisms  in  the  STG  include  at  least  in  part 
altering  membrane  conductances  and  ion  channel  expression 
patterns  of  individual  neurons  (Khorkova  and  Golowasch 
2007;  M  izrahi  et  al.  2001;  Thoby-Brisson  and  Simmers  2002). 
This  kind  of  functional  recovery  following  removal  of  com¬ 
mand  inputs  is  not  unique  to  the  STG .  Similar  phenomena  have 
been  reported  in  recovery  of  stepping  function  in  spinal  cord 
transected  cats  (Barbeau  and  Rossignol  1987;  Rossignol  etal. 
2004a, b),  as  well  as  recovery  of  function  in  the  vestibular 
system  following  deafferentation  (Darlington  and  Smith  2000). 

Although  some  properties  of  these  networks  have  been 
studied  subsequent  to  recovery,  little  is  known  about  the 
functional  state  of  these  networks  during  the  process  of  recov¬ 
ery.  in  the  pyloric  network  of  the  STG,  loss  of  neuromodula¬ 
tion  is  sufficient  to  change  relationships  in  individual  neurons 
at  both  the  level  of  ionic  conductances  and  expression  of  ion 
channel  mRNA  prior  to  recovery  taking  place,  detectable  even 
after  24  h  (Khorkova  and  Golowasch  2007;  Temporal  et  al. 
2012).  Y  et  it  is  unclear  whether  these  changes  prior  to  recovery 
have  an  immediate  impact  on  the  ability  of  the  network  to 
generate  neuromodulation-induced  activity.  Our  focus  in  this 
study  is  to  better  understand  whether  the  changes  that  are 
initiated  as  a  result  of  removal  of  modulatory  inputs  influence 
the  ability  of  the  system  to  respond  appropriately  if  descending 
inputs  are  ultimately  restored.  This  has  particular  implication 
for  nerve  injury  and  disease,  because  if  potentially  compensa¬ 
tory  changes  are  initiated  in  networks  deprived  of  afferent 
input,  then  the  effects  of  restoring  input  may  be  unpredictable. 
In  this  study,  we  investigate  the  impact  of  removal  and  resto¬ 
ration  of  descending  inputs  on  functional  output  of  the  pyloric 
network  of  the  STG  over  a  3-  to  4-day  time  frame  of  deaffer¬ 
entation.  Recovery  of  rhythmic  activity  in  the  STG  can  occur 
over  a  time  course  between  1  and  7  days  (Luther  et  al.  2003; 
M  izrahi  et  al.  2001),  although  the  most  common  that  we  have 
observed  is  on  the  order  of  3-5  days,  similar  to  other  reports  in 
the  literature  (Thoby-Brisson  and  Simmers  1998,  2002).  The 
time  frame  of  our  experiments  therefore  likely  coincides  with 
the  initiation  of  cellular  mechanisms  involved  in  the  recovery 
process.  We  then  determine  whether  changes  in  the  output  of 
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these  restored  networks  result  from  changes  in  presynaptic 
inputs,  target  neurons  of  the  pyloric  network,  or  both. 

MATERIALS  AND  METHODS 

Animals  and  the  stomatogastric  preparation.  Adult  crabs,  Cancer 
borealis,  were  obtained  from  The  Fresh  Lobster  Company  (Glouces¬ 
ter,  MA)  and  maintained  in  artificial  seawater  at  12°C  until  used. 
Crabs  were  anesthetized  by  keeping  them  on  ice  for  30  min  before 
dissection.  The  complete  stomatogastric  nervous  system  (STNS)  was 
dissected  out  of  the  animal  and  pinned  out  in  a  Sylgard-coated  (Dow 
Corning)  dish  containing  chilled  ( 12-  13°C )  physiological  saline. 

The  stomatogastric  nervous  system  consists  of  the  STG  motor 
networks  and  their  associated  inputs  (Fig.  1).  The  output  of  the  pyloric 
network  of  the  STG  can  be  most  directly  monitored  via  projections  of 
motor  neurons  through  two  nerve  types,  the  paired  ivns  and  mvns.  The 
STG  receives  central  inputs  by  only  a  single  nerve,  the  stomatogastric 
nerve  ( stn ),  connecting  the  STG  with  the  rest  of  the  crustacean  central 
nervous  system  (Fig.  1),  including  in  this  preparation  the  paired 
commissural  ganglia  (CoGs)  and  the  oesophageal  ganglion  (OG).  The 
stn  of  the  crab,  Cancer  borealis,  contains  approximately  60  large 
axons,  as  well  as  a  bundle  of  smaller  fibers  (Coleman  et  al.  1992).  Of 
the  60  large  fibers,  about  40  descend  from  the  paired  CoGs,  10  from 
the  OG,  and  10  represent  ascending  interneuron  axons  from  the  STG 
(Coleman  et  al.  1992;  Goldberg  et  al.  1988).  These  stn  inputs  are 
responsible  for  initiating,  maintaining,  and  altering  the  activity  of  the 
pyloric  network  (Stein  2009). 

Solutions.  The  physiological  saline  solution  consisted  of  the  fol¬ 
lowing  (in  mM):  440  NaCI,  11  KCI,  13  CaCI2,  26  MgCI2,  and  10 
FIEPES  buffer,  pH  7.45.  In  most  experiments,  the  stomatogastric 
nerve  (stn;  see  Fig.  1)  was  blocked  with  a  solution  composed  of 
isotonic  sucrose  (750  mM ).  In  some  experiments,  neuromodulators 
were  superfused  to  the  STG  to  determine  their  effects  on  preparations 
that  previously  had  stn  activity  blocked  with  sucrose.  These  modula¬ 
tors  include  octopamine  (Sigma  Chemical,  St.  Louis,  MO),  pilo¬ 


carpine  (Acros  Organics,  Geel,  Belgium),  and  proctolin  (American 
Peptide  Company,  Sunnyvale,  CA).  Unless  otherwise  specified, 
chemicals  were  obtained  from  Fisher  Chemical.  When  the  STG  was 
superfused  with  modulators,  the  stn  was  additionally  blocked  with 
10~6  M  tetrodotoxin  (Alomone  Labs,  J erusalem,  Israel)  in  keeping 
with  common  practice  in  STG  preparations,  as  insurance  that  no 
propagation  of  signals  up  and/or  down  the  stn  could  occur. 

Electrophysiology.  For  electrophysiological  recordings,  petroleum 
jelly  wells  were  placed  on  the  nerves  from  which  recordings  would  be 
made  (Ivn,  lateral  ventricular  nerve;  mvn,  medial  ventricular  nerve; 
dgn,  dorsal  gastric  nerve;  stn,  stomatogastric  nerve;  son,  superior 
oesophageal  nerve;  ion,  inferior  oesophageal  nerve).  Extracellular 
recordings  from  the  nerves  were  made  by  placing  stainless  steel  pin 
electrodes  in  the  wells.  Signals  were  amplified  and  filtered  using  a 
differential  AC  amplifier  (A-M  Systems,  Sequim,  WA).  Throughout 
the  experiments,  both  during  recordings  and  incubation  times  be¬ 
tween,  the  preparations  were  maintained  in  chilled  (12-13°C)  physi¬ 
ological  saline.  During  incubation  periods,  penicillin  and  streptomy¬ 
cin  (Sigma)  were  added  to  the  saline  to  prevent  bacterial  infection  of 
the  culture  (50  U/mL  penicillin  and  50  /u.g/mL  streptomycin).  Data 
were  acquired  using  a  Digidata  1322  data  acquisition  board  (Axon 
Instruments,  Sunnyvale,  CA).  In  some  preparations,  the  descending 
projection  neurons  were  stimulated  extracellularly  via  the  stomato¬ 
gastric  nerve  (10  H z,  90  s).  The  stn  was  stimulated  using  a  Grass  S88 
stimulator  (Grass  Instruments,  West  Warwick,  Rl).  Extracellular 
nerve  stimulation  was  achieved  by  placing  the  wires  used  to  record 
nerve  activity  into  a  stimulus  isolation  unit  (Grass  Instruments)  that 
was  connected  to  an  S88  stimulator. 

The  pyloric  rhythm  was  monitored  by  recording  the  activities  of  the 
lateral  pyloric  (LP),  pyloric  constrictor  (PY),  pyloric  dilator  (PD), 
ventricular  dilator  (VD),  and  inferior  cardiac  (1C)  neurons  on  the  Ivn 
and  mvn  (Fig.  1).  Pyloric  burst  frequencies  were  derived  from  the 
cycle  period  as  determined  by  the  time  between  two  consecutive 
bursts  of  PD  neurons.  The  cycle  period  of  the  VD  neuron  was  also 


Fig.  1.  Schematic  of  the  stomatogastric  ner¬ 
vous  system  showing  locations  of  extracel¬ 
lular  recording  sites  as  well  as  representative 
recordings  obtained.  Each  dark  circle  repre¬ 
sents  a  V aseline  well  on  a  given  nerve  from 
which  extracellular  recordings  of  spiking  ac¬ 
tivity  were  made.  Corresponding  labeled 
representative  traces  (made  simultaneously) 
are  shown  from  a  control  preparation  with 
intact  modulatory  projections.  CoG,  com¬ 
missural  ganglion;  OG,  oesophageal  gan¬ 
glion;  ion,  inferior  oesophageal  nerve;  son, 
superior  oesophageal  nerve;  stn,  stomatogas¬ 
tric  nerve;  STG,  stomatogastric  ganglion; 
mvn,  medial  ventricular  nerve;  Ivn,  lateral 
ventricular  nerve;  dgn,  dorsal  gastric  nerve. 
In  each  experiment,  at  least  two  wells  were 
placed  on  the  stn  to  facilitate  recordings 
and/or  blockade  of  descending  projections. 
The  individual  neurons  that  make  up  the 
pyloric  rhythm  as  detected  on  the  Ivn  and  the 
mvn  (LP,  PY,  PD,  VD,  1C)  are  labeled  in  a 
representative  pyloric  burst.  Inset:  location 
of  the  soma  of  the  anterior  gastric  receptor 
(AGR)  neuron  and  its  axonal  projections 
both  into  the  dgn  and  ascending  the  stn 
toward  the  CoGs.  With  this  projection  pat¬ 
tern,  AGR  can  be  detected  in  recordings  of 
the  dgn,  stn,  and  son.  A  series  of  three  AGR 
spikes  detected  on  these  nerves  are  shown  on 
the  representative  recording  with  *. 
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determined  as  the  time  between  VD  bursts  measured  on  the  mvn.  We 
also  used  the  spiking  activity  of  the  anterior  gastric  receptor  (AGR) 
neuron  to  monitor  the  propagation  ability  of  the  stn  across  a  time 
course  of  pharmacologic  blockade.  AGR  is  a  sensory  neuron  that  has 
a  cell  body  in  close  proximity  to  theSTG  and  projects  axon  collaterals 
both  centrally,  through  the  stn  to  both  CoGs  (see  Fig.  1,  inset),  and 
peripherally,  through  the  dgn  toward  the  stomach  musculature.  AGR 
action  potentials  can  be  monitored  often  as  the  largest  spike  on  the  stn, 
and  can  also  be  measured  from  the  paired  sons  and  the  dgn  (Smaran- 
dache  and  Stein  2007).  Multisweep  recordings  (Spike2  v6.0,  Cam¬ 
bridge  Electronic  Design,  Cambridge,  UK)  were  used  to  detect  AGR 
spikes  at  all  of  these  recording  sites.  If  AGR  was  detectable  on  the 
dgn,  stn  (both  in  the  blocking  well  and  the  other;  see  Fig.  1),  and  the 
son  after  removal  of  the  sucrose  blockade,  the  block  was  considered 
as  likely  not  having  damaged  the  stn  signaling  capability.  Any 
preparations  for  which  we  could  not  use  AGR  to  confirm  the  propa¬ 
gation  ability  of  the  stn  were  omitted.  Finally,  we  measured  the 
overall  spiking  frequency  of  the  units  contained  within  the  sons.  For 
son  recordings,  we  measured  the  instantaneous  firing  frequency. 
Flowever,  in  some  of  our  recordings,  bursting  activity  of  the  esoph¬ 
ageal  OD1  neuron  was  detected  on  the  son  (Nagy  et  al.  1981).  W  hen 
bursting  was  present,  we  used  average  firing  frequency  of  the  son 
spikes  between  bursts  of  OD1  for  a  minimum  of  10  interburst 
intervals  as  the  data  for  son  spiking  frequency. 

Data  analysis.  Pyloric  activity  was  measured  as  the  average  burst 
frequency  of  consecutive  cycles  for  10  min  in  which  the  range  of 
values  did  not  change  visibly  and  were  assumed  to  represent  the 
steady  state.  Spiking  activity  of  the  AGR  and  son  were  measured  as 
instantaneous  frequency  or  averagefiring  frequency  of  a  section  of  the 
recording  (durations  ranged  from  30  to  240  s).  For  experiments  where 
the  block  was  removed  and  pyloric  frequency  measured  (seen  in  Figs. 
3  and  4),  relative  frequency  was  plotted  as  the  ratio  (for  a  given 
preparation)  of  the  frequency  at  timex/control  frequency.  We  chose  to 
express  the  data  this  way,  when  possible,  to  ensure  that  differences 
across  groups  after  restored  modulation  were  not  simply  due  to 
differences  in  starting  frequency  (Spitzer  et  al.  2008).  Differences  in 
mean  firing  rates  or  bursting  frequencies  were  compared  using  un¬ 


paired  t-tests.  Correlations  were  analyzed  using  Pearson's  Product 
Moment  tests.  All  statistical  tests  were  performed  using  SigmaPlot 
vll.O  (Systat  Software,  Inc.,  San  Jose,  CA). 


RESULTS 

W  e  used  a  simple  blocking  regime  with  isotonic  sucrose  to 
examine  reversible  neuromodulatory  blockade  effects  on  the 
output  of  the  pyloric  rhythm.  Our  data  show  that  sucrose 
blockade  virtually  eliminates  all  pyloric  activity  on  both  the/vn 
(LP,  PY,  and  PD  activity)  and  the  mvn  (VD  and  1C  activity) 
usually  within  1  h  of  application  (Fig.  2).  Simultaneous  re¬ 
corded  activity  made  from  the  blocking  well  ( stn-a )  and  pos¬ 
terior  to  the  blocking  well  ( stn-b )  show  a  dissociation  of 
spiking  activity  with  the  sucrose  block  in  place  (Fig.  2,  bottom 
traces),  stn-a  recordings  with  the  sucrose  block  in  place  show 
activity  for  units  that  largely  originate  above  the  block  and 
descend  through  the  stn,  whereas  stn-b  recordings  show  units 
that  largely  originate  below  the  block  that  ascend  the  stn. 
M  ultisweep  analyses  of  A  G  R  before  and  after  the  placement  of 
the  block  show  the  inability  of  AGR  spiking  to  propagate  past 
the  sucrose  block  and  into  the  son;  AGR  spiking  is  absent 
on  the  son  when  the  sucrose  block  is  in  place  (Fig.  3C).  Taken 
together,  these  data  indicate  that  isotonic  sucrose  is  sufficient 
to  block  propagation  of  action  potentials  through  the  stn,  likely 
resulting  in  a  loss  of  modulatory  drive  to  the  STG  during  the 
course  of  the  blockade. 

Our  hypothesis  was  that  removal  of  descending  neuromodu¬ 
latory  inputs  to  theSTG,  and  the  subsequent  change  in  activity, 
would  alter  the  properties  of  the  STG  neurons,  or  the  descend¬ 
ing  inputs,  or  both  and  result  in  impaired  ability  of  the  de¬ 
scending  projections  to  appropriately  drive  rhythmic  activity 
following  subsequent  reconnection  of  the  modulatory  drive. 
We  tested  this  hypothesis  using  the  reversible  sucrose  block- 
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Fig.  2.  Effects  of  blocking  the  stn  with  iso¬ 
tonic  sucrose  solution  on  activity  in  the  sto- 
matogastric  nervous  system.  A  cute  effects  of 
sucrose  blockade  on  burst  frequency  of  the 
pyloric  rhythm  are  shown.  Left  panel:  the 
pyloric  burst  frequency  (mean  ±  SD)  as 
measured  by  the  activity  of  the  Ivn.  Right 
panel:  the  activity  of  the  mvn  as  represented 
by  VD  burst  frequency  (mean  ±  SD).  *** 
represents  a  significant  difference  relative  to 
control  (P  <  0.01;  t-test).  Representative 
recordings  are  provided  for  stomatogastric 
nervous  system  activity  under  conditions  be¬ 
fore  (control)  and  1  h  after  sucrose  block 
(block  intact). 
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Fig.  3.  Effects  of  stn  blockade  over  time  on  the  ability  of  restored  connectivity  (when  the  block  is  removed)  to  generate  pyloric  activity.  A:  representative 
recordings  of  stn  (ascending  and  descending  projections)  and  Ivn  (pyloric  activity)  for  two  preparations  over  the  course  of  4  days  in  culture.  All  recordings  shown 
have  the  stn  unblocked.  CONTROL  preparations  were  maintained  in  physiological  saline  for  the  entirety  of  the  experiment.  SUCROSE  preparations  had  the  stn 
blocked  with  isotonic  sucrose  solution.  The  block  was  removed  once  every  24  h,  and  a  series  of  recordings  (as  shown  here)  made  to  determine  whether  rhythmic 
pyloric  activity  was  restored.  Following  these  recordings,  the  block  was  put  back  into  place  for  the  next  24  h.  B:  relative  frequency  of  Ivn  bursting  and  AGR 
spiking  over  time  with  and  without  sucrose  blockade.  Data  for  the  sucrose  group  were  collected  while  the  block  was  removed.  Relative  frequency  for  each 
preparation  was  calculated  as  a  proportion  of  the  intact  (control,  day  0)  activity  of  that  preparation,  and  expressed  here  as  mean  ±  SD .  *  represents  a  significant 
difference  relative  to  control  for  a  given  day  (P  <  0.05;  t-test).  C :  multisweep  recordings  of  AGR  activity  showing  the  reversibility  of  the  sucrose  block  over 
time  in  culture.  AGR  propagation  into  the  son  is  blocked  in  sucrose,  but  restored  when  sucrose  is  removed  from  the  blocking  well,  here  shown  at  24  and  48  h 
following  the  initial  sucrose  block,  stn- a  was  the  blocking  well  in  this  experiment.  During  sucrose  blockade,  AGR  action  potentials  (initiated  in  the  stn  and 
propagated  both  anteriorly  and  posteriorly)  sometimes  can  propagate  into,  but  not  beyond,  the  blocking  well  (see  24-h  blocked  recording). 
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ade.  Descending  inputs  were  blocked  over  the  course  of  4  days. 
Once  per  24-h  period  we  removed  the  stn  block  for  1  h  and 
recorded  the  activity  of  the  pyloric  network  in  response  to  this 
reestablished  modulatory  connectivity.  Representative  record¬ 
ings  of  stn  (descending  and  ascending  inputs)  and  Ivn  (pyloric 
output)  for  a  single  preparation  in  each  treatment  group  (con¬ 
trol  and  sucrose)  across  the  time  course  of  an  experiment  are 
provided  in  Fig.  34.  All  of  the  recordings  in  Fig.  34  and 
measurements  summarized  in  Fig.  38  were  made  during  the 
time  when  the  sucrose  blockade  had  been  removed  for  that  day; 
thus,  these  preparations  aii  feature  potentially  intact  connec¬ 
tions  from  the  descending  projections  to  the  STG.  None  of  the 
preparations  used  in  these  experiments,  or  any  subsequent 
experiments,  had  yet  undergone  full  recovery  of  the  pyloric 
rhythm  in  the  absence  of  descending  projections.  Thus,  in  most 
cases,  there  was  no  pyloric  activity  in  the  blocked  preparations 
prior  to  removal  of  the  block. 

Other  than  an  initial  slight  decrease  in  pyloric  burst  fre¬ 
quency,  control  preparations  showed  consistent  and  robust 
pyloric  output  across  the  time  course  of  the  experiment  (Fig. 
38,  left).  Conversely,  as  the  preparations  spend  more  time  in  a 
decentralized  state  via  sucrose  block,  they  are  subsequently 
less  likely  to  reestablish  pyloric  output  when  the  block  is 
removed  (Fig.  38,  left).  After  24  h,  sucrose-blocked  prepara¬ 
tions  were  still  capable  of  generating  a  pyloric  rhythm  that  was 
not  significantly  different  from  the  control  when  the  blockade 
was  removed,  but  after  48  h  of  blockade  there  was  a  significant 
decrease  in  pyloric  frequency.  After  4  days,  virtually  all  pyloric 
activity  was  eliminated  even  with  the  sucrose  block  removed 
(Fig.  38).  This  effect  was  notthe  resultof  damage  to  the  stn  by 
the  sucrose  blockade.  Using  the  AGR  neuron  as  one  indicator 
of  the  propagation  ability  of  the  stn,  there  was  no  significant 
difference  in  AGR  spiking  frequency  on  the  stn  across  the  time 
course  of  the  experiment  between  control  and  sucrose  blocked 
preparations  (see  Fig.  3,  4  and  8,  right).  Furthermore,  using 
AGR  as  an  indicator  of  stn  propagation  ability,  we  determined 
that  sucrose  likely  was  not  resulting  in  the  death  of  axons  in  the 
stn  over  the  time  course  of  the  experiment.  Figure  3C  shows  a 
representative  multisweep  analysis  demonstrating  the  efficacy 
of  the  blockade  at  24  h,  as  well  as  the  reestablishment  of  AGR 
action  potential  propagation  into  the  son  following  removal  of 
the  blockade  at  24  and  48  h.  By  48  h,  we  already  saw  a 
substantial  decrease  in  pyloric  frequency  with  the  blockade 
removed  (Fig.  38),  yet  stn  propagation  ability  appears  to  have 
remained  intact. 

Our  results  demonstrate  that  there  is  a  decrease  in  the  ability 
of  restoring  connectivity  in  the  stn  and  the  ability  of  the  STG 
to  generate  pyloric  activity  following  a  period  of  removal  of 
the  inputs.  One  possible  mechanism  for  the  decrease  in  func¬ 
tional  output  in  reconnected  stomato gastric  networks  is  a 
change  in  the  activity  of  the  descending  projections  influencing 
pyloric  activity.  Therefore,  we  measured  changes  in  activity  of 
units  found  on  the  other  nerves  that  project  from  the  commis¬ 
sural  ganglia,  son  and  ion.  The  son  and  ion  contain  the 
projections  of  many  of  the  modulatory  neurons  from  the  CoGs 
that  project  down  the  stn  (Coleman  et  al.  1992).  The  ion  has 
only  three  known  units  running  through  it  (H edrich  et  al. 
2011).  Ion  activity  in  our  preparations  largely  consisted  of 
spiking  of  an  esophageal  motor  neuron  (OM  N 1),  with  variable 
activity  of  two  projection  neurons,  MCN1  and  MCN5.  Our 
analysis  of  the  ion  revealed  no  consistent  relationships  in  the 


firing  of  these  three  units  with  blockade  (data  not  shown).  We 
did  not  obtain  all  of  the  appropriate  recordings  to  unambigu¬ 
ously  identify  M  CN 1  and  M  CN5  spikes  on  the  ion.  Flowever, 
our  analysis  showed  there  was  substantial  variability  in  the 
presence  or  absence  of  M  CN  1  and  M  CN5  activity  from  prep¬ 
aration  to  preparation  in  both  control  and  sucrose-blocked 
conditions.  In  other  words,  we  sometimes  saw  only  OMN 
activity,  whereas  in  some  preparations  we  saw  either  one  or 
two  smaller  units  on  the  ion  that  could  correspond  to  M  CN  1 
and/or  M  CN5.  When  these  units  were  present  in  our  record¬ 
ings,  they  also  varied  between  simply  spiking  or  bursting. 
Overall,  although  wecould  not  definitively  quantify  theM  CN1 
and  MCN5  activity,  we  detected  no  consistent  qualitative 
relationships  between  ion  unit  activity  (present  or  absence)  or 
spiking  pattern  and  pyloric  frequency  in  this  study  (data  not 
shown),  although  previously  such  a  relationship  was  reported 
(FI edrich  et  al.  2011).  Flowever,  we  did  see  a  relationship  with 
son  spiking  activity  across  the  time  course  of  our  blockade 
experiments.  The  son  contains  many  more  fibers  than  the  ion, 
including  both  ascending  and  descending  units,  making  iden¬ 
tification  of  individual  units  far  more  challenging.  To  examine 
whether  more  global  patterns  of  activity  in  the  son  may  relate 
to  the  activity  of  the  pyloric  rhythm,  we  measured  the  overall 
spike  frequency  in  the  son  as  it  related  to  the  STG  activity.  The 
spike  frequency  of  the  son  significantly  decreased  over  the  time 
course  of  the  experiment  with  sucrose  blockade  (Fig.  44),  in  a 
pattern  reminiscent  of  that  seen  for  the  pyloric  burst  frequency 
in  our  earlier  experiments  (see  Fig.  38).  We  therefore  deter¬ 
mined  whether  son  spiking  and  pyloric  bursting  activity  were 
correlated.  Our  results  show  that  overall  pyloric  burst  fre¬ 
quency  and  son  spiking  frequency  were  significantly  correlated 
across  blocking  conditions  in  this  study  (Fig.  48).  Flowever, 
these  analyses  do  not  fully  explain  the  relationships  seen  in  our 
data.  When  we  plotted  the  relationships  among  son  spike 
frequency,  AGR  spike  frequency,  and  the  pyloric  burst  fre¬ 
quency  in  individual  preparations,  we  saw  a  surprising  shift  in 
these  relationships.  In  control  preparations  maintained  in  phys¬ 
iological  saline  for  3  days,  we  see  a  significant  correlation 
between  son  spike  frequency  and  pyloric  burst  frequency  (Fig. 
4C,  top),  but  no  relationships  between  AGR  and  the  pyloric 
output,  orson  and  AGR  activity  (Fig.  4C,  middle  and  bottom). 
Conversely,  preparations  in  which  the  stn  has  been  blocked 
with  sucrose  over  the  course  of  3  days  have  no  relationship 
between  son  and  the  pyloric  output  (Fig.  4C,  top  and  middle), 
but  a  significant  correlation  between  AGR  spiking  and  the 
pyloric  burst  frequency  (Fig.  4C,  bottom).  These  results  show 
that  there  is  not  simply  a  scaling  of  all  activity  in  the  prepa¬ 
ration  over  time,  but  rather  that  there  are  distinct  alterations  to 
spiking  activities  of  projection  neurons  that  correspond  to 
changes  in  pyloric  network  output  in  the  STG. 

The  overall  incremental  effects  of  sucrose  blockade  and 
removal  are  summarized  in  Fig.  5.  All  of  the  recordings  in  Fig. 
5  have  been  made  while  the  sucrose  block  has  been  removed 
and  replaced  with  physiological  saline,  effectively  restoring 
communication  between  the  descending  and  ascending  neurons 
in  the  stn.  The  control  recording  (Fig.  5,  control)  shows  the 
activity  of  an  intact  preparation  freshly  dissected  from  the 
animal.  When  projections  are  restored  following  24  h  of 
sucrose  blockade  of  the  stn,  the  pyloric  network  generates 
rhythmic  activity  that  is  reduced  in  bursting  frequency  from 
controls  (Fig.  5,  24-h  block).  At  this  time  point  we  see  con- 
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served  AGR  spike  propagation  via  multisweep  analyses  across 
the  son,  dgn,  and  both  stn  recording  sites  (Fig.  5,  24-h  block, 
right),  ensuring  the  stn  remains  capable  of  spike  propagation. 
At  24  h  we  can  also  see  a  reduction  in  activity  on  the  son  as 
well.  After  48  h  deprived  of  inputs,  removal  of  the  blockade 
results  in  further  weakened  activity  on  the/vn  and  mvn,  as  well 
as  reduced  activity  of  the  son.  After  72  h  deprived  of  inputs, 
removal  of  the  blockade  is  insufficient  to  restore  activity  to  the 
pyloric  network,  although  the  individual  pyloric  units  can  still 
be  seen  to  be  present  on  the  Ivn  recording.  Again,  spikes  are 
seen  propagating  through  the  blocking  site  (Fig.  5,  72-h  block, 
right )  confirming  the  viability  of  the  projections,  and  son 
activity  is  further  decreased  (Fig.  5,  72-h  block).  However,  in 
this  preparation  we  also  saw  a  substantive  change  in  the 
activity  of  the  AGR  neuron.  It  begins  to  fire  two  independent 
spikes  (instead  of  a  single,  tonic  spike),  one  of  which  propa¬ 
gates  through  the  son,  stn,  and  dgn  that  appears  to  originate  in 


the  stn,  whereas  a  second  independent  spike  propagates  only  in 
the  dgn  in  our  recordings  (Fig.  5,  72-h  block,  far  right ;  see  also 
Dauretal.  2009). 

Although  our  measurements  of  son  spiking  suggest  that  the 
activity  of  modulatory  and  interneuron  projections  may  be 
altered  by  the  sucrose  blockade,  another  possible  explanation  is 
that  decreased  pyloric  activity  could  be  the  result  of  a  loss  in 
the  ability  of  modulatory  projections  to  release  neuromodula¬ 
tor.  We  tested  whether  the  descending  modulatory  projections 
are  still  capable  of  releasing  neuromodulator  after  sucrose 
blockade  by  stimulating  the  stn  both  above  and  below  the  site 
of  the  blockade  after  3  days  of  blockade  with  sucrose.  If 
sucrose  blockade  has  eliminated  the  ability  of  the  modulatory 
projections  to  influence  the  pyloric  rhythm,  then  we  predict 
stimulation  will  be  ineffective  in  producing  output  changes  in 
the  STG.  This  was  not  the  case;  after  3  days  of  sucrose 
blockade,  electrical  stimulation  either  below  the  blocking  well 


I  Neurophysiol  •  doi :  10.1152/j n.00183.2012  •www.jn.org 


Downloaded  from  http://jn.physiology.org/  at  University  of  Missouri-Columbia  on  August  8,  2012 


ISOLATED  NETWORKS  FAIL  TO  RESPOND  TO  RESTORED  INPUTS 


877 


sln-e  >'f  hhmH  ♦'  t»HH  ' 


-V/1 - - 


CONTROL 


sln-b  t"H'- 

dgn  '''I'M  j 


H"l-"  I  w.-v 


Pf~ 


— 


T 


~ 


After  24h  Block 


After  48h  Block 


ivn  HM(l  HiHiw  IIH'H»#l  |W»«i|  WM»  W*HHH 

mvn  '#)'■))■  III!  )!■  'HW  ■jt  'W'  |(  »  ffH”  )(■  '  til  I 


son  44li»»*»<*fH4ki»W»^ 

sin- a  t 

sln-b  'f  1  '  j!**  |i||"||  4- ft  j|  i|«  4 

dgn  1  i  1  |  1  H  I  1 


son  . . .  Ill  H-4t-k-m‘AltHm|Xl^i  lil'i  |>  I 

sln-a  >"  I-  »  >  nK  i»,-  >■  t  "■»  «f  . . .  >  . 

stn-b  1 1 1  I"  ♦  t  . . .  »■  H‘-»f  »"»■  »■  f" 


■  stn-b 


Mm 


<  dgn 


-\P*  -—  ! 


dgn 


Ivn  4ff- 


~wt— H 


— — 


— 


-H 


After  72h  Block 


4\1 


•  sln-b 


1  dgn 


Second  AGR  spike  trigger 


First  AGR  spike  tngger 


sln-b 


dgn 


•4Mm 


i  son 
sln-a 
sln-b 
'  dgn 

i  dgn 


Fig.  5.  Representative  recordings  of  a  single  sucrose-blocked  preparation  over  four  time  points  to  illustrate  the  effects  stn  blockade  has  on  activity  of  the 
stomatogastric  nervous  system,  stn- a  is  the  blocking  well  in  this  experiment.  All  recordings  were  made  with  the  sucrose  block  removed  for  30  min.  Simultaneous 
recordings  of  the  son,  stn  (two  sites),  dgn,  Ivn,  and  mvn  provided  over  72  h  of  sucrose  blockade.  In  the  72-h  panel,  representative  time-locked  AGR  spikes  are 
denoted  by  arrowheads.  Right  panel-,  multisweep  recordings  for  the  son,  dgn,  and  the  two  stn  recording  sites  (stn-a  and  stn-b)  triggered  from  AGR  spikes  recorded 
on  the  dgn.  In  the  72-h  trace,  the  multisweep  results  show  differences  from  triggering  on  either  the  second  AGR  spike  in  the  spike  pairs  of  the  dgn  (left)  or  the 
first  AGR  spike  (right),  dgn  multisweeps  triggered  from  the  first  AGR  spike  are  also  provided  at  a  longer  time  scale  to  show  the  independence  of  the  first  and 
second  AGR  spikes  on  the  dgn. 


(with  the  block  still  intact;  Fig.  64)  or  above  the  blockade  (with 
the  block  removed;  Fig.  68)  resulted  in  similar  ability  to  drive 
the  pyloric  frequency.  This  stimulation  effect  is  reversible,  as 
20  min  after  stimulation  ceased,  the  pyloric  rhythm  had  largely 
returned  to  its  previous  baseline  levels  (Fig.  6,  A  and  8).  These 
results  demonstrate  that  the  stn  fibers  are  still  capable  of 
releasing  neuromodulator  to  influence  pyloric  activity.  Further¬ 
more,  these  results  provide  strong  support  that  we  are  not 
damaging  stn  axons  via  sucrose  blockade,  because  we  get  the 
same  stimulation  effect  whether  we  stimulate  above  or  below 
where  the  sucrose  blockade  was  placed. 

A  third  mechanism  by  which  reconnected  projections  might 
fail  to  appropriately  generate  network  output  is  through  altered 
sensitivity  to  neuromodulation  of  the  target  neurons  of  the 
pyloric  network.  We  tested  this  by  analyzing  the  responses  of 
isolated  pyloric  networks  (subsequent  to  sucrose  blockade)  to 
exogenously  applied  neuromodulators.  After  3  days  either  with 
blockade  via  sucrose,  or  as  an  unblocked  control  preparation, 
the  stn  of  each  preparation  was  blocked  with  tetrodotoxin  to 
ensure  no  interaction  of  the  superfused  modulator  with  de¬ 
scending  neuromodulation.  We  then  perfused  a  range  of  in¬ 


creasing  concentrations  of  the  neuropeptide  proctolin,  the  bio¬ 
genic  amineoctopamine,  or  the  muscarinic  agonist  pilocarpine, 
all  of  which  are  known  to  exogenously  activate  pyloric  activity 
(Goaillard  et  al.  2004;  Flooper  and  M  arder  1987;  Swensen  and 
M  arder  2001),  on  to  the  isolated  ST G  preparation  and  recorded 
pyloric  burst  output.  We  found  a  significant  increase  in  re¬ 
sponse  to  proctolin  in  preparations  that  had  previously  been 
blocked  with  sucrosefor  3  days  relative  to  control  (Fig.  7,  top). 
At  three  of  four  concentrations  tested  (10~8  M  ,  10~'  M  ,  and 
10~6  M  )  there  was  a  significantly  higher  frequency  of  pyloric 
burst  output  in  the  sucrose-blocked  preparations  relative  to 
controls.  Conversely,  there  were  no  significant  differences  in 
response  to  octopamine  between  sucrose  and  control  prepara¬ 
tions  (Fig.  7,  middle).  Finally,  we  found  a  significant  decrease 
in  response  to  pilocarpine  in  sucrose-blocked  preparations 
(Fig.  7,  bottom).  Although  burst  frequency  responses  of  the 
pyloric  rhythm  were  altered  in  response  to  different  modula¬ 
tors,  there  appeared  to  be  no  substantial  difference  in  the 
pattern  of  the  pyloric  rhythm  elicited  by  modulators  between 
control  and  sucrose-blocked  preparations.  Representative  re¬ 
cordings  time-scaled  to  compare  one  complete  pyloric  burst 
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Fig.  6.  Effects  of  stn  stimulation  after  3  days  of  blockade  on  pyloric  burst  frequency.  In  these  experiments,  three  wells  were  made  on  the  stir,  the  middle  was 
used  for  blockade,  whereas  those  above  (upper  stim)  and  below  (lower  stim)  were  sites  of  extracellular  stimulation.  Preps  either  were  maintained  with  the  block 
intact  and  the  stn  stimulated  below  the  blocking  well  (filled  bars;  A)  or  unblocked  and  stimulated  above  the  blocking  well  (open  bars;  B).  A;  pyloric  frequency 
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blocking  well.  Recordings  were  made  in  parallel  fashion  to  sucrose  block  experiments.  *  represents  a  significant  difference  relative  to  control  for  a  given  group 
(P  <  0.05;  t-test).  Sample  size  indicated  in  each  bar.  Individual  pyloric  units  (LP,  PY,  and  PD)  are  labeled  as  recorded  on  the  Ivn  for  reference. 


pattern  suggest  that  although  proctolin,  octopamine,  and  pilo¬ 
carpine  elicit  somewhat  distinct  patterns  (e.g.,  LP  activity 
levels),  no  differences  were  apparent  between  control  and 
sucrose  rhythms  for  a  given  modulator  (Fig.  7).  Together,  these 
results  suggest  that  the  sucrose  blockade  does  indeed  result  in 
a  change  in  neuromodulator  sensitivity  of  pyloric  neurons, 
perhaps  in  particular  for  those  neurons  involved  in  the  gener¬ 
ation  and  frequency  of  the  bursting  of  the  pyloric  rhythm, 
specifically  AB  and/or  PD. 

DISCUSSION 

It  is  only  somewhat  recently  that  the  occurrence  of,  and 
mechanisms  underlying,  homeostatic  compensation  in  neural 
circuits  have  been  more  widely  recognized  and  studied 
(Marder  and  Tang  2010;  Turrigiano  2008).  Such  studies  in¬ 
clude  examples  of  homeostatic  plasticity  as  compensation  for 
gene  knockouts  (Bergquist  et  al.  2010),  synaptic  scaling  (Tur¬ 
rigiano  2008),  sensory-modality  specific  compensation  (Deeg 
and  Aizenman  2011),  and  recovery  of  function  in  networks 
deprived  of  sensory  input  (Cullen  et  al.  2009)  and  descending 
command  fibers  (Golowasch  et  al.  1999;  Tillakaratne  et  al. 
2010).  Y  et  despite  the  substantial  amount  that  we  have  learned 
about  homeostatic  plasticity  in  these  systems,  very  few,  if  any, 
studies  have  examined  the  impact  of  these  homeostatic  pro¬ 
cesses  on  the  subsequent  ability  of  the  network  to  appropriately 
respond  to  further  perturbation  or  even  reversal  of  conditions 


that  initially  triggered  compensation.  Our  study  is  one  of  the 
first  to  examine  effects  of  removal  and  restoration  of  inputs  to 
a  motor  network  over  the  time  course  of  homeostatic  compen¬ 
sation.  Our  data  demonstrate  that  following  blockade  of  mod¬ 
ulatory  inputs,  the  system  is  altered  such  that  removal  of  the 
block  is  insufficient  to  restore  network  activity.  These  effects 
do  not  seem  to  simply  be  due  to  damage  of  the  descending  stn 
fibers.  Using  AG  Ft  activity  as  a  monitor  of  stn  propagation 
ability  (Daur  et  al.  2009),  we  found  no  significant  differences 
in  AGFt  spike  propagation  ability  over  time  in  both  control  and 
sucrose  preparations,  indicating  the  propagation  ability  of  the 
stn  is  likely  maintained  throughout  these  experiments. 

W  e  investigated  whether  these  effects  are  due  to  changes  in 
the  descending  projections  themselves,  responses  of  the  pyloric 
network,  or  both.  T o  determine  whether  blockade  causes  a  loss 
of  modulator  release  from  descending  fibers,  we  electrically 
stimulated  the  stn  both  above  and  below  the  blocking  well  after 
sucrose  blockade.  Stimulation  was  equally  effective  in  eliciting 
pyloric  output  regardless  of  whether  the  stimulation  was  below 
the  blockade  (with  the  block  intact)  or  above  the  blockade 
(with  the  block  removed).  Thus,  although  we  are  unable  to 
determine  whether  quantitative  changes  occur  in  modulator 
release  as  a  result  of  our  blocking  manipulation,  at  least 
qualitatively  we  determined  that  the  ability  to  release  neuro¬ 
modulator  is  not  lost  from  presynaptic  terminals  of  descending 
projections,  nor  is  the  ability  of  theSTG  to  respond  to  released 


I  Neurophysiol  •  doi :  10.1152/j n.00183.2012  •www.jn.org 


Downloaded  from  http://jn.physiology.org/  at  University  of  Missouri-Columbia  on  August  8,  2012 


ISOLATED  NETWORKS  FAIL  TO  RESPOND  TO  RESTORED  INPUTS 


879 


Proctolin 


CONTROL 

LVN  . . |  . . f 

SUCROSE 


CONTROL 


CONTROL 


LVN 


SUCROSE 

* — m 


2  sec 


Fig.  7.  Effects  of  exogenously  applied  neuromodulators  on  pyloric  burst  frequency  (mean  ±  SD)  after  3  days  of  stn  blockade.  After  3  days  of  sucrose,  or  time 
in  culture  control,  the  stn  was  completely  blocked  with  tetrodotoxin  (TTX)  to  ensure  elimination  of  any  endogenous  modulator  release.  Increasing  concentrations 
of  proctolin,  octopamine,  and  pilocarpine  were  superfused  on  preparations  and  pyloric  burst  frequency  measured  at  steady  state  for  a  given  concentration.  * 
represents  a  significant  difference  relative  to  control  preparations  for  a  given  concentration  (P  <  0.05;  t-test).  Representative  Ivn  recordings  for  control  and 
sucrose  blocked  preparations  are  shown  during  the  perfusion  of  10~6  M  proctolin  (top),  10~4  M  octopamine  (middle),  and  10~5  M  pilocarpine  (bottom).  A  single 
pyloric  burst  is  scaled  for  each  recording  to  compare  the  pattern  elicited  by  modulator  application  in  each  treatment  group.  Individual  pyloric  units  (LP,  PY ,  and 
PD)  are  labeled  for  reference. 


modulator  lacking.  These  results  also  provide  support  for  the 
maintenance  of  stn  propagation  ability  throughout  the  study, 
because  stimulation  above  a  removed  block  was  able  to  trigger 
effects  that  traversed  the  area  of  blockade. 

In  addition  to  blocking  descending  fibers,  the  blockade  may 
also  influence  any  ascending  feedback  from  the  STG  to  the 
commissural  ganglia,  as  represented  in  our  studies  by  AGR 
activity  (Daur  et  al.  2009).  In  this  study  we  consider  both 
possibilities,  but  do  not  distinguish  between  the  two.  We 
analyzed  the  activity  of  the  son  nerve,  which  predominantly 
contai  ns  descendi  ng  projecti ons  from  the  C  oG  s  (C  ol eman  et  al . 


1992),  and  found  a  significant  decrease  in  firing  rate  on  the  son 
with  increasing  duration  of  blockade  that  mirrored  changes 
seen  in  pyloric  frequency.  Indeed,  there  was  a  strong  quanti¬ 
tative  correlation  overall  between  son  firing  frequency  and 
pyloric  burst  frequency  among  control  and  sucrose  blocked 
preparations.  Furthermore,  there  is  a  substantial  change  in  the 
relationship  of  son  activity  to  the  pyloric  output  over  the  time 
course  of  the  blockade,  in  control  preparations,  son  spiking  and 
the  output  of  the  pyloric  network  are  coupled.  Y  et  as  a  result  of 
blocking  the  stn  projections,  this  relationship  completely  dis¬ 
appears,  and  pyloric  output  becomes  strongly  correlated  with 
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the  spiking  activity  of  AGR.  We  have  even  seen  dramatic 
qualitative  changes  in  the  activity  of  neurons  like  AGR  in  this 
system  (Fig.  5,  72-h  block )  that  may  be  further  indicative  of 
alterations  in  projections  neurons  that  influence  STG  activity. 
We  cannot  determine  from  these  experiments  whether  these 
changes  in  son,  AGR,  and  pyloric  network  activity  are  caus¬ 
ative  or  correlative.  When  combined  with  the  fact  that  stimu¬ 
lation  experiments  reveal  these  projections  have  not  lost  the 
capacity  to  release  modulators,  the  data  overall  suggest  that 
disruption  of  interactions  between  the  STG  and  the  higher 
ganglia  (CoG  and  OG)  may  lead  to  a  change  in  activity 
underlying  modulation,  but  not  the  capacity  for  modulator 
release,  over  time. 

Finally,  we  determined  whether  alteration  of  descending 
projections  affects  response  to  exogenously  applied  neuro¬ 
modulators.  Proctolin  is  a  neuropeptide  that  is  a  major  constit¬ 
uent  of  releasate  of  presynaptic  fibers  of  theSTN  (Billimoria  et 
al.  2005)  and  is  able  to  exogenously  restore  network  activity  to 
deafferented  preparations  (Golowasch  and  M  arder  1992; 
Flooper  and  M  arder  1987).  There  is  a  significant  increase  in  the 
effect  of  proctolin  on  pyloric  frequency  in  preparations  that 
have  been  blocked  with  sucrose  relative  to  control  prepara¬ 
tions.  It  is  not  possible  to  determine  from  our  data  whether  this 
is  a  result  of  increased  receptor  signaling  activity,  or  down¬ 
stream  targets  of  modulation  such  as  voltage-gated  conduc¬ 
tances,  or  both.  Flowever,  there  was  not  an  increase  in  pyloric 
cycle  frequency  sensitivity  for  all  neuromodulators;  octo- 
pamine  application  (Goaillard  et  al.  2004)  did  not  result  in 
significant  differences  in  pyloric  cycle  frequency  between 
control  and  sucrose-blocked  preparations,  whereas  response  to 
pilocarpine  (Swensen  and  M  arder  2001)  was  significantly 
lower  in  sucrose-blocked  preparations  than  in  controls.  These 
data  demon strate  that  there  is  not  a  uniform  change  in  response 
to  all  modulators  in  the  system,  but  rather  distinct  responses  to 
different  classes  of  modulatory  compounds.  This  could  be  the 
result  of  changes  in  modulator  receptors,  or  perhaps  intracel¬ 
lular  pathways  involved  in  mediating  modulator  response. 
However,  these  data  suggest  that  the  convergent  target  of 
several  modulators  in  this  system,  a  mixed-inward  current 
elicited  by  modulator  exposure  (Golowasch  and  M  arder  1992), 
may  not  be  the  source  of  the  differing  responses.  Both  proc¬ 
tolin  and  pilocarpine  influence  this  current  (Swensen  and 
M arder  2000)  in  similar  fashion,  yet  our  preparations  show 
divergent  alterations  in  the  responses  to  these  two  modulators. 

The  changes  we  see  in  the  ability  of  restored  projections  to 
elicit  output  following  different  blockade  regimes  must  repre¬ 
sent  a  change  at  the  cellular  level,  both  in  the  ability  of  the  STG 
neurons  to  respond  to  neuromodulatory  substances,  and  in 
presynaptic  activity  and  subsequent  release  of  modulator.  It  is 
already  established  that  pyloric  neurons  alter  some  intrinsic 
properties  in  response  to  decentralization.  These  include 
changes  in  the  magnitude  of  ionic  currents  such  as  / H  (Thoby- 
Brisson  and  Simmers  2002),  / A  (Khorkova  and  Golowasch 
2007;  Mizrahi  et  al.  2001),  lBa,  and  /KCa  (Khorkova  and 
Golowasch  2007),  and  disruption  of  correlations  among  ionic 
conductances  as  well  as  channel  mRNA  levels  (Khorkova  and 
Golowasch  2007;  Temporal  et  al.  2012).  Our  results  suggest 
that  one  possible  impact  of  these  changes  is  that  subsequently 
restored  neuromodulation  may  be  ineffective  at  generating 
typical  output  because  the  downstream  targets  of  modulators 
have  been  altered  by  decentralization.  Indeed,  due  to  the  fact 


that  these  changes  have  been  documented  previously,  it  is 
possible  that  such  changes  in  intrinsic  excitability  alone  are 
sufficient  to  account  for  our  results. 

Less  is  known  about  changes  in  receptor  expression  follow¬ 
ing  decentralization  in  the  STG.  Our  data  demonstrate  changes 
in  responses  to  proctolin  and  pilocarpine  in  blocked  prepara¬ 
tions,  but  at  this  time  we  cannot  determine  whether  this  is  an 
effect  of  receptor  expression,  activity,  or  signaling,  or  due  to 
changes  in  downstream  targets.  In  other  systems,  changes  in 
receptor  expression  following  deafferentation  have  been  re¬ 
ported,  including  5HT-2C  receptor  expression  in  mammalian 
spinal  cord  (M  urray  etal.  2010).  However,  not  all  systems  that 
undergo  homeostatic  compensation  have  corresponding 
changes  in  receptor  expression.  For  example,  vestibular  deaf¬ 
ferentation  leads  to  plasticity  and  compensation,  but  changes  in 
expression  and  distribution  of  •y-aminobutyric  acid  (GABA) 
(Gliddon  et  al.  2005),  glutamate  (King  et  al.  2002),  glucocor¬ 
ticoid  (Lindsay  et  al.  2005),  and  cannabinoid  (Ashton  et  al. 
2004)  receptors  have  not  been  detected,  despite  changes  in 
responses  to  these  compounds. 

The  interaction  of  neuromodulation  and  neural  network 
output  is  a  complex  interplay  that  occurs  over  multiple  time 
scales.  A  good  deal  is  known  of  how  neuromodulators  affect 
plasticity  in  these  networks  over  shorter  time  scales,  as  well  as 
influence  properties  of  cells  and  synapses  over  more  organiza¬ 
tional  timescales  (Harris-Warrick  2011).  However,  as  progress 
is  made  toward  understanding  regenerative  processes  in  nerve 
injury,  it  becomes  more  imperative  that  potential  impacts  of 
changes  in  isolated  or  deafferented  neural  networks  as  a  result 
of  loss  of  modulatory  inputs  be  more  completely  understood  as 
well. 
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Ransdell  J L ,  Temporal  S,  West  NL,  Leyrer  ML,  Schulz  DJ. 

Characterization  of  inward  currents  and  channels  underlying  burst  activ¬ 
ity  in  motoneurons  of  crab  cardiac  ganglion.]  Neurophysiol  110:  42-54, 
2013.  First  published  April  10,  2013;  doi:10.1152/jn.00009.2013.- 
Large  cell  motoneurons  in  the  Cancer  borealis  cardiac  ganglion 
generate  rhythmic  bursts  of  action  potentials  responsible  for  cardiac 
contractions.  While  it  is  well  known  that  these  burst  potentials  are 
dependent  on  coordinated  interactions  among  depolarizing  and  hyper- 
polarizing  conductances,  the  depolarizing  currents  present  in  these 
cells,  and  their  biophysical  characteristics,  have  not  been  thoroughly 
described.  In  this  study  we  used  a  combined  molecular  biology  and 
electrophysiology  approach  to  look  at  channel  identity,  expression, 
localization,  and  biophysical  properties  for  two  distinct  high-voltage- 
activated  calcium  currents  present  in  these  cells:  a  slow  calcium 
current  (lCaS)  and  a  transient  calcium  current  (lCaT).  Our  data  indicate 
that  CbCaVl  is  a  putative  voltage-gated  calcium  channel  subunit  in 
part  responsible  for  an  L -type  current,  while  CbCaV2  (formerly 
cacophony)  is  a  subunit  in  part  responsible  for  a  P/Q-type  current. 
These  channels  appear  to  be  localized  primarily  to  the  somata  of  the 
motoneurons.  A  third  calcium  channel  gene  (CbCaV3)  was  identified 
that  encodes  a  putative  T-type  calcium  channel  subunit  and  is  ex¬ 
pressed  in  these  cells,  but  electrophysiological  studies  failed  to  detect 
this  current  in  motoneuron  somata.  In  addition,  we  identify  and 
characterize  for  the  first  time  in  these  cells  a  calcium-activated 
nonselective  cationic  current  (Zc AN),  as  well  as  a  largely  noninacti¬ 
vating  TTX -sensitive  current  reminiscent  of  a  persistent  sodium 
current.  The  identification  and  further  characterization  of  these  cur¬ 
rents  allow  both  biological  and  modeling  studies  to  move  forward 
with  more  attention  to  the  complexity  of  interactions  among  these 
distinct  components  underlying  generation  of  bursting  output  in  mo¬ 
toneurons. 

stomatogastric  ganglion;  calcium  current;  calcium  channel;  CAN 
current;  cardiac  ganglion 


bursts  of  action  potentials  are  vital  to  appropriate  motor 
output  in  central  pattern  generator  networks.  The  organization 
of  burst  potentials  often  relies  on  calcium-  and  sodium-medi¬ 
ated  depolarizing  currents  that  have  an  influence  on  neuronal 
function  beyond  eliciting  a  direct  change  in  membrane  poten¬ 
tial.  For  example,  calcium  influx  via  voltage-gated  calcium 
currents  is  known  to  affect  second  messenger  systems  involved 
in  transcriptional  regulation  (Clapham  2007),  activate  other 
ionic  currents  (Sah  1995;  W  isgirda  and  Dryer  1994;  Y  amoah  et 
al.  1994),  as  well  as  trigger  neurotransmitter  release  at  the 
synapse  (M  ulkey  and  Zucker  1993).  Sodium  has  been  shown  to 
activate  sodium-dependent  potassium  currents  and  influence 
neuronal  excitability  (Hage  and  Salkoff  2012;  Rose  2002). 
These  higher-order  cation  effects  are  also  known  to  act  as  the 
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mechanism  behind  different  forms  of  plasticity  in  the  ongoing 
output  of  a  cell  and/or  network  (Cummings  et  al.  1996;  George 
et  al.  2012;  Turrigiano  2008).  However,  before  we  can  fully 
understand  the  role  of  inward  currents  and  resulting  cation  influx 
during  a  burst  potential,  we  must  have  a  thorough  understanding 
of  what  depolarizing  components  exist  in  a  given  cell  type. 

The  large  cell  (LC)  motoneuron  of  the  Cancer  borealis  cardiac 
ganglion  (CG)  has  been  a  successful  model  for  studying  how 
underlying  conductances  interact  to  mediate  burst  potentials  (Ball 
et  al.  2010;  Franklin  et  al.  2010;  Tobin  et  al.  2009),  and  new 
evidence  shows  there  is  rapid  calcium-dependent  homeostatic 
compensation  among  potassium  currents  that  impacts  both  LC 
intrinsic  excitability  as  well  as  appropriate  network  bursting 
(Ransdell  et  al.  2012).  However,  there  is  little  characterization  of 
what  depolarizing  properties  are  active  in  the  C.  borealis  LC 
burst,  hindering  a  more  complete  understanding  of  how  the  LC 
bursts  are  generated  and  maintained.  The  goal  of  this  study  was  to 
characterize  the  depolarizing  components  intrinsic  to  the  LC  mo¬ 
toneuron  in  the  C.  borealis  CG. 

Despite  the  extensive  use  of  crustacean  motor  systems  as 
models  for  understanding  neuronal  network  function  and  mod¬ 
ulation  (Marder  and  Bucher  2007),  the  inward  currents  have 
been  less  thoroughly  studied  than  their  outward  counterparts. 
In  neurons  of  the  stomatogastric  ganglion  (STG),  total  voltage- 
activated  calcium  current  has  been  successfully  measured  by 
using  sensitivity  to  cadmium  chloride  (CdCI2)  as  an  effective 
blocker  (Golowasch  and  M  arder  1992).  Further  pharmacolog¬ 
ical  dissection  of  STG  calcium  currents  suggests  there  are  two 
distinct  high-threshold  calcium  currents  present:  an  L-type 
(CaVl)  and  a  P/Q-type  ( C aV 2 )  current  (Hurley  and  Graubard 
1998;  Johnson  et  al.  2003).  M  olecular  cloning  studies  support 
the  presence  of  these  two  calcium  channel  types  in  neurons  of 
the  STG  (French  etal.  2002).  Because  the  calcium  dependence 
of  synaptic  transmission  in  STG  neurons  has  a  more  hyperpo- 
larized  threshold  (-60  mV)  than  these  two  high-threshold 
currents,  there  is  thought  to  be  an  additional  low-threshold  or 
T-type  ( C aV 3 )  current  involved  in  neurotransmitter  release 
(Graubard  et  al.  1983;  Harris-Warrick  et  al.  1992).  In  the  CG, 
some  analysis  of  calcium  currents  has  taken  place  in  the  LCs 
(Tazaki  and  Cooke  1979,  1986),  particularly  in  the  lobster 
Homarus  americanus  (Tazaki  and  Cooke  1990).  However,  a 
more  thorough  analysis  of  these  currents  is  necessary  to  un¬ 
derstand  their  role  in  generating  burst  output  in  LC  motoneu¬ 
rons. 

Beyond  calcium  currents,  evidence  for  persistent  sodium 
currents  [tetrodotoxin  (TTX)  sensitive]  has  also  been  reported 
in  both  the  STG  (Elson  and  Selverston  1997;  T urrigiano  et  al. 
1995)  and  the  CG  LCs  of  the  lobster  H .  americanus  (Berlind 
1993).  The  somata  of  these  crustacean  motoneurons  do  not 
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contain  spike-mediating  sodium  currents  (Golowasch  and 
M  arder  1992).  However,  a  characterization  of  persistent  so¬ 
dium  currents  has  yet  to  be  performed  in  these  crustacean 
systems.  In  addition,  plateau  potentials  in  the  dorsal  gastric 
neurons  of  the  STG  depend  on  activation  of  a  calcium-depen- 
dent  voltage-independent  cationic  current,  or  "CAN  current" 
(Zhang  and  Harris-Warrick  1995;  Zhang  etal.  1995),  but  this 
current  has  yet  to  be  identified  or  characterized  in  LC  motoneu¬ 
rons. 

With  the  C.  borealis  LC  motoneuron  being  the  focus  of 
several  recent  studies  exploring  how  intrinsic  conductances 
and  channels  are  organized  and  conserved  during  burst  gener¬ 
ation  (Ball  et  al.  2010;  Franklin  et  al.  2010;  Ransdell  et  al. 
2012;  Tobin  etal.  2009),  a  more  thorough  characterization  of 
the  inward  currents  in  this  species  and  cell  type  is  both 
important  and  valuable  to  moving  research  in  this  area  forward. 
Here  we  utilize  multiple  techniques,  including  two-electrode 
voltage  clamp,  immunohistochemistry,  and  single-cell  PCR 
analysis  of  channel  gene  expression,  to  explore  these  depolar¬ 
izing  components  in  the  LCs  of  C.  borealis. 

MATERIALS  AND  METHODS 

C.  borealis  crabs  were  purchased  and  shipped  overnight  from  The 
Fresh  Lobster  Company  (Gloucester,  MA).  Crabs  were  kept  between 
24  h  and  2  wk  in  artificial  seawater  at  12°C  before  use.  Crabs  were 
anesthetized  in  ice  for  15  min  prior  to  dissection.  The  dissection  took 


place  in  chilled  physiological  saline  comprised  of  (in  mM  )  440  NaCI, 
26  M  gCI2, 13  CaCI2, 11  KCI,  and  10  mM  HEPES  (pH  =  7.4).  Details 
of  the  CG  dissection  can  be  found  in  Cruz-Bermudez  and  M  arder 
(2007).  All  experiments  were  conducted  on  the  three  anterior  LCs  (see 
Fig.  1,  A  and  C).  To  isolate  LC  somata  from  action  potentials  and 
network  activity,  individual  strands  of  bulking  nylon  (made  by  Gut- 
ermann  Creative)  were  used  to  ligate  the  nerve  at  the  anterior  branch 
point  of  the  CG  (Fig.  1).  To  impale  LCs,  each  cell  was  desheathed 
with  a  desheathing  pin.  Anterior  LCs  on  the  same  branch  were 
isolated  from  the  rest  of  the  network  in  a  pair. 

Electrophysiology.  All  experiments  were  done  in  physiological 
saline  cooled  to  12°C.  Two-electrode  voltage  clamp  (TEVC)  and 
two-electrode  current  clamp  (TECC)  protocols  were  carried  out  by 
impaling  a  LC  with  two  glass  electrodes  filled  with  3  M  KCI  (8-  to 
17-MO  resistance)  with  an  Axoclamp  2A  amplifier  (Axon  Instru¬ 
ments,  Union  City,  CA).  In  TEVC  experiments  when  all  potassium 
currents  needed  to  be  blocked,  the  current- injecting  electrode  was 
filled  with  1  M  tetraethylammonium  (TEA)  +  1  M  cesium  chloride 
(CsCI)  (17-  to  18-M  O  resistance).  All  recordings  were  made  from  LC 
somata;  action  potential  conductances  were  blocked  (unless  otherwise 
noted)  by  tightening  athread  ligature  past  the  anterior  branch  point  on 
the  nerve  in  which  the  LC  was  located  (Fig.  L4).  TEVC  and  TECC 
protocols  were  created,  driven,  and  recorded  with  Clampex  9.2  soft¬ 
ware  (M  olecular  Devices).  Current  and  voltage  recordings  were  ana¬ 
lyzed  with  Clampfit  9.2  software  (Molecular  Devices).  Current  and 
voltage  traces  were  filtered  with  a  digital  low-pass  boxcar  filter  using 
seven  smoothing  points  at  a  sample  frequency  of  5  kHz. 

Pharmacology.  Pharmacological  blockers  were  dissolved  in  phys¬ 
iological  saline  and  perfused  onto  the  CG  with  a  Rabbit  peristaltic 


Fig.  1.  The  crustacean  cardiac  ganglion  as  a  model  system.  A:  schematic  showing  the  location  of  the  5  large  cell  (LC)  motoneurons  (black  ovals),  with  the  anterior 
cells  used  in  this  study  marked  with  black  arrows.  The  somata  of  these  anterior  large  cells  were  isolated  from  the  remainder  of  the  ganglion  via  thread  ligatures 
placed  at  the  branch  points  labeled  by  "X ."  Small  cell  (SC)  pacemakers  are  indicated  as  well.  6:  electrophysiological  recording  showing  the  endogenous  network 
activity,  with  a  simultaneous  intracellular  recording  of  a  LC  motoneuron  (top)  and  an  extracellular  recording  of  the  entire  network  output  (bottom).  C :  dye  fills 
(Alexa  Fluor  568)  of  the  5  LC  motoneurons  of  an  intact  cardiac  ganglion,  showing  the  relative  locations  as  well  as  the  clarity  with  which  LC  somata  can  be 
identified  visually. 
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pump  (Rainin  Instrument)  at  a  rate  of  1.5  ml/min  or  added  from  a 
concentrated  stock  solution  via  pipette.  Prior  to  TEVC  or  TECC 
protocols,  perfusion  was  stopped  for  at  least  5  min  to  let  the  prepa¬ 
ration  temperature  stabilize. 

Voltage-clamp  protocols.  To  voltage  clamp  inward  currents,  potas¬ 
sium  currents  were  blocked  by  perfusion  of  50  mM  TEA  and  1  mM 
4-ami  nopyridine  (4-AP)  and  injection  of  1  M  TEA  +  1  M  CsCI  by 
iontophoresis  (300-ms,  1.2-nA  pulses  at  2  Hz)  for  a  minimum  of  1  h.  To 
isolate  calcium  currents,  10~6  M  TTX  also  was  added  for  a  minimum  of 
45  min.  When  necessary,  calcium  currents  were  blocked  with  250  juM 
CdCI2  perfusion  for  a  minimum  of  45  min.  Because  of  difficulty  in  fully 
blocking  outward  currents,  we  often  carried  out  trace  subtractions  before 
and  after  CdCI2  (45  min)  and  before  and  after  TTX  (45  min)  to  isolate 
Cd-sensitive  currents  and  TTX -sensitive  currents,  respectively. 

To  separate  inactivating  and  noninactivating  high-voltage-activated 
calcium  currents,  two  TEVC  protocols  were  carried  out:  one  from  a 
holding  potential  of  -40  mV  with  680-ms  depolarization  steps  from  -50 
mV  to  +15  mV  in  5-mV  increments  (holding  potential  =  -40  mV)  and 
a  second  protocol  identical  except  from  a  holding  potential  of  -80  mV. 
Trace  subtractions  were  used  to  isolate  the  additional  inward  current 
clamped  with  the  -80-mV  holding  potential;  this  current  is  labeled  /CaT, 
while  the  inward  calcium  current  clamped  from  a  -40-mV  holding 
potential  is  labeled  lCaS  (Liu  et  al.  1998;  Prinz  et  al.  2003).  lCaT 
inactivation  was  measured  in  TEVC  with  a  800-ms  prestep  (-80  mV  to 
-40  mV,  5-mV  increments)  prior  to  a  500-ms,  -20-mV  depolarizing 
step.  To  remove  lCaS,  the  current  generated  with  the  -40-mV  prestep 
was  subtracted  from  each  sweep.  lCaS  inactivation  was  measured  in 
TEVC  with  2-s  presteps  (-50  mV  to  +5  mV,  5-mV  increments) 
followed  by  a  500-ms,  0-mV  depolarizing  step. 

Inward  tail  current  was  isolated  with  the  same  pharmacological  agents 
as  voltage-gated  calcium  currents,  and  its  magnitude  was  measured  at 
-80-mV  holding  potential  after  680-ms  depolarizing  presteps  between 
-50  and  +15  mV.  To  isolate /CaS  contribution  to  tail  current,  a  -40-mV 
holding  potential  was  used  prior  to  the  depolarizing  voltage  step;  in  this 
protocol  the  tail  current  was  still  measured  at  -80  mV  (see  Fig.  5B). 
Reversal  of  /CAN  was  measured  after  a  1-s,  +20-mV  depolarizing  voltage 
step,  M  agnitude  of  ZCAN  was  always  measured  at  peak  tail  current  after 
a  depolarizing  voltage  step. 

TTX -sensitive  current  was  isolated  with  the  same  pharmacological 
agents  used  to  isolate  calcium  currents,  with  the  addition  of  250  /xM 
CdCI2.  Depolarizing  step  TEVC  protocols  (steps  =  680  ms,  -50  mV 
to  +15  mV;  holding  potential  =  -40  mV)  and  voltage  ramp  TEVC 
protocols  (4-s  ramp,  -80  mV  to  +20  mV )  were  carried  out  before  and 
after  30-min  exposure  to  10~6  M  TTX,  and  trace  subtractions  were 
used  to  isolate  the  TTX -sensitive  current. 

Sequencing  of  calcium  channels.  We  obtain  de  novo  partial  ORF 
sequence  corresponding  to  three  calcium  channel  a-subunit  cDNAs  to 
aid  in  characterizing  the  channels  potentially  underlying  calcium  currents 
in  LC  motoneurons.  The  methods  for  this  degenerate  PCR  approach  are 
as  described  in  detail  previously  (Schulz  et  al.  2006,  2007).  Briefly, 
channel  transcripts  were  amplified  from  a  cDNA  template  derived  from 
total  RNA  extracted  from  mixed  nervous  system  tissue  from  C.  borealis. 
We  first  designed  and  employed  degenerate  primer  pairs  based  on  con¬ 
served  amino  acid  sequence  compared  across  multiple  invertebrate  and 
vertebrate  species,  including  Lymnaea  stagnalis,  Drosophila  melano- 
gaster,Apis  mellifera,  D  aphnia  pulex,  M  us  musculus,  and  H omo  sapiens. 
Weused  these  primers  in  RT-PCR  reactions  (GoTaq  Green;  Promega)to 
obtain  PCR  products  of  predicted  length,  which  were  subsequently  li¬ 
gated  into  pGEM-T  easy  plasmid  vector  (Promega)  and  sequenced  by 
dye  terminator  cycle  sequencing  (DTCS  Quick  Start  Kit;  Beckman 
Coulter).  Sequences  obtained  were  compared  to  orthologous  sequences 
with  BlastX  (NCBI).  Initial  sequence  obtained  was  extended  by  aprocess 
of  additional  RT-PCR  reactions,  each  using  one  gene-specific  primer 
within  the  newly  elucidated  sequence  and  one  degenerate  primer  outside 
the  known  sequence  that  would  allow  us  to  extend  the  sequence.  Once  we 
obtained  enough  sequence  information  to  feel  that  we  definitively  could 
identify  a  calcium  channel  ortholog,  a  final  PCR  (Platinum  Taq  DNA 


polymerase;  Invitrogen)  was  performed  with  gene-specific  primers  that 
spanned  the  entire  assembled  sequence  to  ensure  that  a  contiguous  cDN  A 
existed  in  our  template  pool.  Accession  numbers  for  these  new  sequences 
are  as  follows  (Cb  =  Cancer  borealis ):  CbCaVl;  JN 809809,  CbCaV2: 
JN809808,  CbCaV3:  J N 809810.  Primer  sequences  that  will  result  in 
PCR  products  that  contain  the  entirety  of  the  assembled  sequence  thus 
far  are  as  follows:  CbCaVl- Forward  (F)  5'-CAATGCATTACGAT- 
GGAGGGG-3',  CbCaVl- Reverse  (R)  5'-AAACCTCCTGAAGTA- 
GTCTTG-3';  CbCaV2-F  5'-TGGAACATTATGGACTTCGTTG-3', 
CbCav2-R  5'-CCAGCCCCACCAAGCCTCCTG-3';  CbCav3-F  5'- 
CAAGGATGGCTGGGTCAACAT-3',  CbCav3-R  5'-GGCTCTAA- 
TGCCCTTGGC  CAT-3'. 

Single-cell  PCR.  Details  regarding  the  harvesting  of  individual  LC 
motoneurons  for  single-cell  molecular  analysis  are  provided  in  a 
previous  publication  (Tobin  et  al.  2009).  Briefly,  ganglia  were 
desheathed  around  the  area  where  LCs  were  visually  identified,  and  to 
facilitate  removal  of  the  neurons  the  ganglia  were  exposed  to  protease 
(Sigma)  for  several  minutes.  The  protease  was  replaced  with  a  cold 
( — 0°C )  solution  of  70%  ethylene  glycol  and  30%  saline,  and  the 
ganglia  were  put  in  a  -20°C  freezer  for  up  to  1  h.  The  cells  were  then 
manually  removed  with  fine  forceps,  and  each  cell  was  placed  in  a 
cryogenic  tube  with  350  ml  of  lysis  buffer  (buffer  RLT ;  Qiagen)  and 
1%  j3-mercaptoethanol.  The  tubes  were  immediately  placed  on  ice  and 
then  stored  at  -80°C  until  RNA  extraction. 

Quantitative  RT-PCR  was  performed  as  in  Schulz  et  al.  (2006)  and 
modified  as  in  Tobin  et  al.  (2009).  Primers  specific  for  real-time  PCR 
detection  of  CbCaVl,  CbCaV2,  and  CbCaV3  with  SYBR  Green  were 
developed  and  designed  with  Primer3  software  as  follows:  CbCaVl-F 
5'-GATGGAGGGAGACGACAAGA-3',  CbCaVl-R  5'-GGTTGTT - 
GAGCCTCGTCATT-3';  CbCaV2-F  5'-GTATCCGGCGGACAGT- 
AAAG-3',  CbCaV2-R  5'-AAACTTGGTGAGAAATGGCG-3'; 
CbCaV 3-F  5'-CCAGTTGGATGTGTTCATCG-3',  CbCaV3-R  5'- 
TGATGGTTGGGTTGATTGG-3'.  PCR  amplicon  sizes  were  90  bp 
(CbCaVl),  120  bp  (CbCaV2),  and  93  bp  (CbCaV3).  Primer  sets  were 
validated  for  PCR  efficiency,  linear  dynamic  range,  and  limit  of  detection 
with  both  a  plasmid  template  that  contained  the  gene  of  interest  and 
cDNA  generated  from  total  RNA  extracted  from  crab  neural  tissue. 
All  qPCR  reactions  were  carried  outwith  a  final  primer  concentration  of 
0.25  /xM  for  each  primer,  using  RT2  SYBR  Green  qPCR  Mastermix 
(SA Biosciences).  RNA  from  single  neurons  was  extracted  with  a  RN  easy 
micro  kit  (Qiagen)  and  reverse  transcription  with  Superscript  III  (Invit¬ 
rogen).  SeeSchulz  etal.  (2006,  2007)  for  more  details.  Assay  quality  and 
variability  were  monitored  via  18S  rRNA  quantitation  (as  in  Ransdell 
etal.  2010;  Schulz  etal.  2006, 2007;  Tobin  etal.  2009).  Data  in  this  study 
are  presented  as  raw  copy  number  per  cell  (see  Tobin  et  al.  2009). 

Cell  fills  and  immunohistochemistry.  To  visualize  stereotypical 
localization  of  LC  somata  and  projections  in  the  CG,  we  iontophore- 
sed  a  fluorescent  tracer  (Alexa  Fluor  568;  Invitrogen)  with  backfilled 
microelectrodes  by  overriding  the  electrode  capacitance  compensation 
for  5  min.  Dye  was  allowed  to  diffuse  throughout  the  ganglia  for 
24  h  in  physiological  saline  at  4°C,  and  then  preparations  were 
directly  imaged  with  a  Leica  M  205F  stereomicroscope  and  a  Leica 
DFC  345  FX  camera. 

For  immunohistochemistry,  CGs  with  desheathed  LCs  were  pinned 
to  Sylgard  and  fixed  in  0.1  M  phosphate-buffered  saline  (PBS;  in  mM  : 
440  NaCI,  11  KCI,  10  Na2HP04,  2  KH2P04,  pH  7.4)  containing  2% 
paraformaldehyde  for  40-90  min  and  then  washed  2  x  10  min  in  PBS 
containing  0.1%  Triton  X  (PBT).  Rabbit  panpolyclonal  primary  anti¬ 
body  to  CaVal  subunit  (Anti-CaVpan  al;  Alomone  Labs)  was  applied 
overnight  at  a  concentration  of  1:100  in  PBT  at  4°C.  Ganglia  were  rinsed 
4  x  10  min.  in  PBT.  CaVal-likeimmunoreactivity  was  investigated  with 
1:250  dilution  goat  anti-rabbit  IgG  conjugated  to  Alexa  Fluor  568  dye  in 
PBT  (M  olecular  Probes/I nvitrogen)  for  12-20  h.  Ganglia  were  rinsed 
again  in  PBS  at  least  5  x  15  min  and  mounted  in  ProLong  Gold  that 
contained  DAPI  (Invitrogen).  I  mages  were  acquired  with  an  Olympus 
1X81  epifi uorescent  microscope  and  M  icroSuite  software. 
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Statistics.  Analyses  were  done  with  SigmaPlot  11.0.  r  Values  reported 
are  results  of  Pearson  correlation  (see  Fig.  31).  After  passing  Shapiro- 
Wilk  tests  for  normality,  changes  in  current  magnitude  before  and  after 
pharmacological  treatments  were  analyzed  with  paired  t-tests. 


RESULTS 

Calcium  channel  expression  in  LCs.  To  determine  which 
calcium  current  subtypes  may  be  present  in  LC  motoneurons, 
we  decided  first  to  adopt  a  molecular  approach  and  look  for 
expression  of  voltage-gated  calcium  channels  at  the  mRNA 
level.  The  most  thorough  combined  molecular  and  el  ectrophys- 
iological  characterization  of  voltage-gated  calcium  channels  in 
invertebrates  has  been  performed  for  the  snail  L.  stagnalis 
(Senatoreand  Spafford  2010;  Spafford  etal.  2003,  2006).  Thus 
we  used  direct  comparisons  of  our  channels  and  currents  to 
these  well-characterized  invertebrate  examples  as  the  basis  for 
interpretation  of  our  molecular  calcium  channel  data.  We 
obtained  de  novo  sequence,  or  extended  partial  sequences,  for 
three  channel  genes  that  appear  to  be  orthologous  to  single 
representatives  for  each  of  three  calcium  channel  gene  families 
CaVl  (L  type),  CaV2  (P/Q  type),  and  CaV3  (T  type)  found  in 
most  invertebrates  (Jeziorski  et  al.  2000;  Spafford  and  Zam- 
poni  2003).  We  refer  to  each  crab  gene  as  CbCaVl,  CbCaV2, 
and  CbCaV3,  Cb  to  denote  C.  borealis  and  1,  2,  or  3  for  the 
gene  family  to  which  the  sequence  is  most  orthologous.  We 
assigned  gene  identities  to  our  channel  transcripts  by  comparing 
our  three  putative  calcium  channel  genes  to  the  orthologous  L. 
stagnalis  calcium  channel  sequences  (Fig.  2 A).  Obtaining  partial 
sequence  of  the  calcium  channel  genes  allowed  us  to  develop 
primers  to  analyze  single  LC  mRNA  expression  of  each  channel 
gene  (Fig.  28).  LC  motoneurons  expressed  detectable  amounts  of 
all  three  putative  calcium  channel  subunits,  with  CbCaVl  having 
the  lowest  mRNA  copy  number  while  higher  levels  of  both 
CbCaV2  and  CbCaV3  were  detected  (Fig.  28). 

Across  the  respective  stretches  of  these  channel  sequences 
between  C.  borealis  and  L.  stagnalis,  CbCaV  1  sequence  shares 
581  of  801  identical  amino  acids  with  L.  stagnalis  (LCaVl; 
73%),  CbCaV2  (initially  referred  to  as  cacophony:  seeTobin  et 
al.  2009)  shares  970  of  1,483  identical  amino  acids  (65%)  to 
LCaV2,  and  CbCaV3  shares  325  of  428  identical  amino  acids 
(76%)  to  LCaV3.  The  locations  of  these  regions  of  overlap 
between  C.  borealis  and  L.  stagnalis  are  shown  in  Fig.  2C  with 
respect  to  putative  transmembrane  domains  and  IQ-like  bind¬ 
ing  motifs  of  these  channels  (Van  Petegem  et  al.  2005).  These 
binding  domains  and  their  amino  acid  similarity  to  mammalian 
orthologs  are  shown  in  Fig.  2D.  In  addition,  to  confirm  their 
identity  as  calcium  channels  we  identified  and  aligned  the 
EEEE  and  EEDD  loci  of  the  pore-forming  loops  with  consen¬ 
sus  mammalian  sequences.  These  EEEE/EEDD  loci  are  known 
to  confer  divalent  cation  selectivity  to  these  channels  (Cens  et 
al.  2007;  Tang  etal.  1993;  Y ang  etal.  1993)  and  area  hallmark 
for  the  identification  of  voltage-gated  calcium  channels. 

Calcium  currents.  To  isolate  voltage-gated  calcium  currents, 
we  minimized  competing  sodium  and  potassium  currents  with 
pharmacological  blockers  (see  materials  and  methods).  Total 
inward  current  under  these  conditions  was  measured  (Fig.  3A). 
This  revealed  what  appeared  to  be  a  faster  transient  inward 
current  (/CaT;  Fig.  3A),  a  slower,  more  persistent  inward  current 
( / c as i  Fig.  3 A),  a  late  outward  current,  and  an  inward  tail 
current  (Fig.  3A).  Addition  of  CdCI2  (see  Fig.  4F)  completely 


eliminates  the  inward  currents  and  reduces  or  eliminates  the 
tail  current,  leaving  behind  only  a  small  residual  outward 
current  that  is  presumably  due  to  an  incomplete  block  of  / A 
and/or /Kd  T aken  together,  these  data  plus  the  gene  expression 
results  suggest  the  presence  of  multiple  distinct  calcium  and 
calcium-dependent  inward  currents  in  the  LC  motoneurons. 

W  e  then  attempted  to  characterize  these  distinct  currents  by 
a  combination  of  pharmacology  and  voltage-clamp  protocol 
manipulation.  Low-threshold,  T-type  calcium  currents  in  L. 
stagnalis  ( L C aV 3 )  activate  at  hyperpolarized  voltages  around 
-60  mV,  are  transient,  and  are  sensitive  to  block  with  nickel 
(Senatore  and  Spafford  2010).  We  were  unable  to  elicit  inward 
current  in  LC  motoneuron  somata  for  any  voltages  lower  than 
-30  mV  (see  Fig.  3D),  and  application  of  N iC l2  (75-250  p-M  , 
n  =  4)  resulted  in  no  change  in  inward  current  magnitude  or 
kinetics  (data  not  shown).  Therefore,  although  we  detected 
relatively  abundant  transcripts  for  CbCaV  3  in  LC  motoneurons 
(Fig.  28),  we  were  unable  to  detect  any  active  T-type  current 
from  somatic  voltage-clamp  experiments  in  these  cells. 

P/Q-type  calcium  currents  (e.g.,  LCav2)  are  more  rapidly 
transient  than  L -type  calcium  currents  (e.g.,  LCavl)  and  inac¬ 
tivate  at  more  hyperpolarized  holding  potentials  (H  i  1 1  e  2001; 
Tsien  et  al.  1988).  Therefore,  we  set  out  to  separate  these 
potentially  distinct  native  high-threshold  calcium  currents  by 
subtracting  inward  current  clamped  with  a  -40-mV  holding 
potential  from  inward  current  with  a  -80  mV -holding  poten¬ 
tial  in  LC  motoneurons.  The  more  persistent  calcium  current 
(I casi  FiQ.  38)  activates  at  -30  mV  and  reaches  its  peak  near 
0  mV  (Fig.  3D).  The  transient  current  (/CaT;  Fig.  3C)  activates 
at  a  voltage  similar  to  /CaS  (i.e.,  -30  mV)  but  reaches  its  peak 
at  more  hyperpolarized  voltages  near  -20  mV  (Fig.  3D). 
Native  current  levels  for  lCaS  and  /CaT  in  individual  LC  mo¬ 
toneurons  were  variable  from  cell  to  cell  and  ranged  from  ~3 
to  6  nA  for  both  currents  (Fig.  3E).  However,  the  most 
noticeable  distinction  between  lCaS  and  /CaT  in  these  cells  is  the 
voltage  dependence  of  inactivation.  lCaT  shows  an  —30  mV 
more  hyperpolarized  half-inactivation  voltage  than  lCaS  (Fig. 
3F).  In  addition,  lCaS  and  /CaT  are  distinguished  by  differences 
in  time  constant  of  inactivation  (Fig.  3G),  and  time  to  peak 
(Fig.  3 H),  as  a  function  of  voltage.  There  is  a  consistent 
relationship  between  time  to  peak  current  for /CaS  and  /CaT  that 
would  suggest  these  currents  are  able  to  maintain  distinct 
impacts  on  cell  voltage  across  a  range  of  activation  voltages 
(Fig.  31).  These  characteristics  are  consistent  with  differences 
reported  between  high-voltage  activated  calcium  currents  LCaVl 
and  LCaV2  in  L.  stagnalis  (Spafford  et  al.  2006)  and  similar  to 
differences  among  L-  and  P/Q-type  calcium  currents  in  mammals 
(Hille  2001;  Tsien  et  al.  1988).  On  the  basis  of  these  results,  we 
propose  that  lCaS  current  in  crab  LC  motoneurons  is  homologous 
to  CaV  1  current  as  described  in  L.  stagnalis  (Spafford  et  al.  2006) 
and  reminiscent  of  mammalian  L -type  calcium  currents  (Hille 
2001;  Tsien  et  al.  1988),  while  lCaT  is  homologous  to  LCaV2  as 
described  in  L.  stagnalis  (Spafford  et  al.  2003,  2006)  and  remi¬ 
niscent  of  mammalian  P/Q-type  currents  (Hille  2001;  J  eziorski  et 
al.  2000;  Tsien  et  al.  1988). 

To  further  confirm  the  characteristics  of  these  calcium  cur¬ 
rents,  as  well  as  investigate  any  calcium-dependent  aspects  of 
their  kinetics,  we  performed  our  voltage-clamp  experiments 
with  B aC 1 2  substituted  for  C aC 1 2  in  the  saline  as  well  as  with 
addition  of  CdCI2  to  the  saline.  Both  lCaS  and  lCaT  are  carried 
by  Ba2+  and  show  activation  characteristics  similar  to  the 
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Fig.  2.  Sequencing  of  voltage-gated  calcium  channel  a-subunits  from  Cancer  borealis.  A:  %  amino  acid  identity  for  sequence  obtained  of  putative  CbCaV  channel 
genes  and  their  orthologs  in  the  snail  Lymnaea  stagnalis.  B:  mean  (±SD)  single  LC  copy  number  for  CbCaVl,  CbCaV2,  and  CbCaV3  mRNAs.  Sample  sizes 
are  shown  in  each  bar.  C:  graphical  amino  acid  alignment  between  full-length  L.  stagnalis  protein  sequence  for  each  voltage-gated  calcium  channel  subtype 
(amino  acid  position  represented  on  x-axes)  and  the  corresponding  portion  identified  in  C.  borealis.  Identical  amino  acids,  as  quantified  in  A,  are  denoted  with 
a  vertical  bar  at  each  amino  acid  position  along  the  x-axis.  Putative  transmembrane  domains  (D I -DIV),  the  IQ-binding  domain  (IQ),  and  the  epitope  to  which 
the  primary  antibody  was  made  (Ab)  are  denoted  above  each  channel  sequence.  D:  amino  acid  alignments  of  the  IQ-  and  IQ-like  binding  domains  of  CbCaVl 
and  CbCav2  with  human  (accession  no.  2 B E 6_ D )  and  rat  (accession  no.  NP_037050.2)  orthologs,  respectively.  Gray  residues  indicate  differences  in  the  crab 
sequence  relative  to  mammal  sequence.  Arrows  denote  key  residues  identified  in  modulating  properties  of  these  channels  with  respect  to  calcium-dependent 
inactivation  in  mammalian  systems  (described  in  text).  E:  sequence  alignments  of  the  pore-forming  loops  of  consensus  mammalian  sequence  from  a-subunits 
for  Cav2.1,  Cavl.2,  and  Cav3.2  with  the  corresponding  C.  borealis  ortholog.  Gray  boxes  highlight  the  position  of  EEEE  (high-voltage-activated  channels)  or 
EEDD  (low-voltage-activated  channels)  loci  in  each  domain  of  the  channel  subunitthat  confer  divalent  cation  selectivity.  Sequence  from  CbCaV3  (inset  box) 
was  sufficient  only  for  comparison  of  domain  III. 


native  calcium  currents  (Fig.  4,  A  and  C).  Furthermore,  both 
/CaS  and  lCaT  are  eliminated  by  C dC 1 2  administration  (Fig.  4F). 
This  further  supports  the  hypothesis  that  these  inward  currents 
are  indeed  calcium  mediated.  We  also  noted  that  the  slow 
inactivation  of  /CaS  is  eliminated  or  significantly  reduced  (P  < 
0.01;  n  =  5)  in  the  corresponding  /Ba  (Fig.  4,  B  and  E), 
suggesting  that  the  inactivation  of  this  current  is  not  purely 
voltage  dependent,  but  rather  a  calcium-dependent  inactivation 
(CDI).  This  is  also  a  characteristic  consistent  with  mammalian 
and  other  invertebrate  L -type  calcium  currents  (Hille  2001). 
We  were  able  to  identify  in  our  CbCaVl  sequence  an  IQ- 
binding  motif  that  has  been  shown  to  strongly  influence  CDI  in 
mammalian  L-type  channels  (Zuhlke  et  al.  2000).  Site-directed 


mutagenesis  shows  that  the  IQ  amino  acid  residues  of  this 
motif,  as  well  as  an  FRK  motif  within  the  domain,  are  signif¬ 
icant  mediators  of  CDI,  with  the  mutagenesis  of  K  to  A  having 
the  least  impact  in  the  motif  (Zuhlke  et  al.  2000).  Our  putative 
IQ  domain  shows  this  same  pattern  of  amino  acid  residues 
(Fig.  2D),  differing  only  in  the  K  residue,  which  is  consistent 
with  the  observation  of  C D I  in  our  /CaS  characterization. 
Additionally,  we  sometimes  detected  a  decrease  in  the  magni¬ 
tude  of  /CaT  following  barium  substitution  (Fig.  4D),  which 
could  be  indicative  of  a  calcium-dependent  facilitation  of  the 
native  transient  calcium  current  in  these  cells  (Tazaki  and 
Cooke  1990).  While  these  collective  data  are  consistent  with 
CDI  of  lCaS  as  suggested  by  barium  substitution,  it  is  possible 
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Fig.  3.  Properties  of  voltage-gated  calcium  currents  in  LC  motoneurons.  A:  raw  trace  of  total  current  elicited  at  0  mV  from  a  holding  potential  of  -80  mV  in 
saline  containing  tetraethylammonium  (TEA),  4-ami  nopyridine  (4-AP),  and  tetrodotoxin  (TTX).  Four  distinct  aspects  of  the  total  current  are  labeled,  including 
"transient"  inward  (/CaT),  "slow"-inactivating  inward  (/CaS),  late  outward  (late  current),  and  inward  tail  (tail  current)  current.  The  voltage  trace  used  to  elicit  the 
current  trace  is  shown  above  the  current  recording.  B:  typical  recording  of  /CaS  at  0  mV  from  a  holding  potential  of  -40  mV  in  saline  containing  TEA,  4-AP, 
TTX,  and  Cs+  loaded  to  reduce  K+  currents  to  the  minimum  possible.  The  voltage  trace  used  to  elicit  the  current  is  shown  above  the  current  recording. 
C:  representative  recording  of  /CaT  at  -20  mV  from  a  holding  potential  of  -80  mV  in  saline  containing  TEA,  4-AP,  TTX,  and  Cs+  loaded  to  reduce  K+  currents 
to  the  minimum  possible.  The  voltage  trace  used  to  elicit  the  current  is  shown  above  the  current  recording;  -20  mV  allowed  us  the  best  recordings  of  /CaT  with 
minimal  contamination  of  /CaS,  as  the  current  shown  is  the  result  of  a  subtraction  of  current  elicited  from  holding  potential  of  -40  mV  (i.e.,  /CaS)  from  the  current 
elicited  from  a  holding  potential  of  -80  mV  (/CaS  +  /CaT).  We  also  provide  representative  traces  of  voltage-dependent  activation  of  this  current  at  -25  mV, 
-20  mV,  and  -10  mV.  D:  mean  (±SD)  current-voltage  (l-V)  relationship  for  /CaS  and  /CaT  indicating  the  voltage  dependence  of  activation.  E:  mean  (±SD) 
current  levels  for  /CaS  and  lCaT  across  LC  motoneurons.  Individual  points  represent  1  cell.  These  data  are  from  the  same  cells  used  to  generate  the  l-V  curves 
in  D  (n  as  shown).  /CaS  measurements  were  taken  at  0  mV  and  /CaT  at  -20  mV.  F:  mean  (±SD)  l-V  relationships  for  voltage  dependence  of  inactivation  of  /CaS 
and  /CaT  as  shown  by  normalized  currents  (///max)  measured  at  0  mV  and  -20  mV,  respectively,  from  a  varying  holding  potential  plotted  on  the  x-axis.  G :  mean 
(±SD)  time  constant  of  inactivation  (rinact)  of  /CaS  and  lCaT  as  a  function  of  voltage.  H:  mean  (±SD)  time  to  peak  current  of  /CaS  and  /CaT  as  a  function  of  voltage. 
/:  correlation  of  time  to  peak  current  of  lCaS  and  lCaT  across  voltages  from  -20  to  +10  mV  (data  are  from  cells  used  in  Fig.  2H).  r  Value  is  the  result  of  a  Pearson 
correlation  test. 
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Fig.  4.  Barium  currents  in  LC  motoneurons.  A:  mean  (±SD)  /scaied- V  relationship  for  lCaS  (n  =  8)  and  the  corresponding  current  in  barium-substituted  saline 
(/ Bas»  n  =  4).  Because  post-Ba2+  currents  often  varied  much  greater  in  magnitude  relative  to  native  Ca2+  currents,  we  plotted  normalized  currents  to  allow  direct 
comparison  of  activation  characteristics.  Scaled  currents  are  normalized  relative  to  the  peak  current  in  each  cell  by  dividing  lCa  or  /Ba  at  a  given  voltage  by  the 
absolute  value  of  /peak  for  each  cell.  B:  typical  recordings  of  /CaS  at  -10  mV  from  a  holding  potential  of  -40  mV  before  (pre-Ba2+;  top)  and  after  (post-Ba2+; 
bottom)  barium  substitution.  Black  arrows  show  the  apparent  loss  or  decrease  of  inactivation  in  this  current  after  barium  substitution;  gray  arrows  show  the  loss 
of  tail  current  after  barium  substitution.  C :  mean  (±SD)  /scaied-l/  relationship  for  /CaT  and  the  corresponding  current  in  barium-substituted  saline  (/BaT;  n  =  5 
for  both).  Currents  were  scaled  as  in  A.  D:  typical  recordings  of  /CaT  at  -10  mV  from  a  holding  potential  of  -80  mV  before  (pre-B a2+;  top)  and  after  (post-Ba2+; 
bottom)  barium  substitution.  Gray  arrows  again  show  the  loss  of  tail  current  after  barium  substitution.  E:  quantification  (mean  ±  SD)  of  changes  in  inactivation 
time  constant  for  /CaT  (left)  and  the  ratio  of  peak  to  residual  /CaS  (right)  following  barium  substitution.  Sample  sizes  are  indicated  in  each  bar.  Significant 
differences  are  reported  for  the  results  of  a  paired  t-test.  F:  effects  of  cadmium  chloride  on  /CaS  and  /CaT.  All  4  traces  taken  from  the  same  cell  before  and  after 
cadmium  exposure;  -40  mV  hold  traces  are  unsubtracted,  and  -80  mV  traces  represent  the  difference  current  before  and  after  cadmium,  as  obtained  by 
subtracting  the  current  from  the  -40  mV  hold  from  the  -80  mV  hold  to  get  the  difference  current.  One  representative  cell  is  shown,  but  this  experiment  was 
performed  n  =  4  times.  Gray  arrow  shows  the  loss  of  tail  current  after  cadmium  substitution. 


that  barium  influences  other  currents  that  may  be  incompletely 
blocked  by  our  experimental  protocol,  altering  the  apparent 
inactivation  kinetics  of  /CaS. 

Barium  did  not  cause  a  change  in  inactivation  kinetics  of 
IcaT  (Fig.  4,  D  and  E),  suggesting  that  this  current  does  not 


show  CDI.  There  is  less  known  about  the  molecular  basis  of 
CDI  and  the  IQ  domain  in  P/Q-type  currents.  However,  CDI  in 
these  channels  has  been  linked  to  an  interaction  with  an  IQ-like 
binding  domain  (including  IM  amino  acid  residues;  see  Fig. 
2D)  with  a  calmodulin  binding  domain  within  CaV2  channels 
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in  mammals  (Lee  et  al.  2003).  While  this  interaction  is  not  yet 
well  characterized,  it  is  noteworthy  that  our  CbCaV2  IQ-like 
domain  lacks  the  IM  motif  implicated  in  CDI,  consistent  with 
our  lack  of  evidence  for  CDI  for  /CaT. 

Calcium-activated  nonselective  cationic  current.  We  ob¬ 
served  a  native  inward  tail  current  (Fig.  3/1)  that  is  eliminated 
after  CdCI2  administration  (Fig.  4 F)  and  significantly  de¬ 
creased  after  B aC 1 2  substitution  (pre-Ba2+  mean  =  -5.7  ±  0.5 
nA;  post-Ba2+  mean  =  -2.7  ±  1.0  nA;  n  =  3;  P  =  0.03)  (Fig. 
4,  8  and  C),  suggesting  the  presence  of  a  calcium-activated 
depolarizing  conductance  in  these  cells.  The  observed  charac¬ 
teristics  of  this  current  are  very  similar  to  those  of  a  calcium- 
activated  nonselective  cationic  current  (/CA n )  found  in  STG 
neurons  of  lobster  and  described  in  detail  by  Zhang  et  al. 
(1995).  We  performed  subsequent  voltage-clamp  and  pharma¬ 
cological  manipulations  to  determine  whether  the  current  we 
observe  in  LC  motoneurons  shares  the  same  kinetics  and 
properties  of  /CAN  as  described  in  lobster  (Zhang  et  al.  1995). 

To  analyze /CAN  in  LC  motoneurons,  we  bath  applied  TTX, 
TEA,  and  4-AP  and  iontophoresed  Cs+  i ntracel I ul arly  as  de¬ 
scribed  in  materials  and  methods.  We  then  elicited  inward 
currents  with  increasing  depolarizing  voltage  steps,  followed 
by  clamping  tail  currents  at  -80  mV.  As  seen  in  Fig.  3A  as 
well  as  Fig.  5A  (inset),  these  protocols  resulted  in  initial  inward 
current  activation,  usually  a  late  outward  current,  and  then  an 
inward  tail  current  measured  at  -80  mV  (see  Fig.  5A,  inset). 


The  tail  current  activated  at  voltages  above -30  mV  (concom¬ 
itant  with  lCa;  see  Fig.  3D),  reaching  a  peak  around  +15  mV 
(Fig.  5A). 

W  e  next  investigated  whether  the  prepulse  voltage  and  thus 
activation  of  distinct  voltage-dependent  calcium  currents  (/CaS 
and  lCaT)  is  correlated  to  tail  current  activation  by  clamping  the 
LC  to  depolarized  potentials  that  activate  both  currents  (hold¬ 
ing  potential  =  -80  mV).  Wefound  thatwhileusing  a  holding 
potential  of  -80  mV  (eliciting  both  lCaS  and  /CaT)  results  in  a 
larger  overall  tail  current  (Fig.  58),  holding  at  -40  mV 
(activating  only  lCaS)  still  is  sufficient  to  activate  the  tail 
current  (Fig.  58),  but  reduced  in  magnitude  ~35%.  These 
results  suggest  that  the  tail  current  is  activated  by  calcium 
influx  via  both  lCaS  and  lCaT. 

To  measure  the  reversal  potential  of  the  tail  current,  inward 
current  was  activated  with  a  depolarizing  step  from  a  holding 
potential  of  -40  mV  up  to  +20  mV  and  then  the  tail  current 
clamped  across  a  range  of  voltages  from  -85  mV  to  -20  mV 
(Fig.  5C).  Through  these  experiments  we  obtained  a  current- 
voltage  (l-V)  curve  based  on  measurements  of  the  peak  tail 
current,  from  which  we  extrapolated  the  reversal  potential  of 
the  tail  current  to  be  approximately  -30  mV  (Fig.  5 D). 

Finally,  we  investigated  whether  the  tail  current  is  sensitive 
to  increases  in  intracellular  calcium  release.  Figure  5E  demon¬ 
strates  that  caffeine  exposure  (bath  applied  for  10  min  at  10 
mM ),  known  to  stimulate  intracellular  calcium  release  (Friel 
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Fig.  5.  Properties  of  calcium-activated  nonselective  cationic 
current  (/CAN)  in  LC  motoneurons.  Recordings  were  made 
in  saline  containing  TEA,  4-AP,  and  TTX  and  cells  loaded 
with  Cs+.  A:  mean  (±SD)  i-V  relationships  for  activation  of 
I  can  measured  at  -80  mV  from  a  varying  depolarizing 
pulse  potential  (x-axis)  to  elicit  calcium  currents  (holding 
potential  of  -80  mV,  n  =  5).  Inset,  typical  recordings  of 
the  total  current,  including  tail  current  (black  arrow)  elicited 
from  these  voltage  protocols.  B:  mean  (±SD)  l-V  relation¬ 
ships  for  activation  of  /CAN  from  holding  potentials  of  -40 
mV,  activating  predominantly  /CaS,  or  -80  mV,  which 
activates  both  /CaS  and  lCaT  (n  =  3).  Inset:  representative  tail 
currents  in  the  same  cell  elicited  from  a  holding  potential  of 
-40  mV  (black  trace)  or  -80  mV  (gray  trace).  C:  typical 
recordings  used  to  measure  the  reversal  potential  of  tail  cur¬ 
rent.  Current  traces  were  obtained  by  depolarizing  the  LC 
neuron  from  a  holding  potential  of  -50  mV  to  +20  mV  and 
then  measuring  tail  currents  at  varying  tail  holding  potentials 
from  -80  mV  to  -20  mV  (as  shown  in  voltage  traces). 
D:  mean  (±SD)  l-V  relationships  used  to  interpolate  the  re¬ 
versal  potential  of  these  tail  currents  (n  =  8).  E:  caffeine- 
evoked  increase  in  tail  current  amplitude  in  LCs.  Recording  of 
tail  currents  from  a  holding  potential  of  -80  mV  before  (black 
trace/arrow)  and  after  (gray  trace/arrow)  bath  application  of  10 
mM  caffeine  (this  experiment  was  performed  n  =  4  times). 
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and  Tsien  1992a,  1992b;  Yoshimura  2005;  Zhang  and  Harris- 
Warrick  1995)  as  well  as  an  increased  level  of  calcium-induced 
calcium  release  (Levi  etal.  2003;  M  ironov  and  Usachev  1991), 
causes  an  increase  in  the  inward  tail  current  (n  =  4,  P  = 
0.047).  These  results  suggest  that  the  tail  current  is  indeed 
calcium  activated  and  not  solely  dependent  on  extracellular 
calcium  influx  via  voltage-gated  calcium  currents.  Taken  to¬ 
gether,  these  results  are  very  similar  to  the  characteristics  of  the 
/CAN  in  lobster,  including  sensitivity  to  caffeine  exposure 
(Zhang  and  Harris- Warrick  1995),  indicating  that  the  same 
current  is  indeed  present  in  LC  motoneurons  of  the  crab. 

TTX  application  reveals  a  largely  noninactivating  TTX-sensitive 
current  in  LC  motoneurons.  During  our  experiments,  we  observed 
that  in  the  presence  of  TEA  and  4-AP  C dC 1 2  alone  does  not 
consistently  prevent  all  voltage-activated  inward  current 
(Fig.  6A),  but  rather  CdCI2  and  TTX  combined  are  required  to 
eliminate  inward  current  (Fig.  6 A).  These  observations  suggest 
that  there  is  a  TTX -sensitive  current,  distinctfrom  spike-mediated 
sodium  currents,  involved  in  bursting  in  these  cells.  We  attempted 
to  voltage  clamp  this  current  by  subtracting  the  TTX-sensitive 
inward  current  in  both  step  and  ramp  voltage-damp  protocols, 
with  inconsistent  results.  Because  of  the  small  magnitude  of  the 
TTX-sensitive  current  (on  the  order  of  2  nA  or  smaller,  clamped 
from  the  soma),  voltage-clamp  protocols  using  step  depolariza¬ 
tions  often  were  not  adequate  to  clearly  characterize  this  current 
above  background  noise,  though  they  did  suggest  that  the  current 
is  largely  noninactivating  (Fig.  66).  Voltage  ramp  protocols  used 
before  and  after  TTX  (10  6  M ,  30  min)  to  isolate  the  TTX -sensitive 
current  during  a  -80  to  +20  mV  ramp  demonstrated  that  TTX- 
sensitive  current  activation  began  at  -37.5  mV  ±  4.7  mV  and  peak 
activation  occurred  at  6  mV  ±  3.2  mV  (Fig.  66;  n  =  3). 

Because  we  had  mixed  success  voltage-clamping  the  TTX- 
sensitive  current,  and  because  of  its  often  smaller  magnitude,  we 
set  out  to  determine  whether  this  current  plays  an  active  role  in  the 
burst  dynamics  of  LC  motoneurons.  To  investigate  this,  we  in¬ 
duced  exaggerated  burst  potentials  in  isolated  LCs  by  adding  50 
mM  TEA  to  the  saline  and  generating  large  burst  potentials  with 
depolarizing  current  pulses  (Fig.  6C).  These  potentials  are  intrin¬ 
sically  driven  because  the  depolarizing  potential  lasts  longer  than 
the  brief  current  injection.  After  burst  potentials  were  successfully 
generated,  10~6  M  TTX  was  added,  thus  eliminating  the  TTX- 
sensitive  inward  current.  In  six  of  six  experiments,  TTX  admin¬ 
istration  completely  eliminated  these  burst  potentials  (Fig.  6C), 
supporting  the  presence  of  this  current  and  suggesting  an  active 
role  for  this  current  in  burst  output.  However,  the  exactimpact  of 
this  TTX-sensitive  current  on  burst  dynamics  in  these  cells  is 
beyond  the  scope  of  this  study.  Nevertheless,  it  does  appear  that 
despite  its  small  magnitude,  this  current  is  present  and  plays  a  role 
in  the  output  of  these  cells. 

H yperpolarization-activated  mixed  cationic  current.  We  also  em¬ 
ployed  hyperpolarized  voltage  steps  in  voltage  clamp,  as  well  as 
negative  current  inject  in  current  clamp,  to  investigate  the  pres¬ 
ence  of  hyperpolarization-activated  mixed  cationic  current  (/ H )  in 
the  somata  of  these  cells.  Although  relatively  abundant  mRNAs 
for  the  HCN  channel  have  been  detected  in  these  cells  (Tobin  et 
al.  2009),  we  did  not  see  any  evidence  for  the  H  current  in  the 
somata  of  these  cells  (data  not  shown)  for  voltages  as  low  as 
-120  mV  from  a  holding  potential  of  -40  mV. 

Immunohistochemical  localization  of  voltage-gated  calcium 
channels  in  LC  motoneurons.  The  crustacean  CG  LC  burst 
potential  is  mediated  by  intrinsic  conductances  at  least  in  part 


c 
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Fig.  6.TTX-sensitivepersistentcurrentin  LC  motoneurons.  A:  l-V  relationship 
in  a  LC  motoneuron  indicating  the  presence  of  a  TTX -sensitive  current.  Inward 
currents  were  measured  from  a  holding  potential  of  -40  mV  at  depolarizing 
voltage  steps  under  subsequent  blocker  conditions.  Currents  first  were  elicited 
in  TEA  +  4-AP,  and  then  Cd2+  was  added  to  block  calcium  currents,  revealing 
a  residual  inward  current  that  was  subsequently  blocked  with  TTX.  B:  rep¬ 
resentative  recording  of  TTX-sensitive  current  elicited  in  saline  containing 
TEA,  4-AP,  and  Cd2+  with  a  -80  mV  to  +20  mV  voltage  ramp;  then 
subtracted  from  this  current  was  the  same  voltage  ramp  with  the  addition  of 
TTX  (top;  n  =  3).  Step  depolarizations  also  reveal  a  persistent  TTX  sensitive 
inward  current.  Gray  trace,  inward  current  in  the  presence  of  CdCI2;  black  trace, 
current  after  the  addition  of  1(T6  M  TTX  ( bottom ;  n  =  4).  C:  TTX-sensitive 
current  is  necessary  for  burst  potential  generation  in  LC  motoneurons.  Exaggerated 
driver  potentials  are  elicited  in  50  mM  TEA  and  a  short  depolarizing  current 
injection  (5  nA,  50  ms;  top).  These  driver  potentials  are  eliminated  in  the  presence 
of  TTX,  even  though  calcium  currents  remain  intact  (bottom),  n  =  6. 

isolated  to  the  LC  somata.  In  many  crustacean  species  this 
intrinsic  burst,  termed  a  driver  potential,  can  be  regularly 
elicited  in  an  isolated  LC  somata  with  a  short  depolarizing 
current  injection  (Cooke  2002;  Tazaki  and  Cooke  1983).  How¬ 
ever,  in  C.  borealis  CG  the  isolated  LC  soma  will  rarely  fire  an 
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intrinsic  burst  without  the  application  of  potassium  current 
blockers  such  as  TEA.  While  outward  currents  are  clearly 
present  (Ransdell  et  al.  2012),  our  voltage  clamp  showed 
relatively  low  calcium  current  levels  in  the  LC  somata.  This 
could  be  due  to  the  fact  that  the  calcium  channels  are  localized 
more  distant  to  the  soma  or  reflect  a  low  overall  calcium 
current  in  the  soma  itself.  Therefore,  we  sought  to  find  whether 
and  to  what  degree  calcium  channels  are  localized  to  the  LC 
somata.  We  used  a  panspecies  primary  antibody  selective  for 
high-voltage-activated  calcium  channels  (i.e.,  L  type  and  P/Q 
type)  combined  with  fluorescence  microscopy  to  investigate 
the  distribution  of  CbCaV  l/CbCaV2  channels  in  the  CG.  The 
19-amino  acid  residue  epitope  against  which  the  primary  an¬ 
tibody  was  raised  is  a  100%  amino  acid  identity  match  to 
CbCaV  1  and  matches  17  of  19  amino  acids  in  CbCaV 2  (Fig.  7; 
see  also  Fig.  1C).  No  similarity  of  the  epitope  to  CbCaV3  or 
full-length  LCaV3  was  detected  via  BLASTP  comparison. 

LC  somata  and  their  projections  are  stereotypical ly  distributed 
and  clearly  identifiable  in  this  preparation  (see  Fig.  1C),  making 
rudimentary  subcellular  localization  fairly  straightforward.  We 
found  clear  punctate  staining  localized  to  the  LC  somata  and  to  a 


lesser  degree  staining  in  the  neurites  extending  from  the  LC 
somata  (Fig.  7).  In  all,  nine  whole-mount  CG  preparations  were 
stained,  and  in  each  preparation  we  were  able  to  detect  punctate 
staining  localized  to  the  LC  somata.  Additionally,  we  often  de¬ 
tected  staining  in  the  processes  of  the  anterior  branch  point  of  the 
ganglion  (Fig.  7C),  where  synaptic  connectivity  seems  to  be  most 
prevalent  among  the  cells  of  the  C.  borealis  CG.  We  performed 
additional  DA  PI  staining  within  these  preparations  to  determine 
whether  the  punctate  immunoreactivity  on  the  somata  may  be  the 
result  of  glia  or  other  support  cells  positively  staining.  Immuno¬ 
fluorescence  +  DA  PI  overlays  showed  no  evidence  of  overlap¬ 
ping  fluorescence  to  indicate thiswas the case(e.g.,  seeFig.  7D2). 
We  also  performed  secondary  antibody-only  incubations  that 
showed  no  significant  fluorescence  indicative  of  nonspecific  sig¬ 
nal  attributable  to  the  secondary  antibody.  Finally,  preabsorption 
of  the  primary  antibody  with  the  antigen  peptide  eliminated  im- 
munolabeling  as  well. 

DISCUSSION 

The  crustacean  LC  is  a  valuable  model  for  studying  mo¬ 
toneuron  burst  organization.  The  robust  nature  of  the  LC 
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Fig.  7.  High-voltage-activated  calcium  chan¬ 
nels  (HVACC)  are  localized  to  the  somata  of 
LC  motoneurons.  We  used  an  Anti-CaV 
Pan-al  primary  antibody  to  detect  putative 
HVACC  immunoreactivity  in  the  cardiac  gan¬ 
glion.  The  alignment  of  the  antibody  epitope 
with  the  corresponding  amino  acid  sequence  in 
CbCaV  1  and  CbCaV2  is  shown  at  top.  The 
entire  cardiac  ganglion  is  shown  based  on  as¬ 
sembled  tiled  images  taken  at  x4  magnifica¬ 
tion.  Insets:  magnified  images  of  the  corre¬ 
sponding  portion  of  the  ganglion  marked  with 
dashed  boxes.  We  detected  immunoreactivity 
for  HVACC  in  the  somata  of  all  5  LC  mo¬ 
toneurons.  A  and  B:  HVACC  immunoreactiv¬ 
ity  in  anterior  LCs  appears  to  be  localized  to 
the  somata,  consistent  with  the  ability  to  mea¬ 
sure  these  currents  in  isolated  LC  somata. 
C :  immunoreactivity  was  also  detected  in  pro¬ 
cesses  of  neurons  in  the  cardiac  ganglion,  par¬ 
ticularly  in  the  anterior  branch  point.  D:  the 
paired  posterior  LC  neurons  also  showed  im¬ 
munoreactivity  apparently  localized  to  the  so¬ 
mata  ( D1 ),  which  is  not  the  result  of  glial  or 
support  cell  staining  as  revealed  by  colocaliza¬ 
tion  of  immunoreactivity  (red)  and  DA  PI 
counterstain  (blue)  (D2).  Scale  bars  in  insets 
all  represent  200  ju,m. 
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somata  in  electrophysiology  experiments  and  segregation  from 
action  potential  conductances  as  well  as  other  network  com¬ 
ponents  make  these  cells  ideal  in  experiments  (Cooke  2002; 
Hartline  1967).  However,  a  detailed  understanding  of  the  ionic 
conductances  and  their  underlying  ion  channels,  particularly 
depolarizing  conductances,  is  lacking  in  these  cells.  Here  we 
report  the  results  of  sequence  from  three  voltage-gated  calcium 
channel  a-subunit  genes  in  the  crab  C.  borealis  and  then 
characterize  their  expression  as  well  as  the  depolarizing  cur¬ 
rents  present  in  LC  motoneurons  of  the  CG.  While  we  detect 
mRNA  for  three  different  channel  types  in  individual  LCs,  we 
measured  and  can  distinguish  only  two  distinct  voltage-acti¬ 
vated  calcium  currents  in  LC  somata.  We  additionally  identi¬ 
fied  a  calcium-activated  nonspecific  cationic  (CAN)  current  as 
well  as  a  TTX -sensitive  persistent  inward  current  in  these  cells. 
Finally,  we  performed  immunostaining  that  demonstrates  that 
the  high-voltage-activated  calcium  channels  described  are  at 
least  partially  localized  to  the  LC  somata. 

Sequence  analysis  reveals  that  the  three  calcium  channel 
a-subunits  in  this  study,  namely,  CbCaVl,  CbCaV2,  and 
CbCaV3,  correspond  very  closely  to  identified  calcium  channel 
subunits  of  the  snail  L.  stagnalis  (LCaVl,  LCav2,  and  LCaV3, 
respectively).  Thorough  sequence  and  expression  analyses  of 
these  snail  channels  have  been  performed  and  the  resulting 
current  properties  and  kinetics  of  each  channel  reported  (Sena- 
tore  and  Spafford  2010;  Spafford  et  al.  2003,  2006).  On  the 
basis  of  the  extensive  similarities  between  these  expressed 
snail  currents  and  the  native  crab  currents  characterized  in  our 
study,  we  are  confident  that  we  can  use  this  analysis  to  infer  the 
identity  of  these  channels  and  currents.  The  characteristics  of 
crab  LC  /CaS,  a  slowly  inactivating  high-voltage-activated 
calcium  current,  are  very  similar  to  those  of  the  current  carried 
by  the  LCaVl  channel  (Spafford  et  al.  2006),  which  has  been 
classified  as  an  invertebrate  L-type  calcium  channel  and  cur¬ 
rent.  Therefore,  based  on  the  similarities  in  sequence  between 
CbCaVl  and  LCaVl,  and  the  similarities  between  crab  /CaS 
and  current  carried  by  LCaVl,  we  suggest  that  CbCaVl 
encodes  an  L-type  calcium  current  with  the  properties  de¬ 
scribed  here  as  the  slow  calcium  current  /CaS.  Similarly,  the 
characteristics  of  crab  lCaT,  a  transient  inactivating  high-volt- 
age-activated  calcium  current,  are  very  similar  to  those  of  the 
current  carried  by  the  LCaV2  channel  (Spafford  et  al.  2003), 
which  has  been  classified  as  an  invertebrate  calcium  channel 
with  biophysical  characteristics  found  in  mammalian  P/Q-  and 
N-type  channels.  On  the  basis  of  the  similarities  in  sequence 
between  CbCaV2  and  LCaV2,  and  the  similarities  between 
crab  /CaT  and  current  carried  by  LCaV2,  we  suggest  that 
CbCaV2  encodes  a  P/Q-type  calcium  current  with  the  proper¬ 
ties  described  here  as  the  transient  calcium  current  /CaT. 

Both  of  the  two  voltage-gated  calcium  currents  (/CaS  and 
/CaT)  were  present  in  experiments  where  the  LC  soma  was 
isolated,  indicating  that  these  channels  are  present  in  the  LC 
somata.  This  result  is  supported  by  immunohistochemistry 
experiments  in  which  high-voltage-activated  calcium  channel 
antibody  was  used  to  stain  the  entire  CG,  resulting  in  immu¬ 
nostaining  largely  localized  to  the  LC  soma,  although  some 
staining  also  could  be  seen  on  neuronal  processes  in  different 
regions  of  the  ganglia.  These  results  are  also  consistent  with 
previous  immunolocalization  of  calcium  channel  subunits  in 
the  lobster  STG;  CalA  (P/Q  type)  and  CalD  (L  type)  immu- 


noreactivity  both  were  localized  to  STG  somata,  as  well  as 
primary  neurites  (French  et  al.  2002). 

Calcium  currents  previously  have  been  studied  in  crustacean 
motoneurons,  particularly  in  the  STG.  Total  calcium  current  has 
been  characterized  in  STG  neurons  of  two  species  of  lobster 
(Johnson  et  al.  2003;  Turrigiano  et  al.  1995;  Zhang  and  Harris- 
Warwick  1995)  as  well  as  crab  (Golowasch  and  Marder  1992; 
Hurley  and  Graubard  1998).  However,  across  these  studies  there 
are  species-  and  cell-specific  differences  in  the  calcium  current 
types  reported.  For  example,  Hurley  and  Graubard  (1998)  specif¬ 
ically  searched  for,  butfailed  to  find,  two  distinct  calcium  currents 
with  differing  kinetics  in  the  somata  of  STG  neurons  of  the  crab 
Cancer  productus.  However,  Turrigiano  et  al.  (1995)  measured 
calcium  currents  in  isolated  STG  neurons  of  the  lobster  Panu//rus 
interrupts  in  culture  and  noted  both  transient  and  persistent 
components  to  these  currents,  ultimately  deciding  to  model  cal¬ 
cium  current  in  these  cells  as  two  separate  components.  M  any  of 
these  discrepancies  also  may  be  explained  by  the  fact  that  these 
inward  currents  have  been,  and  continue  to  be,  difficult  to  char¬ 
acterize  in  voltage  clamp  (Golowasch  and  Marder  1992;  Turri¬ 
giano  et  al.  1995).  Here,  with  the  additional  insight  of  molecular 
biology  and  a  scaffolding  of  comparison  provided  by  work  in  L. 
stagnalis,  we  believe  we  have  made  substantial  progress  toward  a 
better  understanding  of  these  complex  native  currents  in  crusta¬ 
cean  motoneurons. 

Beyond  the  voltage-activated  calcium  currents,  LCs  also 
possess  a  calcium-activated  inward  current  represented  by 
/CAN.  The  calcium-activated  nonspecific  cationic  current  is 
found  in  a  variety  of  excitable  cell  types  (Partridge  and  Swan- 
dulla  1988)  but  is  often  implicated  in  neuronal  bursting  (Rubin 
etal.  2009;  Swandulla  and  Lux  1985;  Wilson  etal.  1998).  This 
current  in  invertebrates  has  been  best  characterized  in  crab 
STG  (dorsal  gastric)  neurons.  In  these  cells,  the  current  is 
important  in  coordinating  bistability  of  output;  voltage-gated 
calcium  currents  allow  calcium  influx  that  activates  /CAN, 
inducing  plateau  potentials  indicative  of  the  output  of  these 
cells  (Zhang  and  Harris-Warrick  1995;  Zhang  et  al.  1995).  LC 
motoneurons  of  the  CG  do  not  undergo  plateau  potentials 
during  normal  network  output,  so  this  is  unlikely  to  be  the  role 
of  this  current  in  the  cardiac  output.  W  hile  we  do  not  yet  know 
the  role  of  /CAN  in  the  cardiac  motoneurons,  a  model  for  similar 
rhythm  generation  in  mammalian  pre-B otzi nger  complex  in¬ 
vokes  /CAN  as  vital  in  triggering  networkwide  bursts  after 
recurrent  synaptic  excitation— a  role  that  may  have  parallels  to 
/CAn  in  LCs  (Del  Negro  and  Hayes  2008). 

Finally,  we  also  report  for  the  first  time  in  crab  LCs  a  largely 
noninactivating  TTX -sensitive  current  that  is  necessary  for 
driver  potential  generation  in  these  cells,  which  could  indicate 
the  presence  of  a  persistent  sodium  current  (/NaP).  /NaP  are 
often  implicated  in  burst  potential  generation.  In  mammalian 
pre-B otzinger  neurons,  /NaP  presence  is  known  to  be  important 
in  burst  generation  and  its  relationship  with  leak  conductances 
is  thought  to  dictate  whether  the  neuron  spontaneously  bursts 
or  requires  excitatory  input  (Del  Negro  et  al.  2002a,  2002b). 
There  is  evidence  that  DG  neurons  of  the  STG  contain  /NaP 
important  to  plateau  potential  generation  (Elson  and  Selverston 
1997),  and  a  TTX  -sensitive  persistent  current  has  been  voltage 
clamped  in  cultured  STG  cells  (Turrigiano  et  al.  1995).  While 
there  is  no  channel  gene  known  to  directly  mediate  this  current 
(indeed,  there  is  only  one  voltage-gated  sodium  channel  known 
in  invertebrates;  Loughney  etal.  1989;  Olson  etal.  2008),  there 
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is  the  possibility  that  the  / NaP  is  mediated  by  the  same  channel 
population  as  the  fast- inactivating  sodium  current  as  seen  in  rat 
tuberomammillary  neurons  (Taddese  and  Bean  2002).  How¬ 
ever,  because  there  are  no  spike-mediating  sodium  currents  in 
these  crustacean  motoneuron  somata  (Golowasch  and  M  arder 
1992),  it  seems  unlikely  that  these  persistent  sodium  currents 
represent  a  window  in  fast-spiking  sodium  currents.  Rather,  it 
seems  more  likely  that  a  distinct  channel  isoform  could  medi¬ 
ate  these  currents,  perhaps  as  a  result  of  a  splice  variant  (Lin  et 
al.  2009)  or  RNA  editing  event  (Liu  et  al.  2004)  of  the 
previously  identified  CbNaV  channel  (Dai  etal.  2010). 

The  specific  roles  of  the  inward  currents  in  these  cells, 
especially  those  that  are  calcium  mediated,  are  still  to  be 
experimentally  determined.  The  crustacean  cardiac  motoneu¬ 
rons  share  extensive  similarities  with  other  excitatory  central 
pattern  generator  circuits.  For  example,  >90%  of  rat  pre- 
Botzinger  complex  neurons  elicit  burst  potentials  only  when 
driven  by  excitation,  and  like  LC  motoneurons  these  neurons 
also  express  /CAN,  /NaP,  and  intrinsic  calcium  currents  (Del 
Negro  et  al.  2002b;  M  orgado-V  al  I  e  et  al.  2008;  Pena  et  al. 
2004).  Whether  the  CG  network  LC  burst  potential  is  a  result 
of  summed  excitatory  postsynaptic  potentials  (EPSPs)  or  re¬ 
quires  some  or  all  of  these  intrinsic  inward  components  to 
generate  its  intrinsic  burst  potential  is  currently  under  investi¬ 
gation  (Ball  et  al.  2010;  Franklin  et  al.  2010). 

The  depolarizing  components  found  intheLC  bursting  neurons 
are  prevalent  across  animal  phyla,  especially  in  those  neurons  that 
are  also  bursters.  The  experimental  advantages  of  the  crustacean 
CG  and  the  simplicity  of  the  network  make  it  an  ideal  system  to 
study  how  these  components  interact  with  each  other  and  other 
network  characteristics  to  shape  and  maintain  output,  in  charac¬ 
terizing  the  inward  components  of  C.  borealis  LC  motoneuron, 
we  obtain  critical  information  to  now  move  forward  in  investi¬ 
gating  how  the  LC  intrinsic  components  function  within  the  rest 
of  the  cardiac  network,  and  how  this  model  can  be  applied  broadly 
to  understand  neural  network  function  in  general. 
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Biological  and  theoretical  evidence  suggest  that  individual  neurons  may  achieve  similar  outputs  by  differentially  balancing  variable 
underlying  ionic  conductances.  Despite  the  substantial  amount  of  data  consistent  with  this  idea,  a  direct  biological  demonstration  that 
cells  with  conserved  output,  particularly  within  the  same  network,  achieve  these  outputs  via  different  solutions  has  been  difficult  to 
achieve.  Here  we  demonstrate  definitively  that  neurons  from  native  neural  networks  with  highly  similar  output  achieve  this  conserved 
output  by  differentially  tuning  underlying  conductance  magnitudes.  Multiple  motor  neurons  of  the  crab  ( Cancer  borealis)  cardiac 
ganglion  have  highly  conserved  output  within  a  preparation,  despite  showing  a  2-  4-fold  range  of  conductance  magnitudes.  By  blocking 
subsets  of  these  currents,  we  demonstrate  that  the  remaining  conductances  become  unbalanced,  causing  disparate  output  as  a  result. 
Therefore,  as  strategies  to  understand  neuronal  excitability  become  increasingly  sophisticated,  it  is  important  that  such  variability  in 
excitability  of  neurons,  even  among  those  within  the  same  individual,  is  taken  into  account. 


Introduction 

Neurons  with  similar  and  conserved  output  across  individuals 
can  exhibit  2- 6-fold  variability  in  their  intrinsic  conductances 
(Golowasch  et  al.,  2002;  Swensen  and  Bean,  2005;  Schulz  et  ah, 
2006;  Goaillard  et  ah,  2009;  Ransdell  et  ah,  2012)  as  well  as  in  ion 
channel  mRNA  levels  (Schulz  et  ah,  2006, 2007).  This  observation 
has  motivated  several  computational  studies  that  used  popula¬ 
tions  of  model  cells  with  variable  intrinsic  properties  and  con¬ 
strained  output  features  to  demonstrate  that  model  neurons 
could  achieve  similar  voltage  output  with  disparate  intrinsic  pa¬ 
rameters  (Goldman  et  ah,  2001;  Prinz  et  ah,  2004;  Taylor  et  ah, 
2009;  Ball  et  ah,  2010).  This  approach  has  complemented  and 
extended  accumulating  biological  evidence  that  is  consistent  with 
the  idea  that  output  features  across  identified  cells  may  be  con¬ 
served  by  both  balance  and  compensation  among  key  intrinsic 
conductances,  rather  than  by  tightly  regulating  the  magnitude  of 
individual  conductances  (MacLean  et  ah,  2003;  Swensen  and 
Bean,  2005;  Grashow  et  ah,  2010;  Ransdell  et  ah,  2012).  However, 
most  of  these  studies  involve  compensation  that  leads  to  con¬ 
servation  of  output  following  perturbation  of  the  cellular  ex¬ 
citability  (MacLean  et  ah,  2003;  Bergquist  et  ah,  2010),  often 
via  gene  knock-out  (Van  Wart  and  Matthews,  2006;  Nerbonne 
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et  ah,  2008)  or  pharmacological  manipulation  (Swensen  and 
Bean,  2005;  Ransdell  et  ah,  2012).  Although  these  studies  dem¬ 
onstrate  that  variability  could  arise  as  a  result  of  such  homeo¬ 
static  plasticity,  these  studies  do  not  address  intrinsic 
variability  found  in  these  neurons  with  all  conductances  in¬ 
tact,  or  whether  this  variability  represents  differential  tuning 
to  conserved  output  in  vivo. 

We  use  the  crustacean  cardiac  ganglion,  which  contains 
multiple  large  cell  (LC)  motor  neurons  that  produce  synchro¬ 
nous  bursting  output,  to  test  whether  individual  neurons 
achieve  conserved  output  by  differentially  balancing  conduc¬ 
tance  magnitudes.  We  examine  multiple  motor  neurons 
within  the  same  ganglia,  which  are  genetically  identical  and 
have  similar  output  in  the  native  network,  and  find  that  they 
have  varying  underlying  conductances.  Then  by  disrupting  a 
select  subset  of  conductances  in  these  cells  with  pharmacolog¬ 
ical  blockade,  we  demonstrate  a  corresponding  disruption  in 
conserved  output.  We  believe  this  is  the  first  direct  demon¬ 
stration  that  multiple  neurons,  particularly  within  the  same 
neural  network,  naturally  generate  similar  output  by  differen¬ 
tially  balancing  variable  underlying  conductances,  indepen¬ 
dent  of  induced  compensatory  changes. 

Materials  and  Methods 

Animals.  Crabs  (Cancer  borealis)  of  either  sex  were  purchased  and 
shipped  overnight  from  The  Fresh  Lobster  Company.  Crabs  were  kept  at 
12°C  in  artificial  seawater.  Dissections  took  place  in  chilled  physiological 
saline:  440  mM  NaCl,  26  mM  MgC12  13  mM  CaC12, 1 1  mM  KC1,  and  10  mM 
HEPES,  pH  7.4.  Chemicals  were  obtained  from  Fisher  Scientific  unless 
otherwise  noted. 

Electrophysiology.  All  intracellular  recordings  were  made  from  LC3,  4, 
and  5  somata.  LC  structure  consists  of  a  cell  soma  and  proximal  neurite 
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Figure  1.  A,  Schematic  of  C.  borealis  cardiac  ganglion  showing  location  of  LC  somata  (ovals)  and  pacemaker  interneurons  (small  circles).  Extracellular  recordings  are  taken  from  the  central  nerve 
(dashed  circle).  All  experiments  were  conducted  on  three  anterior  LCs  3-5  (colored  ovals).  These  cells  project  a  proximal  neurite  from  their  somata  toward  a  common  spike  initiation  zone  region 
(green  dashed  box),  and  then  axons  continue  posterior.  8,  Activity  of  LC  neurons  during  spontaneous  network  activity.  Extracellular  recording  (bottom  trace)  shows  LC  (tall  spike)  and  pacemaker 
(short  spike)  APs.  Inset,  Red  box  shows  enlarged  recordings  displaying  that  each  pacemaker  AP  (small  spike)  results  in  an  EPSP  in  all  three  anterior  LCs.  Two  pacemaker  APs  and  resulting  EPSPs  are 
labeled  (arrows).  C,  Variability  in  LC  inward  and  outward  current  magnitudes.  0,  /HTK  variability  across  LCs  within  the  same  network.  Lines  connect  LCs  from  same  individual.  Inset,  /HTK  is  blocked  by 
25  mM  TEA  exposure  (30  min).  Current  traces  before  (black)  and  after  (gray)  TEA  were  clamped  at  0  mV  step  from  —40  mV  holding  potential,  f,  TEA  exposure  in  the  intact  network  desynchronizes 
membrane  potential  waveforms  of  LCs  during  ongoing  rhythmic  activity.  F,  Spiking  output  (extracellular  recordings)  is  affected  by  TEA  applied  to  LC4  and  LC5  somata.  Dashed  line  for  comparison 
of  spike  amplitudes.  Traces  were  taken  from  the  same  recording  well  within  10  min  of  one  another. 


that  projects  to  a  common  spike  initiation  area  located  posterior  to  the 
anterior  branch  point  (Hartline,  1967;  Tazaki  and  Cooke,  1983;  Cooke, 
2002;  Ransdell  et  al.,  2013)  (Fig.  1A).  Somata  were  separated  from  action 
potential  (AP)  conductances  (unless  noted  otherwise)  by  tightening  a 
thread  ligature  past  the  anterior  branch  point  on  the  nerve  the  LC  was 
located  in  (Figs.  1A,  2A,B).  Isolated  somata  carry  ionic  conductances 
that  generate  driver  potential  output  (Tazaki,  1972;  Tazaki  and  Cooke, 
1983, 1986;  Cooke,  2002),  even  with  ligatures  as  close  as  200  fxm  from  the 
soma  (Tazaki  and  Cooke,  1983). 

Experiments  were  performed  in  physiological  saline  at  12°C.  Extracel¬ 
lular  recordings  were  made  with  stainless  steel  pin  electrodes  in  petro¬ 
leum  jelly  wells  placed  on  the  nerve.  Signals  were  amplified  and  filtered 
using  a  differential  AC  amplifier  (A-M  Systems).  Intracellular  recordings 
were  made  with  glass  electrodes  containing  3  M  KC1  (8-17  Mfl)  using  an 
Axoclamp  2A  amplifier  (Molecular  Devices).  Two-electrode  voltage- 
clamp  (TEVC)  and  two -electrode  current-clamp  (TECC)  protocols  were 
created,  driven,  and  recorded  with  Clampex  9.2  software  (Molecular 


Devices).  Recordings  were  analyzed  with  Clampfit  9.2  (Molecular  De¬ 
vices).  R ^  was  measured  in  TECC  with  10  negative  current  injection 
steps  (f  =  2  s,  —  1  to  — 10  nA). 

TEVC  measurements  for  IA,  IHTK,  /Ca,  and  fNaP  were  made  using  pro¬ 
tocols  described  in  detail  in  our  previous  work  (Ransdell  et  al.,  2012, 
2013).  Outward  currents  were  measured  at  0  mV,  inward  currents  at 
peak  conductance. 

Stimulus  protocols.  To  generate  current  injection  waveforms,  we  re¬ 
corded  voltage  waveforms  from  LC  somata  during  intact  network  activ¬ 
ity.  For  Stimulus  Protocol  1,  this  consisted  of  20  s  recordings  from  LC3 
that  included  four  burst  potentials  with  both  pacemaker  EPSPs  and  LC 
back-propagating  APs  (BPAPs).  During  this  recording  Vrest  =  — 55  mV. 
For  Stimulus  Protocol  2,  to  remove  LC  BPAPs,  we  hyperpolarized  the  LC 
soma  to  —65  mV  to  eliminate  spiking.  This  enabled  us  to  record  the 
voltage  change  caused  by  a  single  burst  of  pacemaker  EPSPs.  We  then 
isolated  LC3  with  ligature  and  used  the  voltage  waveforms  recorded  as 
the  command  voltage  in  TEVC.  The  average  current  injection  over  five 
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Figure  2.  Realistic  current  injection  protocols  used  to  test  intrinsic  excitability  of  isolated  LCs.  A,  Stimulus  protocol  1  reproduces  both  synaptic  currents  and  back  propagating  APs  (BPAPs).  LCs  3, 
4,  and  5  are  isolated  from  the  network  via  thread  ligatures  (slash  marks),  and  the  stimulus  protocol  is  applied  before  (Al,  control,  red  box)  and  after  (A2-TEA  exposed,  blue  box)  exposure  to  25  mM 
TEA.  Representative  traces  show  the  response  to  the  entire  stimulus  protocol,  and  overlaid  traces  show  the  fourth  burst  in  the  sequence.  Scatterplots  show  pairwise  correlation  of  time  matched 
voltages  (sampled  at  500  Hz)  of  the  waveforms  shown  in  the  representative  traces.  R  values  are  results  of  Pearson's  correlation  test  for  these  two  cells.  Loss  of  conserved  output  is  demonstrated  by 
the  lack  of  correlation  in  the  scatterplot  as  well  as  the  decrease  in  R  value.  B,  Stimulus  Protocol  2  reproduces  synaptic  currents  LCs  receive  from  pacemaker  neurons.  Experiments  and  analyses  were 
performed  as  in  Stimulus  Protocol  1. 


sweeps  necessary  to  generate  each  of  the  voltage  waveforms  was  used  for 
Stimulus  Protocols  1  and  2. 

Each  stimulus  protocol  was  applied  at  Vrest  and  four  starting  voltages 
(  —  60,  —50,  —40,  — 30mV)  controlled  with  DC  current  injection.  Block¬ 
ers  (25  mM  TEA  or  1  mM  4AP)  were  applied  to  the  preparation  for  30  min, 
after  which  stimulus  protocols  were  repeated  at  each  starting  voltage. 

Statistics.  Analyses  used  SigmaPlot  11.0  and  Clampfit  9.2  software.  R 
values  are  results  of  Pearson  correlation  (Fig.  2)  or  cross-correlation 
function  estimate  (Figs.  3,  4)  with  a  10  ms  lag;  in  cross-correlation  anal¬ 
ysis,  the  peak  R  value  across  the  lag  period  was  used.  Correlation  analyses 
were  conducted  only  between  LCs  of  the  same  ganglion  using  the  entire 
voltage  trace  from  each  stimulus  protocol.  After  failing  Shapiro-Wilk 
tests  for  normality,  changes  in  Rin,  VTtxt,  and  differences  between  R  values 


before  and  after  pharmacological  blockade  were  analyzed  with  Mann- 
Whitney  U  tests. 

Results 

In  the  intact  cardiac  ganglion  (Fig.  1A),  bursts  in  LCs  are  initiated 
via  excitatory  input  from  pacemaker  interneurons  (Fig.  IB)  that 
result  in  synchronous  burst  and  AP  generation  in  all  five  LC 
motor  neurons  (Fig.  IB).  In  addition,  LC  somata  contain  the 
intrinsic  conductances  that  influence  burst  potentials  (Fig.  1C). 
However,  the  magnitudes  of  these  currents  vary  not  only  across 
animals  (Fig.  1C),  but  even  among  cells  ofthe  same  ganglion  (Fig. 
ID).  If  variability  in  current  magnitude  suggests  that  individual 
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Figure  3.  Changes  in  output  after  TEA  exposure.  A,  B,  Box  plots  show  distributions  of  R  values  from  cross-correlation  analysis  of  LC  voltage  waveforms  generated  by  Stimulus  Protocol  1  (21)  and 
Stimulus  Protocol  2  (B)  before  (red  boxes)  and  after  (blue  boxes)  25  mM  TEA  exposure.  Lines  within  boxes  mark  the  median,  box  boundaries  represent  25 th  and  75 th  percentiles,  and  whiskers 
represent  maximum  and  minimum  values.  All  comparisons  were  made  between  LCs  of  the  same  individual  (IV  =  12  comparisons  across  6  individuals).  Each  stimulus  was  injected  into  an  isolated 
LC  at  its  natural  l/rest  and  four  starting  voltages.  Significant  differences  in  median  R  value  were  tested  with  Mann-Whitney  U  tests,  and  are  denoted  by  asterisks  (***p  <  0.005,  **p  <  0.01,  *p  < 
0.05).  C,  D,  Comparison  of  change  in  /lvalue  with  difference  in  /?in  between  LCs  after  TEA  exposure  for  Stimulus  Protocol  1  (0  and  Stimulus  Protocol  2  (0).  E,  F,  Comparison  of  change  in  /(value  with 
change  in  /?in  after  TEA  exposure  for  Stimulus  Protocol  1  (£)  and  Stimulus  Protocol  2  (F). 


cells  differentially  balance  their  conductances  to  achieve  a  com¬ 
mon  output,  then  blocking  a  subset  of  conductances  should  re¬ 
veal  differential  tuning  of  the  remainder,  as  demonstrated  by  a 
shift  away  from  conserved  output.  Our  initial  data  were  consis¬ 
tent  with  this  hypothesis:  when  LCs  were  exposed  to  high- 
threshold  K+  current  UHtk)  blocker  (TEA;  Fig.  ID,  inset),  their 
normally  synchronous  waveforms  in  spontaneously  active  net¬ 
works  became  desynchronized  and  variable  (Fig.  IE).  TEA  appli¬ 
cation  localized  to  LC  somata  is  able  to  affect  spiking  output, 
causing  changes  in  the  spiking  pattern  as  well  as  apparent  desyn¬ 
chronization  of  spiking  as  revealed  by  a  decrease  in  spike  ampli¬ 
tude  that  we  interpret  as  loss  of  summed,  synchronous  AP 
activity  (Fig.  IE ) .  This  decrease  in  amplitude  cannot  be  attributed 
to  TEA  block  of  APs,  as  the  somata  localization  of  the  TEA  in  this 


experiment  left  axonal  conductances  intact.  These  data  demon¬ 
strate  that  LC  somatic  conductances  play  a  significant  role  in 
shaping  waveform  as  well  as  spiking  output  in  these  cells. 

Because  cells  in  intact  networks  are  subject  to  many  factors 
that  influence  output,  we  used  an  experimental  strategy  to  test 
solely  the  influence  of  variable  conductances  on  output  in  neu¬ 
rons  isolated  from  their  constituent  network  connections.  To 
reliably  demonstrate  a  shift  away  from  conserved  output  at  the 
single  cell  level  following  conductance  block,  we  devised  stan¬ 
dardized,  biologically  relevant  current  injection  protocols.  We 
isolated  the  somata  of  the  three  anterior  LCs  from  pacemaker 
synaptic  input,  electrical  coupling,  and  BPAPs,  effectively  isolat¬ 
ing  their  burst  generating  conductances  (Fig.  2 A,B).  We  then 
used  current  injection  to  characterize  cell  output  before  and  after 
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Figure  4.  Effects  of  4-AP  on  conserved  output.  A,  Overlaid  representative  traces  of  LC  responses  to  Stimulus  Protocol  1  before 
and  after  4-AP  exposure.  R  values  from  cross-correlation  analyses  as  in  Figure  3.  B,  Box  plots  show  distributions  of  the  /lvalues  from 
cross-correlation  analysis  of  isolated  LCs  voltage  waveforms  generated  by  Stimulus  Protocols  1  and  2  before  (red  boxes)  and  after 
(blue  boxes)  1  mM  4AP  exposure.  The  line  within  the  box  marks  the  median,  box  boundaries  represent  the  25 th  and  75 th  percen¬ 
tiles,  whiskers  above  and  below  represent  the  maximum  and  minimum  values.  All  comparisons  made  between  LCs  of  the  same 
individual  (N  =  7  and  N  =  6  for  Stimulus  Protocols  1  and  2,  respectively),  as  described  in  Figure  3.  Differences  in  median  R  value 
were  tested  with  Mann-Whitney  U  tests,  and  are  denoted  by  asterisks  (***p  <  0.005,  **p  <  0.01,  *p  <  0.05). 


blocking  conductances.  In  designing  the 
protocol,  which  employs  somatic  current 
injections,  we  considered  that  somata  nat¬ 
urally  receive  two  kinds  of  current  input: 
synaptic  and  BPAPs.  We  recorded  these 
somatic  input  currents  from  cells  in  intact 
networks,  and  subsequently  used  these  re¬ 
cordings  as  the  basis  for  realistic  current 
injections  applied  to  isolated  cell  somata 
(Fig.  2 A,B;  see  Materials  and  Methods). 

Stimulus  Protocol  1  consists  of  both  bio¬ 
logically  realistic  synaptic  currents  plus 
BPAPs  (Fig.  2A).  Stimulus  Protocol  2  con¬ 
sists  solely  of  current  injection  of  synaptic 
drive  to  the  soma  via  pacemaker  inputs 
(Fig.  2 B).  We  first  determined  that  our 
current  injection  protocols,  when  applied 
to  isolated  somata  with  all  intrinsic  con¬ 
ductances  intact,  produced  virtually  iden¬ 
tical  voltage  responses  across  LCs  as 
measured  by  cross-correlation  analysis. 

Under  control  conditions,  current  injec¬ 
tion  elicited  conserved  membrane 
waveforms  across  LCs  in  the  same  net¬ 
work  (Fig.  2  Al,  Bly  controls).  Thus  the 
basis  for  our  experiments  is  whether 
these  conserved  output  responses  are 
maintained  subsequent  to  pharmaco¬ 
logical  blockade  of  conductances. 

We  considered  four  key  points  in  test¬ 
ing  the  hypothesis  of  conserved  output 
from  variable  underlying  conductances. 

First,  we  ensured  that  individual  conduc¬ 
tance  levels  do  indeed  vary  across  isolated 
cells.  We  measured  current  levels  carried 
by  four  ionic  conductances  in  LC  somata: 

IHtk  (high-threshold  K  +  current,  a  com¬ 
bination  of  7Kd  and  LKca)>  I  a  (transient 
K  +  ),  7Ca,  and  7NaP  (persistent  Na+).  Each 
current  measured  varied  at  least  2- 4-fold 
across  individual  cells  (Fig.  1C).  Second,  if 
output  simply  were  constrained  by  elec¬ 
trical  synapses  between  LCs,  then  isolated 
neurons  should  no  longer  produce  con¬ 
served  output  from  standardized  current 
injection.  This  was  not  the  case;  isolated 
LCs  with  all  conductances  intact  pro¬ 
duced  remarkably  consistent  output  (Fig. 

2,  controls).  Third,  if  conserved  output  is 
simply  a  result  of  isolated  cells  passively 
following  current  injection,  then  blocking  conductances  should 
not  alter  this  output.  This  also  was  not  the  case;  when  7HXK  is 
blocked  with  TEA  (Fig.  ID),  cell  output  as  a  result  of  the  stan¬ 
dardized  current  injection  was  substantially  altered  (Fig.  2  A2,  B2, 
blue  boxes).  Finally,  and  perhaps  most  pertinent  to  the  hypothe¬ 
sis,  if  cells  are  not  differentially  tuned  to  produce  similar  output, 
then  blocking  a  conductance  might  still  be  expected  to  alter  out¬ 
put,  but  should  produce  uniform  changes  across  cells.  This  was 
not  the  case;  when  TEA  was  applied  to  isolated  LCs,  and  stan¬ 
dardized  current  injection  applied,  the  once  similar  outputs  di¬ 
verge  between  cells  (Fig.  2)  as  revealed  by  significantly  lower  (and 
highly  variable)  R  values  from  cross-correlation  analyses  follow¬ 
ing  TEA  exposure  (Fig.  3  A,B).  To  ensure  that  these  differences 


were  not  due  to  changes  in  Vresl  after  TEA  exposure,  current 
injections  were  performed  from  native  Vrest,  as  well  as  standard¬ 
ized  baseline  voltages  ranging  from  —  60  mV  to  —  30  mV.  TEA  did 
not  cause  a  significant  change  in  Vrest  (mean  ±SD,  control: 
-49  ±  7.03  mV,  TEA:  -46.3  ±  10  mV,  p  =  0.23,  N  =  12). 
However,  TEA  did  cause  significant  reduction  in  mean  R  value 
between  cells  at  Vresti  —50  mV,  —40  mV,  and  —30  mV  starting 
voltages,  but  not  at  —60  mV,  in  both  Stimulus  Protocols  (Fig. 
3 A,B).  At  —60  mV,  the  depolarization  in  membrane  potential 
due  to  current  injection  did  not  reach  the  activation  voltages  of 
most  of  the  conductances  (Ransdell  et  al.,  2012, 2013).  Therefore, 
the  —  60  mV  data  actually  demonstrate  what  happens  when  these 
cells  are  simply  passive  followers  of  current  injection,  and  sup- 
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port  the  result  that  at  voltages  where  conductances  are  activated, 
there  is  a  true  loss  of  tuning  when  a  subset  of  currents  are  blocked. 

To  ensure  these  results  are  not  due  to  changes  in  input  resis¬ 
tance  ( Rin )  of  cells  in  TEA,  we  analyzed  relationships  between  Rin 
and  activity.  TEA  did  not  significantly  or  predictably  alter  Rin 
(mean  ±SD,  control:  5.03  ±  2.4  Mil,  TEA:  5.25  ±  3.28  Mfl,p  = 
0.795,  N  =  12).  However,  we  predicted  that  if  loss  of  conserved 
activity  across  LCs  with  blockade  were  a  result  of  variable  impact 
of  the  TEA  on  Rin,  then  the  largest  changes  in  activity  between 
cells  should  occur  with  the  largest  difference  in  Rin  between  those 
two  cells.  We  found  no  such  relationship  between  change  in  Rin  of 
individual  cells,  or  the  difference  in  Rin  between  pairs  of  cells, 
with  the  change  in  R  value  in  the  corresponding  cross-correlation 
(Fig.  3C-F). 

Finally,  if  this  is  a  general  phenomenon  and  not  just  specific  to 
fHXK,  then  blocking  another  conductance  in  this  system  also 
should  result  in  loss  of  conserved  output.  We  blocked  a  second 
K+  current,  IA,  using  4-aminopyridine  (4-AP).  4-AP  caused  the 
same  loss  of  conserved  output  in  LCs  (Fig.  4A)  with  both  Stimu¬ 
lus  Protocols,  as  well  as  an  overall  increase  in  variability  of  R 
values  (Fig.  4B). 

Discussion 

Evidence  for  compensation  of  intrinsic  excitability  that  stabilizes 
neuronal  output  has  accumulated  from  nervous  systems  across 
phyla,  ranging  from  the  cardiac  neurons  used  in  this  study  (Rans¬ 
dell  et  al.,  2012)  to  mammalian  systems  (Desai  et  al.,  1999;  Ner- 
bonne  et  al.,  2008).  This  homeostatic  plasticity  of  intrinsic 
excitability  is  likely  a  critical  mechanism  underlying  neural  net¬ 
work  dynamics  (Davis,  2006;  Turrigiano,  2011),  and  maybe  in¬ 
volved  in  disorders  of  neuronal  excitability,  such  as  epilepsy 
(Echegoyen  et  al.,  2007;  Howard  et  al.,  2007).  These  compensa¬ 
tory  responses,  and  the  changes  and  variability  in  ionic  conduc¬ 
tances  they  induce,  are  usually  the  result  of  pathology,  injury,  or 
experimental  manipulations  that  substantially  and  drastically  al¬ 
ter  the  parameter  landscape  of  these  neurons.  Yet  few  of  these 
model  systems  capture  native  variability  of  neurons  within  un¬ 
manipulated  networks.  So  the  question  remains,  do  neurons  in 
native  networks,  that  are  genetically  identical  (i.e.,  cells  within  the 
same  network),  achieve  similar  output  at  least  in  part  by  differ¬ 
entially  tuning  the  underlying  intrinsic  conductances  via  distinct 
mechanisms?  It  is  this  question,  independent  of  compensation, 
that  our  study  addresses. 

We  addressed  several  key  components  of  this  question  in  this 
study.  First,  we  demonstrated  that  motor  neurons  of  the  cardiac 
ganglion  have  similar  membrane  potential  waveforms  within  the 
same  network,  both  during  intact  network  activity  (Fig.  IB)  and 
as  isolated  neurons  responding  to  biologically  relevant  current 
injection  protocols  (Fig.  2).  This  ability  to  record  from  unambig¬ 
uously  identified  neurons  both  in  intact  networks  and  under  ex¬ 
perimentally  controlled  conditions,  both  of  which  are  crucial  to 
directly  test  the  hypothesis  in  this  paper,  represents  a  key  strength 
of  this  model  system,  and  how  it  can  contribute  uniquely  to  this 
question.  Second,  we  measured  native  current  magnitudes  of 
multiple  ionic  conductances,  to  ensure  that  these  cells  also 
showed  variability  in  magnitudes  of  conductances  as  reported  in 
other  systems  (Fig.  1C).  Not  only  did  these  cells  show  2-4-fold 
variability  in  their  underlying  conductances,  but  even  cells  within 
the  same  network  showed  substantial  variability  in  current  mag¬ 
nitudes  (Fig.  ID).  Third,  we  ruled  out  the  possibility  that  the 
voltage  waveforms  of  isolated  cells  were  passively  following  ex¬ 
ternal  current  injection  (Fig.  2 C,D).  Finally,  these  preconditions 
allowed  us  to  directly  pose  the  experimental  question:  does  nat¬ 


ural  variability  in  underlying  current  magnitudes  result  from  dif¬ 
ferentially  tuned  conductances  converging  on  similar  output  via 
distinct  solutions?  We  tested  this  with  a  relatively  straightforward 
experimental  design.  Namely,  if  output  is  the  result  of  finely  bal¬ 
anced  ratios  of  ionic  conductances,  then  disrupting  one  conduc¬ 
tance  should  cause  each  cell  to  experience  a  cell-specific 
detuning,  resulting  in  loss  of  conserved  output.  Our  results  were 
entirely  consistent  with  the  hypothesis:  blocking  either  /HTK  or 
independently  results  in  loss  of  similar  output  across  three  LC 
neurons  within  the  same  network  (Figs.  2,  3  A,B). 

Similar  manipulations  in  pacemaker  neurons  of  the  stomato- 
gastric  ganglion  (STG)  have  yielded  different  results  (Sziics  and 
Selverston,  2006).  In  the  STG,  4AP  used  to  disrupt  A- type  current 
in  paired  pyloric  dilator  neurons,  also  reported  to  have  variable 
conductance  magnitudes  (MacLean  et  al.,  2005),  showed  little 
impact  on  synchronized  activity  between  these  cells  (Sziics  and 
Selverston,  2006).  These  data,  although  seemingly  contradictory, 
are  intriguing  as  Sziics  and  Selverston  (2006)  left  chemical  and 
electrical  synaptic  connectivity,  as  well  as  neuromodulatory  in¬ 
puts  to  these  networks  (Marder  and  Thirumalai,  2002)  intact 
during  4AP  application.  We  carefully  selected  cardiac  LCs  as  a 
model  where  we  have  control  over  these  higher-order  network 
properties,  allowing  us  to  focus  solely  on  intrinsic  excitability. 
Interpreted  in  the  context  of  the  Sziics  and  Selverston  (2006) 
work,  convergent  output  as  a  result  of  differentially  tuned  con¬ 
ductances  likely  is  only  one  layer  of  regulation  in  these  systems 
ensuring  synchronized  output  (Prinz  et  al.,  2004).  Work  is  ongo¬ 
ing  to  determine  whether  and  how  higher-order  processes,  such 
as  electrical  coupling  and  neuromodulation,  ensure  that  cellular 
and  network  output  are  maintained  in  patterns  necessary  for 
robust  function  (Grashow  et  al.,  2009,  2010;  Marder,  2012). 

The  most  parsimonious  explanation  for  our  results  is  that  cells 
with  variable  underlying  conductances  are  indeed  individually 
tuned  to  produce  conserved  output.  We  believe  these  findings 
provide  a  clear  and  unambiguous  demonstration  of  this  phe¬ 
nomenon,  helping  unify  a  substantial  amount  of  compelling  and 
consistent  evidence  in  multiple  model  systems,  both  biological 
and  theoretical  (Marder  and  Taylor,  2011;  Marder,  2011;  Turri¬ 
giano,  2011).  It  is  striking  that  differential  tuning  is  not  simply 
found  across  different  individuals,  but  is  seen  even  among  cells 
with  in  the  same  individual.  One  compelling  question  that  arises 
from  these  findings  is  the  nature  of  the  developmental  and  post- 
developmental  origin  of  this  variability.  Specifically,  are  individ¬ 
ual  neurons  produced  with  common  underlying  conductance 
relationships  that  then  diverge  over  time  (via  compensation),  or 
is  variability  inherent  to  the  initial  development  and  tuning  of  a 
neuron,  and  maintained  across  the  lifetime  of  the  animal?  An¬ 
swers  to  these  questions  will  rely  on  model  systems  that  allow  for 
identified  neurons  to  be  characterized  effectively  over  both  devel¬ 
opmental  and  postdevelopmental  time  scales. 
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